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PREFACE 


The  Department  of  Aeronautics ,  Naval  Postgraduate  School 
offers  several  courses  at  differ*  't  levels  on  the  dynamics  of 
flight  vehicles.  These  courses  include  studies  of  fixed  wing 
stability,  rotary  wing  stability,  and  the  stability  of 
missiles.  In  addition,  there  are  general  courses  that  cover 
the  dynamics  of  bodies  in  six  degrees  of  freedom,  including 
the  development  of  the  equations  of  motion. 

There  is  currently  no  signle  resource  available  that 
covers  the  material  for  the  airplane  and  missile  courses,  as 
well  as  the  introductory  dynamics  classes,  and,  as  a  result, 
the  instructors  must  make  use  of  large  numbers  of  locally 
prepared  handouts  and/or  assign  several  textbooks,  only 
portions  of  which  will  be  used. 

In  order  to  provide  a  solution  to  this  problem,  the 
following  material  has  been  developed  for  use  in  the  dynamics 
classes  in  the  Department  of  Aeronautics.  Although  this 
material  is,  in  places,  of  greater  depth  than  what  would  be 
covered  in  the  courses,  it  is  felt  that  the  value  of  this 
material  is  enhanced  by  its  degree  of  thoroughness. 
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CHAPTER  1 


DYNAMICS  OF  MOTION 


1.1  INTRODUCTION 

This  chapter  is  intended  to  be  a  formal  review  of 
mathematical  and  physical  concepts  that  will  be  used  to 
develop  the  theory  of  dynamic  stability  of  ^flight  vehicles. 

The  starting  point  is  the  basic  laws  of  Newtonian 
Mechanics,  with  a  useful  application  to  spring-mass-damper 
systems.  We  shall  find  later  how  the  motion  of  any 
conventional  flight  vehicle  is  similar  to  those  systems. 

We  then  proceed  to  rigid  body  equations  of  motion, 
described  in  non- inertial  reference  frames. 

1.2  NEWTON'S  LAWS  OF  MOTION 

In  any  system  of  units,  basic  quantities  and  units  are 
related  through  Newton's  laws  of  motion. 

NEWTON'S  1ST  LAW  OR  THE  LAW  OF  COORDINATE  SYSTEMS: 

If  there  is  no  force  acting  on  a  particle,  its  velocity 
remains  constant. 

This  law  is  often  called  the  law  of  inertia  and  applies 
to  inertial  frames  of  reference. 

NEWTON'S  2ND  LAW  OR  THE  LAW  OF  DYNAMICS: 

The  time  rate  of  change  of  the  linear  momentum  of  a 
particle,  is  proportional  to  the  net  force  acting  on  it. 


The  vector  quantity 


■pff)  -  'Wl  (i)  v(i) 


(1-1) 


is  defined  as  the  linear  momentum  of  the  particle,  where  m(t) 
is  the  inertial  mass  of  the  particle  and  v(t)  is  its  velocity 
relative  to  the  reference  frame  in  question. 

The  2nd  law  is  mathematically  expressed  as 


£  l  dpV 

hn,i  =  &  cA 


(1-2) 


F,H  S(t)J 


If  the  mass  m  remains  constant,  we  can  write: 


(1-3) 


*  l  *  £  Sm 


(1-4) 


3  T  5  (i] 


(1-5) 


where  a(t)  is  the  acceleration  caused  by  the  net  applied 
force  Fn«t. 

The  proportionality  factor  is  expressed  as  Vgc  and  its 
value  depends  on  the  system  of  units  being  used. 

In  the  English  Engineering  system,  one  pound  force 
(1-lbf)  will  cause  to  one  pound  mass  (1-lbm) ,  an  acceleration 
of  32.174  ft/sec*  Thus  substituting  into  Equation  1-5 


l-lbf  *  £  flbm)  32  174  (M/sec*) 


(1-6) 


Hence 


gc  =  32  174  Ibm  li/lbf  see2  (1-7) 

A  few  other  units  used  for  velocity,  are  given  below  with 
their  corresponding  conversion  factors: 

1  mile  per  hour  (lmph)  =  1.4666  ft/sec 
1  knot  =  1  nautical  mile  per  hour  =  1.6878  ft/sec 

NEWTON'S  3RD  LAW  OR  THE  LAW  OF  MUTUAL  FORCES: 

If  one  body  exerts  a  force  on  a  second  body,  the  latter 
exerts  an  equal  but  opposite  force  on  the  first. 

This  law  allows  us  to  define  the  gravity  force,  or  the 
weight  of  a  body,  that  has  a  given  mass  (m) .  It  is  the  force 
that  results  due  to  the  earth's  attraction. 

Besides  the  three  fundamental  laws  of  mechanics,  the 
Newton's  Universal  Gravitational  Law  states  that  two  inertial 
masses  ml  and  m2,  separated  by  a  distance  R  are  subject  to  an 
attraction  force  according  to  the  equation: 

r  n  n  'm>  ^ 

h  '  9C  U  K2  (1-8) 

where  G,  is  the  gravitational  constant,  which  in  the  EE 
system  is  G  =  6.52*40'"  Ibf --fiVlt-m2 
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If  now  m L  =  Me  =  1.32  X  10  lbm  (mass  of  the  earth), 
then  near  earth's  surface  where  K  *  Rc  =  2.09  X  107  ft 
(radius  of  the  earth),  we  can  write: 


or 

W  (1-10) 

where  g  *  gtGMe/R \  =*  32  .174  ft/sec2,  called  acceleration  of 
gravity. 

By  Equation  1-10  we  can  see  that  1-lbm,  subject  to  the 
earth's  gravitational  field,  experiences  a  force  1-lbf  and 
hence  it  accelerates  with  32.174  ft/sec2,  near  earth's 
surface. 

We  can  also  define  another  unit  of  mass,  the  slug,  as  the 
mass  that  experiences  an  acceleration  1-ft/sec2,  if  subject 
to  a  force  of  1-lbf. 

By  this  definition 

1  slug  *  32.174  lbm  (1-11) 

Consequently,  we  can  derive  the  density  (p)  ,  as  the  mass 
per  unit  volume  usually  measured  in  slugs/ft. 

dens,l!l  =  'voUmeTv)  (1-12) 

People  in  the  field  use  the  (g) ,  as  a  short  hand  notation 
for  force  or  acceleration.  A  5-g  force  simply  means  a  force 
5  times  the  body's  weight. 

1.3  ANGULAR  MOMENTUM 

A  special  significance  to  the  2nd  law  is  given  by  the 
angular  momentum  approach  derived  simply  by  crossing  the 
radius  vector  7  ,  on  to  the  2nd  law. 
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v. 


(1-13) 


We  define 


•*  -*  j- 

r%  =  r  x  F^et 

at 


G-wt  -  rx 


(1-14) 


the  net  torque  or  moment  of  the  force  about  the  origin  and 

h  =  rxp  (1-15) 


the  angular  momentum  about  the  origin. 


Differentiating  Equation  1-15  with  respect  to  time,  we 


db  dr  -*  -*  dp 

ST  =dTxP’rx3F 

The  first  term  is  zero,  because 


(1-16) 


—  x-p  =  'Jx(wd)  =  'vn(.uxu)*»o 
dt  1 


(1-17) 


Figure  1-1 

Moment  of  Force  About  the  Origin 


So  we  finally  get 


*  G  M<»t 


(1-18) 
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which  tells  us  that  the  time  rate  of  change  of  the  angular 
momentum  of  a  particle,  equals  to  the  net  torque  or  moment 
acting  on  it. 

If  no  torque  is  acting  on  a  particle,  its  angular 
momentum  is  conserved,  likewise  its  linear  momentum  is 
conserved,  if  no  force  is  acting  on  it. 

A  simple  example 

For  a  pendulum  that  swings  in  a  vertical  plane,  having  a 
mass  (m) ,  the  acting  forces  are,  its  weight  W  *  mg  and  the 
tension  of  the  rod  T.  The  rod  is  assumed  to  have  negligible 
mass. 


Figure  1-2 
A  Simple  Pendulum 

If  9  is  the  angle  between  the  rod  and  the  vertical  axis, 
s,  is  the  arc  length  corresponding  to  the  angle  0  and  1,  is 
the  length  of  the  rod,  a  direct  application  to  the  2nd  law 
(Equation  1-5),  gives: 
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w  s 


(1-19) 


-  Trig  sin©  = 

and  equilibrium  condition  in  f  be  direction  of  the  rod,  gives: 


-  Tncj  cos£  =  T 

(1-20) 

Equation  1-19, 

can  be  written 

s  +  g  sm  9=0 

(1-21) 

If  @  is  small. 

we  can  assume 

s  =  29 


(1-22) 


and 


©  S  sin©  (1-23) 

With  this  assumption.  Equation  1-21,  can  be  written 


Q  0  =0 


(1-24) 


Equation  1-24,  is  a  linear  differential  equation,  with 
constant  coefficients,  which  describes  the  motion  of  the 
pendulum. 

The  same  result  can  be  derived,  if  we  use  the  angular 
momentum  approach. 

The  angular  momentum  of  the  pendulum,  about  the  pivot,  is 
simply  (2),  cross  the  linear  momentum  ('ms).  The  time  rate  of 
change  of  the  angular  momentum,  has  to  be  equal  to  the  net 
external  torque. 


—  (?x<ms)  *  -2  mil  sm© 

Since  i  is  constant  (d?/dt  *  0) ,  we  can  write 

im  s  -  -  2  vng.  sm  0 


(1-25) 


(1-26) 
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The  minus  sign  in  Equations  1-19  and  1-25 ,  indicates  that 
if  the  angle  0  decreases,  the  acceleration  s  in  Equation  1-19 
and  the  angular  momentum  in  Equation  1-25,  increase. 

1.4  UNDAMPED  HARMONIC  OSCILLATOR 

In  this  section  and  in  the  following  section,  we  will 
derive  the  equations  and  the  solutions  to  the  spring  mass 
damper  systems. 

The  undamped  oscillator  is  simply  a  spring  mass  system 
without  friction,  shown  in  Figure  1-3. 

Let  k,  denote  the  spring  stiffness,  m,  the  mass  of  the 
body,  x(t),  the  displacement  of  the  body  from  the  equilibrium 
position.  Notice  that  the  weight  of  the  body  stretches  the 
spring  permanently,  so  that  xe  is  the  equilibrium  position, 
after  the  body  has  hanged  to  the  spring. 


Figure  1-3 
Spring  Mass  System 


(1-29) 


The  restoring  force  according  to  Hook's  law,  is 

Fs  =  -  k  x 

Then  according  to  the  2nd  law 


m  x  =  -  k  x 

(1-30) 

X  +  —  X  s  0 
m 

(1-31) 

As  an  example  we  will  take 

m  =  18  lbm 

(1-32) 

k  ■  40.5  lbf/ft 

(1-33) 

If  we  let 


Mn 


(1-34) 


where  is  the  so  called  natural  frequency  of  oscillations. 
Equation  1-31  becomes 

*  +  x  =0  (1-35) 

This  is  a  linear  differential  equation  with  constant 
coefficients.  Its  characteristic  equation  is  clearly 

3*+utf  -  0  d-36) 

which  has  characteristic  roots 

(1-37) 

The  most  general  solution  of  Equation  1-3-5  can  be 
expressed  as 


xttMM 


JUJnt 

e 


(1-38) 
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where  A,.  and  A_  are  constants  depending  on  the  initial 
conditions. 

Expanding  the  exponential  terras  of  Equation  1-38,  we 
obtain 


x(-t)  =  a  COS  (Uni  +P  smuJnt 

(1-39) 

where 

<x  =  A  +  +  A_ 

(1-40) 

and 

fi.j(A.-A.) 

(1-41) 

If 

we  now  express 

smuj„i  =  cos  (uint  *  90°) 

(1-42) 

we  can 

give  the  following  useful  form  to  our 

solution 

x(i)  *  X*,ax  COS  ( UJ,,i  4-0.) 

(1-43) 

where  is  the  maximum  deviation  from  the  equilibrium 

position  and  B.  ,  the  initial  phase. 

We  can  relate  the  quantities  xWBX  and  ©.  ,  with  the 
initial  conditions,  as  follows: 

At  time  t=0,  the  initial  displacement  x.  ,  is  given  by 
Equation  1-43. 

x„  =  xWQx  cos  0O  (1-44) 

and  the  initial  velocity  x„  ,  is  given  by  differentiating 
Equation  1-43  with  respect  to  time  and  evaluating  at  t=0. 

x.  -  -  sin  ©„  (1—45) 

As  an  example  we  will  take 


*«.  »  4  ft 
*<>  *  32  ft/sec 


(1-46) 

(1-47) 


From  Equations  1-44  and  1-45,  the  initial  phase  and  the 
maximum  displacement  can  be  found  to  be 


0o  =  4an ' ( 


-uax„ 


-) 


\'/2 


X  -  +  /  xz  +  ) 

*ma.%  ~  -  l  (A?  / 


(1-48) 


(1-49) 


Note:  In  Equation  1-48,  keep  the  sign,  so  that  the 
quantity  in  brackets  will  give  the  right  quadrant  for  the 
angle  0* . 

For  our  example 


,  40  5 

t^n  *  gc  -  »  32.174  >  =  79  39  sec'2 


m 
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(1-50) 


0„  =  W(-~~)  =  456  76°  =  2.59 
»  *  (<6  * =  ±  5.49  « 


(1-51) 


(1-52) 


so  that  the  response  is 

x (1 )  -  -5.49  cos  (  8  3i  i  +  2.59  rad)  (1-53) 

1.5  DAMPED  HARMONIC  OSCILLATOR 

Consider  the  spring  mass  damper  system  shown  in 
Figure  1-4.  Assume  that  the  friction  force  exerted  by  the 
damper,  is  a  linear  function  of  the  mass  velocity. 


1 


*•_ 
xft)  l 


///////// 

s 


k 


I 


7/77777? 


Figure  1-4 

Spring  Mass  Damper  System 

Let  k  the  spring  stiffness,  m  the  mass  of  the  body,  x(t) 
the  displacement  of  the  body  from  the  equilibrium  position 
xe  ,  as  it  was  defined  in  section  (1.4) ,  and  c,  the 
coefficient  of  friction. 

The  restoring  force  is  again 

Fs  =  -k  x  (1-54) 

and  the  friction  force  is 

Tf  =  '  cx 

Then  according  to  the  2nd  law 

ofnx  =  -  c  x  ~  k  x 

or  c 

x  +•  —  x  +  —  x  =  0 

'in  m 


I 


(1-55) 

(1-56) 
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As  an  example  we  will  take 
m  *  18  lbm 


(1-58) 


k  *  40.5  lbf/ft 


(1-59) 


If  we  let 


(1-60) 


=  3*  * 


(1-61) 


where  5  is  the  so  called  damping  coefficient,  then  the 


Equation  1-57  is  written 


v  +  2  5u)„  x  +•  (Ufl  x  =  O 


(1-62) 


This  is  a  linear  differential  equation  with  constant 
coefficients.  Its  characteristic  equation  is  clearly 

32  +  25Wt,3  ♦  «*,*  =0  (1-63) 

Three  cases  arise  in  obtaining  solutions  to  the 
differential  Equation  1-62. 


As  an  example,  we  will  take 
■  4  ft 

*.  =  32  ft/sec 


As  an  example  we  will  take 


(1-64) 

(1-65) 


so  that 


c  =  9.52  Ibf/^.sec'1 


^  =  ^  2m<Vr,  =  1 


(1-66) 


(1-67) 


In  this  case,  the  characteristic  equation  has  a  double 


root 


=  -  tOr 


(1-68) 


I 


(1-69) 


and  the  most  general  solution  is 

xU)  =  (A*  +  A.i)  €_a,wt 

where  A  +  and  A_  are  constants  determined  by  the  initial 
conditions  as  follows 

A,  =  *o  (1-70) 

and 

A.  =  *„  *  x0u)„  (1-71) 

so  that  the  solution  in  terms  of  the  initial  conditions  can 
be  written 

v  r  -  -U)nt 

x  (4)  =  [xo((fa)n±)  +  x.  ij  e  (1-72) 

For  our  example 

,  .  r  -  8.51  i 

x  U)  r  [  4  (  If  8.51 1)  +  32  t]  G  (1-73) 

or 

,  ,  -8.511 

x(i)  -  (4  >  66.04  i)  C  (1-74) 

Case  2; _ Over  damped  notion.  (5>1) 

As  an  example  we  will  take 

c  =  14.28  lbf/ft-sec'1  (1-75) 

so  that 

5  =  &  *2^7  =  4  5  d-76) 

In  this  case  the  roots  of  the  characteristic  equation  are 
real  and  negative 

3  -  -Bu) v  ±  (Vr,  (5*~l)  /z  (1-77) 

and  the  most  general  solution  is 

-5u)nl  r  ,  -b/StT*  _a)n /5Vi  t  i 
xU)  =  e  [^4.2  +A.  c?  J  (1-78) 
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where  A+  and  A_  are  determined  by  the  initial  conditions  as 
follows 


x0  =  A+  +  A. 

(1-79) 

=  -5tun  (A+  +  A.)  4-cu (X-A.) 

(1-80) 

which  yields 

2  /s’T 

(1-81) 

,  Ui'-i -5)x.-3; 

2/5M 

(1-82) 

For  our  example 

A  +  =»  6  .36  ft 

(1-83) 

A _  *  -2.36  ft 

(1-84) 

and 

9.514  -9.54  4  \ 

x(-l)  =  e  (6.56  e  -  2.56  e  J 

(1-85) 

or 

-i  ,  ,  -  3.26 1  -2228  t 

M  =  6.56  <2  -  2.36  e 

(1-86) 

Case  3; _ Underdamped  Motion  (S<1) 

As  an  example  we  will  take 

c  =  4.76  lbf/ft-sec'1  (1-87) 

so  that 

6  =  2vnuJn  =  °'  (1-88) 

In  this  case  the  roots  of  the  characteristic  equation  are 


imaginary  with  negative  real  parts 


If  we  define 
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and 


5ui„x0+  x0  =  -  Xcud  sin 0„ 


(1-100) 


So  the  constants  are  determined  as 


--bn'1  [  5a)nx0  4-  — “-] 


(1-101) 


x  .i[ 


2  n  i/2 


(1-102) 


Note:  In  Equation  1-101,  keep  the  sign,  so  that  the 
quantity  in  brackets  will  give  the  right  quadrant  for  the 
angle  &  . 


For  our  example 


UJd  =  8.5  {  ( 1  -  0.5* )  2  *  7  57  sec' 
©o  *  12102°  *  2.11  rad 


y.  *  ±  7.76  ft 


so  that  the  response  is 


-42C*t 

x(t)  =-7.76  e  COS  (?3?i  +  2.K  rad) 


(1-103) 

(1-104) 

(1-105) 


(1-106) 


The  roots  of  the  characteristic  equation  can  be  plotted 
in  the  complex  plane  as  shown  in  Figure  1-5. 

In  cases  1  and  2,  the  damping  is  so  large  that  no 
oscillations  take  place  and  the  mass  simply  returns  gradually 
to  the  equilibrium  position  x=0.  In  those  cases  the  roots 
are  real  and  negative. 

In  case  3,  the  damping  has  been  reduced,  so  that 
oscillations  about  the  equilibrium  position,  do  take  place, 
the  amplitude  of  which  dies  out  as  t^®. 
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Figure  1-5 

Plot  of  the  Characteristic  Roots  in  the  Complex  Plane 


Figure  1-6  shows  the  responses  in  the  above  three 
described  cases  where  it  is  assumed  that  the  initial 
displacement  is  ,  and  the  initial  velocity  is  *<>=0. 


critically  dauiped  motion  5*1 

ourdnrwptd  motion  5  >1 

underdo**  pari  motion  5<1 


envelope  C  e 


Figure  1-6 

Response  of  Displacement  (x)  of  a  Spring  Mass  Damper  System 


■  -  . 


Notice  that  in  the  case  where  the  damping  5  is  zero,  we 
have  a  spring  mass  system,  without  friction,  described  in 
section  1.4.  In  this  case,  the  roots  are  imaginary  and  is 
the  same  as  ««  (Figure  1-5) . 

It  is  clear  that  the  system's  physical  constants  (in  our 
case  5 , m  ,  c  ,  k  )  ,  uniquely  define  the  roots  of  the 
characteristic  equation.  Those  roots  are  given  the  name: 
normal  modes  because  one  can  determine  from  those  roots,  the 
corresponding  characteristic  motion  (in  our  case,  undamped 
critically  damped,  overdamped,  underdamped) . 

We  have  seen  so  far  how  an  harmonic  damped  oscillator 
responds,  if  initially  disturbed  from  its  equillibrium 
position.  Then  knowing  the  initial  conditions  we  can  always 
find  the  response  of  the  displacement  or  velocity  versus 
time. 

Now  we  are  interested  in  what  happens  if  the  oscillator 
is  driven  by  an  external  input  force  F(t). 

In  this  case  the  solution  will  consist  of  two  parts: 

xft)  -  ft)  «■  *pari  ft)  (1-107) 

where  is  the  homogenous  part  of  the  solution 

discussed  so  far,  i.e.  considering  the  input  force  zero  and 
*part  (i)  is  the  particular  solution  due  to  this  particular 
applied  force,  i.e.  considering  the  initial  conditions  zero. 

To  obtain  particular  solutions  one  can  use  the  method  of 
undetermined  coefficients,  the  method  of  variation  of 
parameters  or  the  Laplace  Transform  method. 
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In  illustrating  the  above  methods  we  are  considering  the 


following  cases: 

(a)  Sinusoidal  Inputs 

The  spring  mass  damper  is  driven  bythe  sinusoidal  force 
F ( t )  -  F0  cos  a?  t  (1-108) 

As  an  example,  we  will  take 

F (t)  =  24  cos  3.5  t  (1-109) 

The  method  of  undetermined  coefficients  assumes  a 
superposition  of  sinusoidal  solutions  of  the  form 

*par*  W  3  COS Ul t  +•  A,  smo)t  (1-110) 

or  equivalently 

xpart  3  X  COS  (<l)l  +  0.)  (1-111) 

where  X  and  9.  are  the  coefficients  to  be  evaluated. 

Substituting  into  the  differential  equation 

x  +  25cu„  vc  +  co*  x  =  F0  coscot  (1-112) 

we  obtain 

(cW-cu)  cos(u'i+0B)  -  sm  (cot  +  ©<>)  *  —■  coscot  (1-113) 

■A 

Expanding  cosine  and  sine  and  collecting  terms,  we  get 
[  (uu«  -go*)  cos  ©o  -  25ufnd3sin0o]  coscot 

-  [  tun?  -  to* )  sin9„  +  25u>noocos0o]  smoot  a  ^  cosuot  (1-114) 

Equating  coefficients 

(<Ut? -w*)  cos  Qo  -  25uJnUJ  sin  0„  =  —■  (1-115) 

and 

((U„2-CD2)  SinOo  +•  25(JOr,U)  C05  ©o  =0  (1-116) 


From  Equation  1-116  we  obtain 


-tan  0O 


-2$  (OnlO 

UJrf-CO2 


(1-117) 

Multiplying  Equation  1-115  by  cos&>.  Equation  1-116  by 
sin  & r  smd  by  adding ,  we  obtain 


C  U)ri2  -  U)1  )  —  COS  Q0 

X 


(1-118) 


The  cos  0d  is  expressed  in  terms  of  tan  as 
cos  0o  ~ 


V lari  Qo  tUn-u)1 


[  l  *  (^)'f 


zVb 


(1-119) 


so  that 


and 


X  • 


h 


[(*-«•)*  +  (2Wu  )2]"2 


(1-120) 


Go  -  -tun 


-i  -25u>r*cu 


tttf.o)*  (1-121) 

For  our  example  with  the  three  different  values  of  5 
considered,  the  responses  are: 

-for  5-1  Xpt.rtU)  =  0.28  COS  (3.5 1  +5.5  rad)  (1-122) 

-for  5=1.5  x  pari  (4)  =0,22  cos  (  351  +  5.31  rad)  (1-123) 

-for  5*0.5  VpuriU)  *  0.56  cos  (351+5.82  rad)  (1-124) 

The  total  responses  correspondingly  are: 


-8  51  i 

x(l)  r  (4+66.04  4)  g 

+  0.28  cos  (  3.5  4  +  5.5  rad) 

xft)=  6.56  c  “‘-US 

+  0.22  cos  (  3.5  t  +  5.31  rad) 


(1-125) 


(1-126) 
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-4  2 6  i 

x Ct)  *  -7.76  e  cos  C?37 -fc  2.H  rad) 

+  0.56  cos  C  5  5  i  *5.32  rad)  (1-127) 

(b)  Singularity-Inputs 

The  Laplace  Transform  method  can  be  very  easily  applied, 
since  we  are  dealing  with  linear  differential  equations 
having  constant  coefficients. 

Appendix  A  gives  a  short  Laplace  Transform  Table  and  some 
Laplace  Transform  Theorems  that  help  to  transform 
differential  equations  to  algebraic. 

In  illustrating  the  above  method,  we  are  considering  unit 
impulse  and  unit  step  input  forces. 

(1)  Onit  Impulse  Response 

The  spring  mass  damper  is  driven  by  the  singularity 


input 


Fft)  =  W 

(1-128) 

The  differential  equation  becomes 

x  *  2  S  uUy,  x  t  u)nz  x  =  S(t) 

(1-129) 

Neglecting  initial  conditions  and 

Laplace 

transforming,  we  obtain: 

Xparl  ts)  =  C  s2  *  25a>ri  S  +  )  =  1 

(1-130) 

or 

(s)  =  i  /  s2*  25ul„S*a?nz 

(1-131) 

For  our  example  with  $*1 


XportC5)  =  */(s*S5i)2 

t .  \  -  8.51  -t 

x  part  M  =  i  « 
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(1-132) 

(1-133) 


(1-134) 
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Xpart  (s)  =  s/s2  +  25uj„s  +  tt>n2 


(1-144) 


For  our  example  with  5*1 


Xpor{(s)  =  Cs+S.Sl)2  =  IT tsT"  Cs! 


85) 

i+8502 


-85U  -SSI  t 

Xpari  U)  *  e  -  8.51  4  e 

and  the  total  response  is 

x(4)  =  (4+66.044.)  e  +  e  -  i  a 

With  5*1.5 


(1-145) 

(1-146) 


(1-147) 


Xpeirt  (s) 


s  _  017  -1.17 

( s+52t)Cs^22  28)  *  S+526  S+22  2& 


Xo„i  H)  «  0.17  e'5'2Gt  -  117  2 


-22.284 


and  the  total  response  is 

ft)  tit  *i26i  ozC  -22.281  -3.2C1  -22.Z8  4 

x(U  =  6.56e  -2.56  e  +  0.17  e  -i(7c 

With  5  s0 .5 


(1-148) 

(1-149) 


(1-150) 


-10M’  30.ll® 

0.58  2  0582 


y  ,  x _ s _ 05 og 

f0r{  *  “  (s+4  26+j  7  57  )(s+ 4  26 -j  7.57)  *  S+ 4.26 +j  7*7  *"  S+4.26 -j  7.37 


xMri  tt)  =  1  16  e4  26i  cos  (7.57  4  0.55  rod) 


and  the  total  response  is 


,, .  -4.26  4 

x0.)  =  -  7.76  e  Cos  ( 7,57  4  +  2  11  rad )  + 


(1-151) 

(1-152) 


-4.26  4 

+ 1.16  e  cos  (7.57  4  +0.55  rad  ) 


(1-153) 


1 . 6  AXES  SYSTEMS 


In  describing  the  motion  of  our  flight  vehicle,  we  will 
make  the  following  two  assumptions: 

First,  we  will  establish  am  inertial  frame  of  reference 
on  the  earth  considering  the  earth  to  be  fixed  in  space.  In 
this  inertial  frame,  Newton's  laws  are  valid. 

Second,  the  vehicle  will  be  assumed  to  benave  like  a 
rigid  body  so  that  translational  motion  can  be  described  by 
considering  it  as  a  particle  with  all  its  mass  located  at  the 
center  of  mass,  and  rotational  motion  can  be  described  by 
considering  moments  about  that  center  of  mass  and  constant 
moments  of  inertia. 

To  describe  the  translational  and  the  rotational  motion, 
one  has  to  establish  a  second  non-inertial  axes  system,  fixed 
on  the  vehicle  with  origin  at  its  center  of  mass,  in  this 
reference  frame  Newton's  laws  do  not  hold,  in  general,  and  we 
have  to  add  additional  terms  to  our  basic  fundamental 
equations. 

It  is  convenient  to  use  cartesian  coordinate  systems  for 
both  reference  frames  used. 

Let's  consider  the  translational  motion  first.  We  will 
assume  that  Oxy 2  is  the  inertial  reference  frame  and  OVyV 
the  non-inertial.  As  in  Figure  1-7  shown,  r  and  r'  are 
position  vectors  pointing  to  an  object  and  originating  from 
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r 


the  origins 
denotes  the 
measured  from 


of  the  coordinate  axes 
position  vector  of  the 
the  inertial  origin. 


correspondingly,  r. 
non-inertial  origin 


Figure  1-7 

Inertial  and  Non-inertial  Reference  Frames 
Translational  Motion 

The  relation  between  those  three  vectors  is  clearly 

-*  / 

r  r  r  +  T0 

By  differentiating  twice,  with  respect  to  time,  we  obtain 

r  =  r'  +  r.  (1-154) 

Now  since  we  have  assumed  that  Newton's  laws  hold  on  the 
frame  ,  we  may  write 

m r  =  F  (1-155) 

or  substituting  Equation  1-154 

mr'  =  F  -  tnf;  (1-156) 

where  'rn  is  the  mass  of  the  object  and  F  the  force  applied 
onto  it. 
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Equation  1-156  merely  says  that  we  need  the  extra  term 
mf.  to  describe  the  translational  motion  of  the  object  in  the 
non-inertial  frame. 

If  Z  -0  then  the  frame  OxyV  is  inertial  and  Newton's  laws 
hold.  In  this  case.  Equation  1-156  takes  our  familiar  form 

m?'  =  F  (1-157) 
which  is  similar  to  Equation  1-155. 

Let's  now  consider  the  rotational  case.  We  will  assume 
that  r*  =0  so  that  o'  is  the  same  as  0,  as  shown  in 
Figure  1-8. 


Figure  i-8 

Inertial  and  Non-inertial  Reference  Frames 
Rotational  Motion 


The  position  vector  r  can  be  expressed  in  the  inertial 
frame  as 


r  c  *  t  + 

or  in  the  non-inertial  frame 

-»  /A 

r  =  x  t  + 


b  4  +  2  k 


A 


(1-158) 

(1-159) 
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where  (*,i4,z)  and  (  x' ,  ,  z' )  are  the  corresponding  axes 

components  of  r . 

Similarly  the  vector  A  can  be  expressed  as 

A  •  A*  t  ♦  A9  $  ♦  A,  k  (1-16  0) 

or 

a  =  a;  i\  Ay  y  *  Az  k'  (1-16D 

Differentiating  Equation  1-160,  1-161  with  respect  to 
time,  we  obtain 


dA 

dAx  A 

d Ay  a  dAt  * 

dl  = 

dt  1  * 

1  4-  -  It 

dt  dt 

(1-162) 

and 

(#)' 

dAx  *  > 
'  dt 

dAy  dAz  *  / 

‘sr*  *irk 

(1-163) 

where  ( Y  denotes  that  we  are  taking  the  time  derivative 
of  A  on  the  primed  non-inertial  system  where  the  unit  vectors 
i  ,  k  are  constants. 


We  are  interested  now  in  the  time  derivative  of  the 
vector  A  relative  to  the  inertial  system,  but  we  would  like 
to  differentiate  Equation  1-161  instead  of  Equation  1-160, 
obtaining  thus: 


dt  cH  dt  j  dt  k  dt.  dt  dt 


(1-164) 


or 


d A  /dA^  .  /  dT  .  ,  dj'  k'  dk 

■ar’lir)  5r*^dr>Ai5r 


dT  i  '  dj'  ^ 


dt  '  dt  T  dt  (1-165) 

It  can  be  shown  that  a  more  convenient  expression  for  the 
three  last  terms  is  given  by 
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where  u)  is  the  angular  velocity  vector  of  the  primed  system 
which  is  assumed  to  rotate  about  an  axis  having  the  direction 
of  u>  ,  through  the  origin,  and  with  angular  velocity  u3, 
relative  to  the  unprimed  inertial  frame. 

For  proof,  assume  for  a  moment  that  ^  is  a  vector  "at 
rest"  on  the  primed  system  so  that 

(c^/cU)'  =0  (1-167) 

We  shall  show  that 

d^t  =ujxA  (1-168) 

Using  the  time  derivative  definition,  we  can  write: 

d  A  t  A  U*  At) -Aft)  -  7 

- - -  -  ».A  (1-169) 

Figure  1-9  shows  two  successive  positions  of  the  vector  A 
which  rotates  with  angular  velocity  w  with  respect  to  the 
inertial  frame,  so  that  the  tip  of  the  vector  traces  a  circle 
perpendicular  to  the  direction  of  (5  . 

To  show  Equation  1-169,  we  must  show  that  the  directions 
and  magnitudes  of  the  vectors  on  the  left  and  right-hand  side 
are  the  same. 

Notice  that  the-  direction  of  the  vector  dA/dt  ,  by 
Equation  1-169,  is  the  same  as  the  direction  of  the  vector 
A(t+At)-  A(i)  ,  which  is  the  same  as  the  direction  of  the 
vector  JxA  by  looking  at  Figure  1-9. 


A(4*M)  -  A(4) 

rc uad 


Figure  1-9 

Time  Derivative  of  a  Rotating  Vector 


The  magnitude  of  the  (left-hand  side  of  Equation  1-169 


.  |A(-UM.)  -  A  Ct)  |  .  fAsmd>)(u>Al) 

tim  77 -  =  “,r”  - 77 - 

Al-0  Ai  M+O  M 


(1-170) 


=  A  u>  sm<$> 


(1-171) 


which  is  the  same  as  the  magnitude  of  the  right-hand  side 
vector,  liox  A| . 

Hence,  both  directions  and  magnitudes  are  the  same,  so 
Equation  1-168  holds  and  if  A  is  not  stationary  with  respect 
to  the  non-inertial  system,  we  have  according  to 
Equation  1-165 
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+  co  x  A 


(1-172) 


cM  idXy 

dt  1 

which  says  that  the  time  derivative  of  a  vector  A  relative  to 
an  inertial  frame  equals  the  time  derivative  of  the  same 
vector  measured  at  the  non-inertial  frame  plus  the  term  u>*A 
where  u>  is  the  angular  velocity  of  the  non-inertial  as  seen 
by  the  inertial  frame. 

Combining  translational  and  rotational  motions,  one  can 
easily  arrive  at  the  equation: 

dA  /ciA\/  -  r  - 

+  a)xA  +  (1-173) 

where  7.  is  the  linear  velocity  of  the  non-inertial  with 
respect  to  the  inertial  frame. 

One  step  further,  differentiating  once  again 
Equation  1-172,  with  respect  to  time,  we  obtain 


which  says  that  the  second  time  derivative  of  a  vector  A 
relative  to  an  inertial  frame  equals  the  second  time 
derivative  of  the  same  vector  measured  at  the  non-inertial 
frame,  plus  some  extra  terms  involving  5,  A  ,  dw/dl  ,  and 
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( dA /di  )  where  uj  is  again  the  angular  velocity  of  the 
non-inertial  as  seen  by  the  inertial  frame. 

Combining  translational  and  rotational  motions  one  can 
easily  arrive  at  the  equation 


(1-176) 


where  r»  is  the  linear  acceleration  of  the  non-inertial  with 
respect  to  the  inertial  frame. 

A  special  significance  of  Equations  1-173  and  1-176  is 
given  by  substituting  A  with  the  position  vector  7  of  an 
object.  They  are  taking  now  the  form 


dr  dr 


+  ouxr  +  r. 


(1-177) 


V  3  (  V)'  +  cux r  + 

where  V  is  the  velocity  of  the  object,  and 


(1-17  8) 


(1-179) 


d  =  (a  )'  +  2  <3  x(\7)  *  u3  x  (uix  r )  +  uJxr  +  (1-180) 

where  a  is  the  acceleration  of  the  object. 

One  can  recognize  that 

2a3x(v’)/  s  corriolis  acceleration  (1-181) 

and  exists  if  the  object  has  a  velocity  (V)/  with  respect  to 
the  non-inertial  frame.  Also, 


u)x(u)xr)  “  centripedal  acceleration 


(1-182) 
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which  is  always  directed  toward  the  axis  of  rotation  and  is 
perpendicular  to  that  axis,  as  shown  in  Figure  1-10. 


Figure  1-10 

Centripetal  Acceleration 

Since  we  have  assumed  that  the  unprimed  system  is 
inertial,  we  may  write 

Tna  =  F  (1-183) 

and  substituting  Equation  1-180 

'rn(d)/  *  F  -  2-m[uix('/)']  -'Yri[un<(i&vr)]-/rn(u)xr)  -  wa0  (1-184) 
where  w  is  the  mass  of  the  object  and  F  the  force  applied  to 
it. 

Equation  1-184  merely  says  that  we  need  four  extra  terms 
to  describe  the  translational  and  the  rotational  motion  in  a 
non- inertial  frame. 

It  is  clear  that  the  forces: 
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2'w[(i3x(V),J  »  corriolis  force 

(1-185) 

'YT>  [  (5x((Sxr)j  s  centripedal  force 

(1-186) 

'■mCuSxr)  a  Euler  force 

(1-187) 

‘wia.  *  D’  Alembert  force 

(1-188) 

are  not  true  forces.  They  only  exist  due  to  the  fact  that 
our  object  is  moving  on  a  non- inertial  frame. 

According  to  our  first  assumption  that  the  earth  is  the 
inertial  frame  of  reference/  we  consider  that  it  is  rotating 
with  fixed  angular  velocity  u>  so  that  we  can  write 
Equation  1-184  as: 

-rtf  (a)'  =  F  -  Zm  a}  x  (V)  -  m  uJxfioxr  )  -  wa,  (1-189) 

We  have  seen,  up  to  nov«  the  effects  of  considering  a 
non-inertial  frame  of  reference  and  our  next  step  will  be  to 
orient  our  flight  vehicle  in  that  frame.  The  reason  for 
selecting  a  non-inertial  frame  will  be  greatly  appreciated 
when  one  will  see  how  easy  it  is  to  formulate  equations  of 
motion  of  a  body  attached  in  this  frame  and  moving  with  it. 

For  flight  control  applications,  it  is  appropriate  to 
select  the  origin  of  the  axis  system  on  the  center  of  mass  of 
the  vehicle. 

From  then  on  we  will  orient  the  nose  of  the  vehicle 
towards  the  positive  x-axis,  which  will  also  be  a  symmetry 
axis.  The  positive  y-axis  will  point  to  the  right  and  the 
positive  z-axis  will  point  downwards  having  the  direction  of 
gravity,  if  one  considers  straight  level  flight  as 
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Figure  1-11  shows,  so  that  the  xz  plane  will  be  a  plane  of 
symmetry. 

Vector  quantities  associated  with  this  axis  system  are 
shown  in  Figure  1-11  and  are  described  as  follows: 


Figure  1-11 
Axes  System 


Velocity  of  the  center  of  mass  with  respect  to  the 


inertial  frame  measured  on  this  axis  system 

V  =  UL4-Yj*wlc 


(1-190) 


Angular  velocity  about  the  center  of  mass 


Aerodynamic  force  acting  on  the  center  of  mass 

F  -  X  t  +  Yj  t-Zk 

Moment  about  the  center  of  mass 

G  =  + 

Directions  as  shown  in  Figure  1-11 
Also,  we  can  define:  (see  Figure  1-12) 


Angle  of  attack: 


ot  »  -tan"1  ( w/u) 


(1-192) 


(1-193) 


(1-194) 


Figure  1-12 

Angle  of  Attack  and  Angle  of  Sideslip 


Angle  of  sideslip 


p  •  sirf'  ( Vv) 


(1-195) 


The  reference  axis  system  has  been  adopted  by  B.  Etkin, 
D.  McRuer  and  others,  and  is  usually  referred  to  as  body  axis 
system  (see  references).  The  notation  used  is  quite  common 
on  most  references. 


Another  axes  system  adopted  by  C.  Perkins  is  shown  in 
Figure  1-13.  It  is  usually  referred  to  as  wind  or  stability 
axes  system. 


Figure  1-13 
Wind  Axes  System 

In  this  system  the  x-axis  points  to  the  opposite 
direction  to  that  of  the  relative  wind.  Notice  that  in  the 
undisturbed  flight  the  body  axes  coincided  with  the  wind 
axes,  but  in  the  disturbed  flight  that  may  or  may  not  happen. 

1.7  RIGID  BODY  DYNAMICS 

Rigid  body  is  a  system  of  particles  in  which  the  distance 
between  any  two  of  them  is  not  changing  regardless  of  the 
acting  forces.  For  our  flight  vehicle  the  rigid  body 


fey  -  \  ...i 
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assumption  implies  that  as  far  as  its  motion  is  concerned  in 
three  dimentions,  only  negligible  errors  are  introduced. 

We  will  consider  first  the  planar  motion  of  rigid  body 
where  we  will  find  the  kinetic  energy  and  the  angular 
momentum  expressions. 

Assume  that  the  rigid  body  in  Figure  1-14  can  move  in  the 
xy-plane  and  can  rotate  about  an  axis  perpendicular  to  that 
plane  passing  through  the  center  of  mass  c. 


■ 


Figure  1-14 

Rigid  Body  -  Planar  Motion 

♦ 

We  denote  70  ,  the  position  vector  of  the  center  of 
mass  c,  7  the  position  vector  of  any  particle  on  the  rigid 
body  and  r'  the  vector  from  the  center  of  mass  to  that 
particle. 

Usually,  the  mass  of  the  rigid  body  is  not  uniformally 
distributed  over  the  entire  volume  V.  It  is  convenient  to 


think  for  a  moment  that  the  entire  volume  of  the  rigid  body 
containing  a  mass  'ffi,  is  projected  on  the  xy-plane  on  an 
area  A.  Let's  denote  the  mass  per  unit  projected  area  of 
the  body.  It  is  clearly  seen  that  the  summation  over  the 
entire  area  of  all  those  masses  gives  the  total  mass  of 
the  body.  This  mathematically  is  expressed  as: 


(1-196) 

From  Figure  1-14  we  have: 

r  =  r.  *  V  (1-197) 

For  a  differential  mass  *idA  at  which  the  f  vector  points, 
the  kinetic  energy  is  expressed  as  j^dAr2  .  Integration  over 
the  entire  area  A  gives  the  total  kinetic  energy  of  the  rigid 
body. 


»i|>ndA 


(1-198) 


The  quantity  f1  according  to  Equation  1-197  is  expressed 


f2  =  rr  *  C r.  +r')-(i \+r')  =  r.2  tr'SZr.r' 
Then  Equation  1-198  becomes 


(1-199) 


(1-200) 


(1-201) 
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We  now  express  the  velocity  r  as  a  sum  of  two 
velocities : 

7'  dr'  ^  A  a 

f  *  dT  7  *  ry  *  (1-202'* 

where  d r'/dt  is  the  radial  velocity,  9  is  the  angular 
velocity  of  the  body,  r'0  the  tangential  velocity  of  the 

A 

particle  ,  t  and  9  denote  the  corresponding  unit  vectors  in 
polar  coordinates. 

Noticing  that  dr'/dtsO,  Equation  1-200  is  written 


K  =  .Lg'J  r'\dA 


(1-203) 


We  define 


I  =  /  r^dA  (1-204) 

JA 

the  moment  of  inertia  of  the  body  about  an  axis  perpendicular 
to  the  plane  of  motion  through  the  center  of  mass.  Then 

H 4  (1-205) 

This  equation  merely  says  that,  in  addition  to  the  term 

^  nn'r.1  that  describes  the  translational  motion,  we  need  the 

extra  term  2.  T9*  which  describes  the  rotational  motion  of 
2 

the  body  around  an  axis  perpendicular  to  the  plane  of  motion 
through  the  center  of  mass. 

The  angular  momentum  of  the  same  differential  mass  ijdA 
about  the  origin  is  expressed  as 

rx^dA?  (1-206) 
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Integration  over  the  entire  A  gives  the  total  angular 
momentum  of  the  rigid  body. 

=jW)n<«  a-207) 

The  cross  product  term  can  be  written 


r  xr  =  (f;  +  r')x  (r,*r')  =  r0  x  i;  +  r'x ?'  +  r'xZ  +  ra  xr' 

Then  Equation  1-207  becomes 


(1-208) 


fl  - 

In  this  expression  the  integral  of  the  last  term  is  zero 
as  Equation  1-201  shows.  By  a  similar  argument  the  third 

•  •  A 

term  is  also  zero.  Noticing  that  r'=r'03  ,  Equation  1-209 
becomes 


r.xr,  f r[dA  +  f ( T \ r)  i|dA  *  r„  v  f  r'rfdA  +  raxf  r  VjdA 
-4  Jk  Jk  -4 


(1-209) 


h  =(x0xr^)<m  +Q  r'(r'x  i))^dA 


(1-210) 

If  we  use  the  center  of  mass  at  the  origin  (^  =0)  we  are 
left  with: 


h  .e£r(T*$)n<n  (i_2n) 

The  cross  product  in  the  integral  indicates  that  the 
angular  momentum  vector  is  perpendicular  to  the  plane  of 
motion.  Since  (rx!ll  =  r,  the  magnitude  of  the  angular 
momentum  is  given  by: 


h  =  0  \t\6K  *10 

■4 
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(1-212) 


v 


So  the  magnitude  of  the  angular  momentum  of  a  rigid  body 
rotating  in  a  plane.  t.bout  an  axis,  through  the  center  of 
mass,  is  given  by  Equation  1-212.  The  angular  momentum 
vector  points  in  the  direction  of  the  axis  of  rotation  that 
can  be  easily  found  by  the  right  hand  rule. 

In  summary,  we  have  seen  that  the  motion  of  a  rigid  body 
in  planar  motion  can  be  described  by  a  translational  motion 
of  the  center  of  mass,  plus  a  rotational  motion  about  an  axis 
through  the  center  of  mass. 

We  will  now  look  at  the  motion  of  a  rigid  body  in  a  three 
dimensional  space.  It  generally  has  six  degrees  of  freedom 
if  not  constrained.  Hence,  we  need  six  coordinates  to 
specify  the  motion.  The  three  of  them  can  describe  the 
translational  motion  of  the  center  of  mass  of  the  body  and 
the  remaining  three,  the  rotational  motion. 

We  consider  two  coordinate  systems.  An  unprimed  system, 
as  shown  in  Figure  1-15,  whose  origin  is  0,  for  simplicity  is 
fixed  on  the  rigid  body  but  its  axes  do  not  rotate  with 
respect  to  an  inertial  frame,  and  a  primed  system  whose 
origin  and  axes  are  fixed  on  the  rigid  body. 

We  also  assume  that  the  rigid  body  has  a  mass 
distribution  per  unit  volume  (density)  which,  in  general,  is 
a  function  of  the  position  vector  r  to  a  particle  on  the 
rigid  body  having  a  differential  volume  civ  .  The  total  mass 
is  given  by: 
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Figure  1-15 

Rigid  Body  -  Space  Rotational  Motion 

m  =I  PCr)civ  (1-213) 

where  p(r)  is  the  density  as  a  function  of  r  and  dv  a  small 
differential  volume  where  r  points  to. 

For  the  differential  volume  dv  the  angular  momentum  for 
rotation  about  the  origin  of  the  unprimed  system  is  given  by 

r x [ pCr)dv  f ]  (1-214) 

Integration  over  the  entire  volume  gives  the  total  angular 
momentum:  „  f 

h  Cf^J  Prr)d'  (1-215) 

where  f,  ?  are  measured  on  the  unprimed  system. 

Let  to  the  angular  velocity  of  the  rigid  body  about  the 
origin  of  the  unprimed  system  and  r  the  position  vector  from 
the  primed  system. 
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It  is  more  convenient  to  express?,  which  can  be  seen  as 
the  velocity  measured  on  an  inertial  frame,  as 

7  s?'  +  a>  x  7  (1-216) 

Notice  that 

'  dt  1  *  (1-217) 

A  ^ 

where  z‘  and  i)  are  unit  vectors  in  polar  coordinates  of  the 
primed  system. 

It  is  clear  that  dr  /di  =  0  according  to  rigid  body 
assumption,  and  9=0  because  there  is  no  relative  rotation  of 
the  differential  volume  dv  that  can  be  seen  by  the  primed 
system.  In  other  words,  ?'  is  a  vector  at  rest  in  the  primed 
system  so  that  we  can  write 

r  -  cio x  r  (1—218) 

Substituting  in  Equation  1-215) ,  we  obtain 
h  -  j"  r  x  ( (ixr )  p(r )  dv 
Using  now  the  tripple  product  identity 
A*($XC)  -  (A  C)  B  -  CA-B  )  C 
we  can  write 

rx(ujxr)  s  rJ(jT  -  (ai-r)  r 
so  that  Equation  1-215  becomes 

h  [  r*  co  -  Twrjr  j  p(?) dr 
In  component  form  we  can  write 

ht  *  £  -Cai-r)  rj  pcf)dv 

where  1=1,  2,  or  3  denotes  the  x,  y,  or  z  component 
respectively.  The  dot  product  is  expressed  as: 


(1-219) 


(1-220) 


(1-221) 


(1-222) 


(1-223) 
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(1-224) 


(1-225) 


(1-226) 


h  "1 1  r*5v  -nol  {1_227) 

as  the  ij  th  component  of  the  inertia  tensor  I  about  the 
origin  of  the  unpriraed  system. 

So  finally,  we  express  the  ith  component  of  the  angular 
momentum  as: 


3 

^ 


(1-228) 


In  the  expanded  matrix  form  this  gives: 


h*  I**  Ixy  Ixr  U,, 

hij  s  Ivty  Jyy  Jyz  Wy  (1—229) 

Hz  I  xz  Jyz  Izz  vz 

This  result  is  compared  with  the  result  of  the  planar 
motion  case  that  is  expressed  by  Equation  1-212,  where  the 
inertia  scalar  term  is  expressed  now  as  a  matrix. 

The  components  of  the  inertia  tensor  are  given  by 
Equation  1-227  in  cartesian  coordinates  as: 
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(1-23  0) 


=  iff  (*2*-y*+z*)  -**]  P^)  dxdydz  = 

=Jjj  (yz  +  z2)  p(r)  dxdydz 


Similarly, 


Also, 


Jyij  -  jjj'  (  xJ4  z*)  p(f)  dxdydz 

(1-231) 

I  zz  *jjj  (  x2*y2)  p C  r  )  dx  dydz 

(1-232) 

I Xij  =  iijx  *-fff  *y  p(r)  dxdydz 

(1-233) 

fxz  -  Izx  *z  pc^)  dxdydz 

(1-234) 

Examples  Calculate  the  inertia  tensor  of  a  cube  of 
side  a,  mass  m  and  uniform  density  p  about  one  of  its  corners 
(Figure  1-16) . 


Figure  1-16 

Inertia  Tensor-Calculation  of  a  Cube 


h*  *  efff  (y2+z2)  <ixe!yeiz  i  tj'+z*)  JzJtjj* 


The  eigenvalues  and  the  eigenvectors  of  an  inertia  tensor 
have  a  very  special  significance.  The  eigenvectors  determine 
the  directions  of  the  so-called  principal  axes  of  a  rigid 
body  that  pass  through  the  point  where  the  inertia  tensor  is 
calculated  (usually  the  center  of  mass).  The  corresponding 
eigenvalues  are  called  the  principal  moments  of  inertia. 

It  can  be  shown  that  if  a  rigid  body  rotates  about  one  of 
its  principal  axes,  its  angular  momentum  about  the  point 
which  the  inertia  tensor  is  calculated,  is  parallel  to  its 
angular  velocity. 

For  proof,  let  a  principal  axis,  and  B  the 
corresponding  principal  moment  of  inertia.  The  eigenvalue 
problem  equation  which  has  to  be  satisfied  is: 

r  £  *  "B  f  d-240) 

or 

Iij  ft  -VPi  (1-241) 

Having  rotation  about  a  principal  axis  simply  means  that 
u>  has  the  direction  of  the  principal  axis  which  is  expressed 
as: 

U)  *  c  f&  (1-242) 

where  C  is  a  constant.  In  fact,  if  p>  is  unit  vector,  C  is 
the  magnitude  of  the  angular  velocity  u5. 

Now  the  ith  component  of  the  angular  momentum  becomes 
=•  I ij  (*> j  =■  ItJ  CcjSj)  «  c  Iij  ft  *  c  S£t  =  SCcjSt )  (1-243) 

and  finally, 

k;  *  B  cot 
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(1-244) 


In  vector  form 


-  h  =  B  3  (1-245) 

*♦ 

which  proves  that  h  is  parallel  to  <3. 

It  can  also  be  shown  that 

1.  Any  plane  of  symmetry  of  a  rigid  body  is 
perpendicular  to  a  principal  axis. 

2.  Any  axis  of  symmetry  of  a  rigid  body  is  a 
principal  axis. 

The  density  p  may  not  be  constant.  For  example,  if  the 
xy  plane  is  a  plane  of  symmetry,  that  simply  means  that 

p(-*,y,z)  «p  C*,y,z)  (1-246) 

and  if  the  z-axis  is  an  axis  of  symmetry,  that  simply  means 
that 

p(-y,-y,z)  -  p(y,y,z)  (1-247) 

3.  Since  the  inertia  tensor  of  a  rigid  body  is 
symmetric,  (i.e.  l>y*Jyx  etc.),  the  principal  axes  are  always 
orthogonal  to  each  other  (Important  result  from  Linear 
Algebra) . 

4.  If  two  or  three  eigenvalues  of  a  symmetric  tensor 
are  the  same,  there  exist  two  or  three  orthogonal 
eigenvectors  associated  with  that  eigenvalue. 

Example:  The  inertia  tensor  of  a  rigid  body  has  been 
calculated  with  the  reference  point  at  the  center  of  mass  and 
has  been  found  to  be 
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(5 

C 

-JO 

0 

25 

0 

-<o 

0 

30 

axes 

and 

the 

Find  the  principal  axes  and  the  corresponding  principal 
moments  of  inertia. 

Solution: 

Let  j)  denote  the  eigenvector  (principal  axis)  and  A  the 
corresponding  eigenvalue  (principal  moment  of  inertia) . 

The  following  eigenvalue  equation  has  to  be  satisfied 


IJ  = 

AI 

(1-248) 

lich 

1  A- 15 

0 

JO  I 

> 

^4 

it 

O 

<0 

o 

A-25 

o 

II 

O 

(1-249) 

1  JO 

0 

A- 50  | 

(A-i5)( A-  25)  [A- 50)  +  J0[. 

JOfA-25)]  *  O 

(1-250) 

+  (A-25)  [  (  A-J5)  (A-50) 

-  400]  =o 

(1-251) 

=?>  (A- 25)  ( A  -  35)  C  A-JO) 

3  0 

(1-252) 

the  principal  moments 

of 

inertia  are: 

A,  =  10  A2  -  25 

A,  r  55 

(1-253) 

To  find  the  eigenvectors: 

For  Aj  »10  Equation  1-249  becomes 


-5 

0 

40 

0 

-(5 

0 

5 

(1-254) 

JO 

0 

-  20  | 
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which  yields 


j 


+  =0 

-  =  0 

=> 

r  & 

-20^  = 0 

& 

X 

II 

(1-255) 


The  eigenvector  that  satisfied  Equation  1-254  is 


2 

o  (1-256) 

i 

It  is  more  convenient  to  express  the  eigenvector  as  a 
unit  vector 

m 

5  «  0  (1“257> 

*/6 


Similarly,  we  find  that 
0 

-  J  ass  os  i  a  led  wi4h  A2  *  25  (1-258) 

0 


and 

Vis 

=  0  assosicrled  w/i4lt  A3  =  55  (1-259) 

It  is  important  to  note  that  the  above  eigenvectors  are 
orthogonal  since  they  are  coming  from  a  symmetric  tensor. 

Since  now  the  principal  axes  of  a  rigid  body  are 
orthogonal  to  one  another,  or  can  be  chosen  orthogonal,  it  is 


i  . 

\  i 


very  convenient  to  have  them  coincide  with  the  primed  system 
axes,  the  origin  of  which  can  be  moved  to  the  origin  0  of  the 
the  unprimed  system.  The  importance  of  this  arrangement  will 
be  shown  a  little  later. 

The  kinetic  energy  of  the  differential  volume  dv  is 
2  (r  r)p(f)dv.  Integration  over  the  entire  volume  gives  the 


total  kinetic  energy: 

K  =  4  [  tr  r) 

*  Jy 


(1-260) 

Substituting 

7  =  uix?  (1-261) 

we  obtain 

K  r  (u>xr)(<uxr)  pCrjdx  (1-262) 

Using  Equation  1-224,  1-226  and  the  vector  identity 

(  A  X  B ) (  C  X 3)  =  CX C)(i5)  -  (X1>)  (  B  C)  (1-263) 

we  can  write 

3  \  /  3 


(i5x? 


)(u>xr)  *  SijWiY^  5ijo)j)-f 

ul  '  J*J  '  '  t*J  /  V  j.j 


so  that 


l-*J  U'J 


ci  1^-  fjr‘5<j-riri]“|j  |“i 


(1-264) 


(1-265) 


and  finally. 
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J 


—  I  3  \  - 

K-J  2.(1  (1-266> 


u-1  >1 


Notice  that  if  the  body  is  rotating  about  a  principal 
axis  (  h  II uj  )  its  kinetic  energy  is  maximum  K  =*  j  hu> .  That 
happens  when  the  rigid  body  rotates  about  a  principal  axis. 

This  compares  with  the  result  of  the  planar  motion  case 
that  is  expressed  by  Equation  1-205  ,  if  we  neglect  the 
translational  motion. 

We  are  now  in  a  position  to  formulate  the  rigid  body 
equations  of  motion 

The  translational  motion  equation  is  given  by  the 
equation  of  motion  of  the  center  of  mass  of  the  rigid  body 

wra  r  F  (1-267) 


f  dV  -  (7 )  -r 

—  +  cv  x  V  i  =  F 
l  / 


(1-268) 


The  rotational  motion  about  the  point  (O)  which  is 
conveniently  chosen  as  the  center  of  mass  (c) ,  as  Figure  1-17 
shows,  is  given  by 


dh  -r 

ST 


(1-269) 


where  G  is  the  net  external  moment  about  the  center  of  mass 
(c)  and  H  is  given  by  Equation  1-228. 
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Figure  1-17 

Rigid  Body  -  Space  Motion 

The  vector  dk/cH  is  measured  in  the  unprimed  system,  but 
it  can  be  expressed  as: 

dk  felCV  -  r 

dt  -IcR  )  +  Vxk  (1-270) 

where  (c^/di)  is  the  same  quantity  measured  in  the  primed 
system,  and  <2  is  the  angular  velocity  about  the  center  of 
mass. 

So  Equation  1-269  is  written  in  component  form  as: 

'dk,V 


(x)  +  ^  =  ^ 


(1-271) 

where  i  refers  to  a  primed  axis  now.  Or  better,  using 
Equation  1-228 
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(1-272) 


=q> 

j-'i 

where  I  tj  is  a  component  of  the  inertia  tensor  calculated 
along  the  primed  axes. 

Notice  that  the  primed  axes  unit  vectors  are  vectors  at 
rest,  in  the  primed  system,  and  obviously  Iij  is  not  changing 
with  time.  Hence, 

X  ^lJ  Hi  +  ^X^1  m  (1-273) 

j-i 

This  result  is  interpreted  as  follows:  Suppose  the  rigid 
body  rotates  about  a  principal  axis  with  constant  angular 
velocity  o> .  Rotation  about  a  principal  axis  implies  ui 
parallel  to  ^  ,  so  that  3*^*0.  This  simply  means  that  the 
moment  required  for  this  rotation  is  zero. 

This  compares  with  the  planar  motion  case  result  which 
comes  by  taking  the  time  derivative  of  Equation  1-212. 

K  =  I  ctJ  *  G  (1-274) 

Equation  1-27  3  is  the  equation  of  rotational  motion  in 
the  general  case  for  rotation  about  a  random  axis,  where 
external  moment  is  required  even  for  constant  angular 
velocity. 

It  is  now  clear  why  we  formulate  the  equations  of  motion 
in  the  primed,  non-inertial  coordinate  system.  We  have 
avoided  the  time  differentiation  of  the  inertia  tensor  in 
Equation  1-272  by  paying  a  little  for  the  extra  term  toxk. 


'  •  *  *  -  . 
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In  summary ,  we  found  the  equations  of  motion  of  the  rigid 
body  described  in  the  convenient  form  which  will  be  used  to 
derive  the  flight  vehicle's  equations  of  motion  in  Chapter  2. 
They  are: 


fdJ  -  r,\ 

W  —  +  u)x  V  *■ 

'  di  I 


(1-275) 


and 


Iij  ^  ♦  Ctuxk)t  *  G, 


dl 


(1-276) 


where 


3 

^  =  Z_  ^lJ  (1-277) 

j-i 

and  Itj  is  calculated  about  the  center  of  mass  with  respect 
to  the  primed  axes. 

Finally,  if  the  primed  axes  are  principal  axes  of  the 
rigid  body  it  is  not  difficult  to  see  that  the  inertia  tensor 
has  only  diagonal  terms  representing  the  principal  moments  of 
inertia  (eigenvalues).  This  result  comes  from  the 
diagonalization  process  of  a  square  matrix.  So,  in  this  case 
Equation  1-276  is  simplified  as: 


h  ST  ♦ 


where 


k  *  Iii 


Gi 


(1-278) 


(1-279) 
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CHAPTER  2 


YSBICLE  IS  EQUATIONS  OF.  MOTION 
2.1  INTRODUCTION 

With  the  methods  of  obtaining  rigid  body  equations  of 
motion,  established  in  the  preceding  chapter,  we  now  proceed 
in  this  chapter  to  formulate  the  equations  of  motion  in  a 
most  generally  useful  form. 

The  equations  are  basically  simple  in  the  sense  that  they 
are  derived  by  making  use  of  the  basic  Newtonian  laws  for 
forces  and  moments.  Great  advantages  are  obtained,  as  we 
have  seen,  by  using  a  non-inertial  frame  of  reference, 
because  with  an  axis  system  fixed  in  the  body,  the  inertial 
properties  are  constant  and  that  simplifies  the  moment 
equations. 

Some  penalties  are  obtained  by  introducing  non-inertial 
frames,  the  most  important  of  which  is  the  addition  of  extra 
terms  in  the  basic  equations  of  forces  and  moments,  in  trying 
to  make  the  basic  laws  valid  in  this  environment. 

In  any  case,  although  the  equations  are  simple,  they  are 
highly  non-linear  and  not  generally  solvable  mathematically. 
The  formulation  is  greatly  simplified  by  the  use  of 
linearization  techniques,  and  small  pertubation  theory. 
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This  method  assumes  an  equilibrium  position  which  is,  in 
fact,  a  solution  to  these  equations.  That  allows  for  the 
changes  in  forces  and  moments  to  be  expressed  as  linear 
functions,  in  terms  of  small  pertubations  in  linear,  angular 
velocities,  and  other  parameters  involved. 

The  equations  of  motion  will  be  derived  for  the  body  axis 
system  established  in  Figure  1-11  and  will  be  left  in  their 
dimentional  form.  Non-dimensional  form  will  be  mentioned  as 
it  applies  to  body  axis  system  (Etkin's)  and  wind  axis  system 
(Perkin's). 

Finally,  the  equations  of  motion  will  be  split  into  two 
separate  groups,  the  longitudinal  and  the  lateral  groups  of 
equations. 

2.2  FORCE  EQUATIONS 

Starting  from  the  vector  force  equation 


IF 

=  m  [ 

dV 

—  H 

■  LUX 

V) 

di 

(2-1) 

:f  = 

1X1  ♦ 

2Y i 

+  11 

A 

Ic 

(2-2) 

v  = 

ul  + 

A 

VJ  + 

A 

W  Ic 

(2-3) 

u)  - 

'p'i  + 

A 

V  + 

A 

r  k 

(2-4) 

we  can  derive  the  force  equations  in  component  form. 
The  cross  product  term  is  expanded  as: 


letting  now  St-»o  and  neglecting  the  higher  order  term  Swd©  , 
we  obtain 

ax  =■  u  +  vwc^  (2-11) 

and  for  the  corresponding  force  equation 

Z  y  =  mi  (b  +  wc^)  (2-12) 

which  is  the  same  as  Equation  2-6  by  neglecting  the  terra  vr  . 
Similarly,  the  acceleration  along  the  z-axis  is  expressed 

as: 

(w+Sw)  cosS©  -  (u+Su)  sinS©  -  w 

2  "  s  <2“13) 

which  is  simplified  as:^ 


or 


a,  = 


'aj+Sw  -  uS8  -  Si>89  -  vo 

St 


(2-14) 


az  -  **  ~<y>  (2-15) 

and  the  corresponding  force  equation  becomes: 

ZZ  -  (2-16) 

which  is  the  same  as  Equation  2-8  by  neglecting  the  term  pY . 

If  we  now  consider  the  motion  of  the  center  of  mass  in 
the  xy  plane  as  Figure  2-2  shows,  in  two  successive  different 
locations  on  the  flight  path  during  time  St,  we  can  derive 
the  acceleration  along  the  y-axis  as: 


as  - 


(vaSv)cosSv  »  (u»Sv)  sm&y>  -  y 

_ 


(2-17) 
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Figure  2-2 

Motion  of  the  Center  of  Mass  in  the  xy-Plane 


For  small  yaw  angle  ip ,  cosStp»l  and  sin&u>=8n>r  so  that 

„  (v+£v  u8a>  tbudv-v) 

a,  *  - -  (2_18) 

letting  now  <M.-*o  and  neglecting  the  higher  order  term  huhv  , 
we  obtain 

ay  -  v  +  v<"  (2-19) 

and  for  the  corresponding  force  equation 

IT  »  m  (v  *  ur)  (2-20) 

which  is  the  same  as  Equation  2-7  by  neglecting  the  term  pw. 

In  summary,  the  force  equations  are  simplified  neglecting 
higher  order  terms  as: 


IX  =  n\  C  u  +  yO 

(2-21) 

IY  =  'Wi(y-t-ut') 

(2-22) 

I Z  *  to  ( w  -  cjv  ) 

(2-23) 

2.3  VEHICLE  SYMMETRIES 

Dp  to  this  point  we  have  made  two  assumptions.  The  first 
assumption  was  that  earth  is  considered  to  be  fixed  in  space 
so  that  we  can  establish  an  inertial  reference  frame  on  it. 
The  second  assumption  was  that  the  vehicle  is  considered  to 
be  a  rigid  body. 

Facts  are  simplified  further  by  considering  the 
symmetries  of  the  vehicle.  For  an  airplane  or  a  missile/  the 
xz-plane  is  a  plane  of  symmetry/  if  one  assumes  that  the 
following  holds  for  the  density  distribution. 

p  (x,  y,z)  -  p  ( Vy ,0  (2-24 ) 

so  that 


Ixtj  ~  lyx  * 


-jjjp  *y  J^Jp(vy>z)*y 


(2-25) 


Notice  that  Equation  2-24  says  that  p  is  an  even  function 
with  respect  to  y  but  the  integrant  function  is  odd  because 
it  is  a  product  of  an  even  and  odd  function.  Since  the  y 
integration  has  symmetric  limits,  the  integral  is  zero. 

For  the  same  reason 


The  inertia  tensor 

Jxx 

0 


-Ixz 


(x,y,z)  dx 

z  therefore 

0  -1x2 

Iyy  0 

0  J« 


*  0 

(2-26) 

z  for  an  airplane  becomes: 

(2-27) 
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Notice  also  that  the  y  axis  is  an  eigenvector  of  the 
inertia  tensor  (Equation  2-27)  because  it  satisfies  the 
eigenvalue  matrix  equation: 

Ixx  0  -IyZ  0  0 

0  Iyy  0  J  =  lyy  t  (2—28) 

“Ixz  0  0  0 

Hence,  the  y-axis  is  a  principal  axis,  with  L/y  being  the 
corresponding  eigenvalue  or  principal  moment  of  inertia. 
This  was  expected  since  y-axis  is  perpendicular  to  the  plane 
of  symmetry  xz. 

Unfortunately,  this  is  the  only  symmetry  that  exists 
generally  in  any  conventional  aircraft. 

In  addition  to  this  symmetry  another  symmetry  can  be 
found  in  gun  projectiles  and  in  most  missiles  and 
space-crafts.  They  have  a  second  plane  of  symmetry,  the  xy 
plane,  so  that  one  can  assume 

p  (x,y,Z  )  =  p  C  x,y  ,-z)  (2-29) 

and  hence, 

Ixz  =  Izx  (2—30) 

Then,  the  inertia  tensor  for  a  two-symmetry  vehicle 
becomes 


Jxx  0  0 

0  Iyy  0  (2-31) 

0  0  Zzi 
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Notice  that  in  addition  to  the  y-axis  being  a  principal 
axis,  the  z-axis,as  well  as  the  x-axis,  are  all  principal 
axes  with  corresponding  principal  moments  of  inertia  Izz,  and 
Ixx ,  respectively. 


2.4  MOMENT  EQUATIONS 


Starting 

from 

the  vector  moment  equation: 

c  -’I  —  +(u3xi!) 

U  1  Al 

(2-32) 

or  in  component  form 

c. 

‘  1  lij  X1  +  Cwx  k)t 
j*j 

(2-33) 

where 

K  *  lij  “Jj 

(2-34) 

5  «  Li  t  Mj  »Nk 

(2-35) 

k  *  kx  ^  ^  J  k2  k 

(2-36) 

A  ’  Plf£U  *rt 

(2-37) 

Equation 

2-34 

is  expanded  as: 

ky  ^  pixy  -  r  Ixz 

(2-38) 

ky  = 

(2-39) 

kz  =  -p  Ixr  +  r  Jzz 

(2-40) 

also,  the  cross  product  terms  is  Equation  2-33  is  expanded  as 
(uJxli)  ®  -rlay )  1  +  (rhx-.pk1)j  ♦  ( pky  (2-41) 
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or 


A 

l 


or 


((2xk)  =(-<^1*2  t  Y  In  -r^IyiJ) 

*  (rp  I*x  -  r2IX2  +  p2Ixz  -rpln)  } 


+  (  Iyy  -  ^xx  t  r<^  Ixz  )  k 

(2-42) 

uixk)  =  [-tj^plxz  t  <jr  ( Iw  - lyy)]  t 

+  [  rp  f  Ixx- J22  )  *  1x2  (p2-r*)]  J 

+  i  (Jyy-lxx)  +  Ixz]  k 

(2-43) 

the  moment  equations  are: 

L  -  Ixx  p  -  I m  f  -  <^p  Ixz  +  CjT  (Izz-Iyy) 

(2-44) 

^  s  \  f  rp  (  Ixx  *  Izz )  +  Ixz  (  p2  - 1'2) 

(2-45) 

I'J  *  Izz  r  -  1x2  p  +  p^  ( Iyy  -Ixx )  +■  rc^  Ixz 

(2-46) 

2.5  EXPANSION  OF  GRAVITY  FORCES,  THE  LINEAR  AND  ANGULAR 
VELOCITY 

The  total  force  acting  on  a  flight  vehicle  can  be  split 
into  a  summation  of  aerodynamic  forces  (lift,  drag  and 
thrust)  and  a  body  force  which  is  the  gravity  force. 


i 
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Figure  2-3  shows  an  accidental  orientation  of  the  vehicle 

-* 

with  respect  to  the  vertical.  The  gravity  force  W  can  be 


expanded  in  components  along  the  body  axes  as  follows: 

Wx  =  -vn%  sin  §  (2-47) 

Wjj  =  mricj  cos  Q  sin  <t>  (2-48) 

W2  *  'nmj  cos®  cos  <J)  (2-49) 


Figure  2-3 

Expansion  of  the  Gravity  Force 


For  the  expansion  of  the  linear  velocity  vector  along  the 
body  axes  consider  in  Figure  2-4  the  inertial  reference  frame 
being  Ox, y,21  .  Transformation  to  the  body  axes  coordinate 
system  is  obtrained  by  rotating  the  reference  frame  in  an 
order-dependent  manner  as  follows: 

(a)  Rotate  around  reference  axis  0z{  by  the  yaw  angle 
which  defines  the  intermediate  reference  system  0*2y,  ^  ■ 

Om**!  • 
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(b)  Rotate  around  intermediate  axis  0y4  by  the  pitch 
angle  which  defines  a  second  intermediate  system  Ox,^,  - 

(c)  Rotate  around  final  axis  Ox  by  the  roll  angle  <j>  *  to 
define  the  desired  body  axes  system  0*yz  =  0x,y2  . 


Figure  2-4 

Expansion  of  the  Linear  Velocity 

The  components  of  the  velocity  vector  in  the  body  axes 
system  (U,V,W)  are  defined  in  terras  of  the  vector  components 
in  the  fixed  inertial  frame  (Uj ,V,  ,W, )  by  means  of  the  angles 
and  the  intermediate  vector  components  (Ui ,V2  ,W, )  and 

(Uj  f  ’vj  t Wj  )  • 
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The  following  relations  hold: 
for  the  first  rotation 

U2  *  Uj  cosip  +  Vi  simp  (2-50) 

V2  =  -Dj  smip  +  Vj  cosip  (2-51) 

W2  =  Wi  (2-52) 

for  the  second  rotation 

Uj~  03  cos  9  -  sin  0  (2-53) 

V»=  V2  (2-54) 

W3=  smS  +  W2  cos  0  (2-55) 

and  for  the  third  rotation 

U  -  l>3  (2-56) 

V  -  V3  cos  $  +•  Wj  sm  <p  (2—57) 


w  =  ~v3  +  W3  cos  (f?  (2-58) 

so  that  the  velocity  components  along  the  body  axes  are 
expressed  as 

U  *  Vi  cos  pcos  0  +  V,  smzpcos  0  -  W,  Sm0  (2-59) 

V  ~  Vj  (cos ip  sin  9  sin  -  siw(p  cos  <p )  *V,  (simp  sm0sin<p -1- cos ipcos<p) 

+  W,  cos  9  smip  (2-60) 

'A/  =  Dj  C  cosip  sm0  cos 4?  +  simp  siv>4>)  +V,  (smip  sm0  cos<p  -  cos  ip  sm  d> ) 

+  VP,  cos©  cos 4>  (2-61) 

Similar  approach  can  lead  us  from  the  body  axes  system  to 
the  inertial  for  which 

*  V  cos 9  cosip  +  V  (  siH(p  sin  0  cos  ip  -  cos4>  sm  ip  ) 

+  W(cos<p  sm0  cosip  +•  sirup  smip)  (2-62) 
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V,  =  DcosSsirny  +  V(  iiyi(|jsm0sm4>  +  cos^c av) 

t  W  (cos4>  sirtS  smy  -  s,n<j7  cos  ip)  (2-63) 

Wj  =  -  U  sinG  v  V  sm<fc>coi9  v  W  cos (2-64) 
Integration  of  Equation  2-62,  2-63  and  2-64  give  the 
expressions  for  the  translational  motion  but  in  other  than 
special  simple  cases  they  are  not  integrable. 

Following  a  similar  procedure  we  can  expand  the  angular 
velocities  p,q,r  into  components  along  the  body  axes. 

With  the  vehicle  in  horizontal  flight,  one  can  see  the 
yawing-velocity  vector  having  the  same  direction  as  gravity, 
the  rolling  velocity  vector  having  the  direction  of  the 
longitudinal  axis  and  the  pitching  velocity  vector, 
perpendicular  to  both  as  shown  in  Figure  2-5.  Those 
orientations  are  obtained  according  to  right-hand  rule. 


Figure  2-5 

Expansion  of  the  Angular  Velocity 
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Then  the  angular  velocity  vector  uj  can  be  expanded  in  the 


body  axes  as  follows: 

•  , 

P  =  <}>  -  Tj)  SW0 
r  9  cosp +•  ip  cos  ©  sin4> 

t  =  -0  sm<$>  +  ip  cos©  cos 43 


(2-65) 

(2-66) 

(2-67) 


One  can  now  solve  the  Equations  2-65,  2-66  ,  2-67  for  9, 


V,  <p 


Multiplying  Equation  2-66  by  cos  <f>  ,  Equation  2-67  by 


sin  <t>  ,  and  subtracting,  we  obtain; 


9  =  <^cos<|> -rsm4>  (2-68) 

Multiplying  Equation  2-66  by  sirJtp,  Equation  2-67  by 
cos  <f>  ,  and  adding,  we  obtain 


sw<t>  cos  db 

IP  =  9  — —  ♦  r  — — 

V  cos  ©  cos  © 


(2-69) 


Then  substituting  Equation  2-69  into  Equation  2-65,  we 


obtain 


4>  =  p  +  <^  sm4>  ian©  +•  rcos<p4ar»© 


(2-70) 


Integration  of  Equations  2-68,  2-69  and  2-70  give  the 
expressions  for  the  rotational  motion  but  in  other  than 
special  simple  cases  they  are  not  integrable. 

We  can  now  modify  the  vehicle's  equations  of  motion  by 
considering  the  gravity  as  follows: 


Force  equations: 


X  =  C  i>  -yr  +  g  sm©  ) 

Y  -  tn  (  v  +  ur  -  pw  _  ^  cos  0  sir1  $  ) 

Z  =  (  \m  +  pv  -cjv  Cos  &  cos  ip) 


(2-71) 

(2-72) 

(2-73) 


where  now  the  X,Y,Z  forces  represent  only  aerodynamic  forces 


(i.e.,  contributions  of  lift,  drag  and  thrust). 


Moment  equations: 

L  =  Jx*  p  -  Ixz  f"  -  C^p  Ixz  t  cyr  ( Izz  -  Iyy  ) 

(2-74) 

M  s  Iyy  +  f  p  (  I«  *  Izz  )  +  Ixz  (  p2  _r2  ) 

(2-75) 

N  *  Izz  f  -  p  +  (I^y-lxx)  f  Ixz 

(2-76) 

with  auxiliary  relations: 

'p  -  4?  -  ip sin0 

(2-77) 

-  ©cos#  +-  ip  cos  ©  sm# 

(2-78) 

r  =  -  0  sin  4>  +  ip  cos  ©  cos# 

(2-79) 

Notice  that  the  gravity  force  contributes 

zero  moments 

because  it  acts  at  the  center  of  mass. 

2.6  LINEARIZATION  OF  FLIGHT  VEHICLE'S  EQUATIONS 

The  three  force  and  the  threee  moment  equations  derived 
in  Section  2.5  must  be  linearized  for  use  in  the  stability 
and  control  analysis.  It  is  assumed  that  the  motion  of  the 
vehicle  consists  of  small  deviations  about  a  reference  steady 
flight  path. 

A  steady  flight  path  is  described  by  the  general 
equations  of  motion,  if  we  consider  all  time  derivatives 
zero.  We  denote  by  capital  letters  and  the  subscript  (o)  , 
the  steady  flight  conditions,  obtaining  thus: 

~  W  (QoWo  -  VoP„  4  ^  sin  0o  ) 

Y0  *  w  (  V0Zo  -  ?oW„  cos  ©o  sm  <£>. ) 

Z„  -  'TT1  Oo  Vo  -  Q.Vo  -  ^  COS  ©o  COS  <b0 ) 
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(2-80) 

(2-81) 

(2-82) 


(2-83) 

(2-84) 

(2-85) 


L0  r  “  QoPo  1x2  +  C  1 22  -  ) 

Mo  =  >U  (Ixx'Ia)  +  Ixz  (P/-*.*) 

No  =  %Qo  (Iyy  -Ixx)  +  P.Qo  Ixz 
By  subtracting  from  the  original  equations  the  above 
steady  flight  equations ,  we  obtain  the  perturbed  equations. 
For  example ,  subtracting  Equation  2-80  from  Equation  2-71  we 
gets 

X-X,  =  <1*  [  6  t  cjW  -  QoWo  -  vr  +VoR»  +  £  (  SW0  -  snr)0o)]  (2-86) 

Notice  that 

(2-87) 


Similarly 

vr  *  VoRo  *V.  Sr  +  R0§v  (2-88) 

Also 

sin0  =■  sm(0o+50)  =■  sm0o cos60 +cosS0  sm56  (2-89) 


For  small  perturbation  quantities.  Equation  2-89  becomes: 

%\r\B  -  sm0o  +( cos  Go)  6©  (2-90) 


Combining  all  these.  Equation  2-86  becomes 

6X  =  om[  VU»5c^  fQoSw  -V«6r  -P05v  +C|gm&)59J  (2-91) 

Although  we  are  deriving  perturbed  equations,  it  is  more 
convenient  to  drop  the  differential  symbol  § . 

Similar  work  leads  us  then  to  the  following  set  of 
linearized  equations: 
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X  *  %  l  0  +  4-Qo  w  -  Vo  r  -  Ro  v  +  ( yjj  0„)  ©J  ( 2-92 ) 

Y  *  on  [  v'  +  IU  +  £oV  -W»p  -P,w  -(geos 0»  cos  <£,)<$> *(gsm6Uir>&)e](2-93) 

Z  -  <m  [  *>  +  V8p  -t-fou  -  0o<^  -Q«o  +(g  cos£>  cos  4)4»  +(^sm€L cos 2-94) 

L  r  Ixx  p  "Ixz  r  -  (Qop  Jxr  +  (Qor  +  (?o^)(I«  *Iyy )  ( 2-95 ) 

M=  I^^  +  (Rftp  +  ?.r)(IyX-Izz)  +I«(2P0p-2Ror)  (2-96) 

Nl  =  Jzz  r  —  l*z  p  +  (P«^  +  Q.p)(Iyy  -Ixx)+  (&>^fQ0r)Jxz  (2-97) 

In  addition,  for  the  auxiliary  relations  described  by 
Equations  2-77,  2-78,  2-79,  the  corresponding  perturbed 
equations  are: 

-p  =■  -  TpstnOo  -  (Wo  cos  0„)@  (2-98) 

=  9  cos  <t>0  -  sin  Go  sm<fe)9  +  (y,ocos0o  cosdi  -  9„  sinta)  4>  + 

+  ip  cos  G0  sin  <&>  (2-99) 

r  s  ij>cos  Go  cos  4>o  -  (ty,  CDS  Go  SID  <k  +  0*  cos  <f>8)  <|>  -  0  sin  d?o 

-  (4>8  sm0o  cos  ^0)0  (2-100) 

Although  we  have  obtained  linear  equations,  they  are 
extremely  complex  and  not  suitable  for  control  analysis. 

A  further  simplification  is  obtained  if  we  consider  the 
reference  flight  path  to  be  straight,  i.e. 

TV.  *  bo  =0  (2-101) 

with  no  angular  velocities 

ye  -  Qo  ■  'Ro  =0  (2-102) 


wings  level  flight,  i.e 


symmetric,  i.e 


IPo  =0  (2-104) 

Vo  »  0  (2-105) 

With  this  assumption,  the  reference  flight  path  is 
described  by  the  equations: 


X„  =  'rncj  sm  0o 

(2-106) 

o 

H 

>-° 

(2-107) 

Zo  -  -  ©o 

(2-108) 

Lo  =  Mo  -  No  —  0 

(2-109) 

while  the  perturbed  equations  become 

X  =  om  [  u  +  \A/0  +  g  cos  Q0  O  J 

(2-110) 

Y  s/m^+Uor  +  ^p-gcosQ.^] 

(2-111) 

2  -  m  [  i)  -  U.c^  +  ^sin&©] 

(2-112) 

L  =  Jvx  p  -  Ixz  r 

(2-113) 

M  *  Jyij  ^ 

(2-114) 

N  =  IZi  f  -  Ixz  f 

(2-115) 

with  auxilliary  relations: 

p  =  cp  -  ip  sin  ©«. 

(2-116) 

v  ® 

(2-117) 

r  =  tp  cos  ©0 

(2-118) 

We  find  that  the  equations  are  separable  into  two  groups. 
The  first  group  associates  the  equations  that  describe  the 
motion  in  the  xz  body  plane,  called  longitudinal  equations 
and  the  second  group  associates  the  equations  that  describe 
the  motion  in  the  xy  body  plane,  called  lateral  directional 
equations.  g7 


We  find  also,  the  force  equations  divided  by  mass  and  the 
moment  equations  divided  by  the  corresponding  moment  of 


inertia. 

In  the  longitudinal  equations  the  independent  variables 

are: 

v  *  x-axis  linear  velocity 
w  =  z-axis  linear  velocity 
q  =  pitching  angular  velocity 
In  the  lateral-directional  equations  the  independent 
variables  are: 

v  »  y-axis  linear  velocity 
p  *  rolling  angular  velocity 
r  *  yawing  angular  velocity 

with  all  these  quantities  being  perturbations  from  the 
initial  flight  conditions. 

We  arrive  then  to  the  following  form  of  equations: 
Longitudinal  equations 


(2-119) 


~  Z  =  w  -  U.  ^  +  g  $in0o@ 

i 

Lateral -directional  equations 

—  Y  =  v  +  U,  r  -  V\Z0  p  -  g  cos  0O  4> 


(2-120) 

(2-121) 

(2-122) 
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r  Iw 


(2-123) 


I XX 

1  Kl  ■  Ixz  • 

irN  =  r'inp 

For  further  simplification/  it  is  convenient  to  introduce 
here  the  wind  axes  system/  which  as  we  have  seen  is  produced 
by  a  rotation  about  the  y-body  axis  so  that  the  x-wind  axis 
points  in  the  opposite  direction  of  the  relative  wind.  This 
axis  is  always  tangent  to  the  flight  path  and  makes  an  angle 
of  flight  path  with  the  horizontal,  as  shown  in 
Figure  2-6. 


horiuiitai 


Direction  of  the  Body  Axes  with  Respect  to  Wind  Axes  in 
Straight  Level  Perturbed  Condition 


Since 


0  =  v  +  « 


(2-125) 


(2-126) 


in  the  steady  situation 

0o  =  * 

where 

(X0  =  °  (2-127) 

In  fact,  in  order  to  maintain  straight  level  flight,  «<> 
has  to  be  small.  The  appropriate  relation  for  expressing  the 
angle  of  attack  <*  in  terms  of  the  w-velocity  if  obtained  by 
Equation  1-194  if  one  considers  small  <*  and  w ,  so  that  one 
can  express: 

a  =  w/()0  (2-128) 

From  Equation  2-127  then  we  deduce: 

=  *■/«.  =0  ^  W„  so  (2-129) 

Notice  that  by  this  alignment,  since  the  wind  axes  are 
fixed  for  fixed  relative  wind,  we  can  measure  the 
perturbation  quantities  about  in  the  wind  frame. 

It  is  further  desirable  to  express  the  Lateral 
Directional  equations  in  terras  of  the  sideship  angle  £ 
instead  of  the  v-  velocity.  The  appropriate  relation  for 
expressing  the  angle  /3  in  terms  of  the  v-  velocity  is 
obtained  by  Equation  1-195,  if  one  considers  small  f?  and  v  , 
so  that  one  can  express 

£  *  V/U0  (2-13  0) 

So  introducing 

from  Equations  2-126  and  2-127 

0,  =  lo  (2-131) 
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from  Equation  2-129 


- 

Wo  =  o 

(2-132) 

from  Equation  2-130 

fiu0 

(2-133) 

we  obtain 

Longitudinal  Equations 

i  *  =  0  +  8costf*e 

*  W  -  VoC^  +  <jsin^e  (2-134) 


— -  M  =  q, 

%  ° 

Lateral  Directional  Equations 

=  if  +  ^r  "  8°*#*  45 


_J _  .  _  *  Ivz 

i*,  L  ~  r 


(2-135) 


r  - 


Ixr 

Izz 


p 


with  auxilliary  relations: 

p  -  <h  -  ?  siu  y. 

^  *  9  (2-136) 

r  *  V  cos 

The  only  thing  that  remains  is  to  express  the  left-hand 
side  forces  and  moments  in  series  of  terms  involving 
aerodynamic  derivatives  and  small  perturbation  quantities, 
that  in  fact  generate  the  forces  and  moments. 
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2.7  EXPANSION  OF  FORCES  AND  MOMENTS 

In  Section  2.6,  it  was  assumed  that  the  motion  of  the 
vehicle  consists  of  small  deviations  around  a  reference 
steady  flight  path. 

Small  disturbance  theory  applied  on  the  perturbed 
equations  of  motion  assumes  that  the  aerodynamic  forces  and 
moments  are  functions  of  the  disturbance  velocities  (angular 
and  linear)  ,  the  control  surface  deflections  and  their 
derivatives. 

In  general,  for  a  function 

Kta,  Ji,  V  ,  ...  )  (2-137) 

the  Taylor  series  expansion  about  a  reference  value 

M3*,  )  (2-138) 

assuming  that  the  function  has  derivatives  of  all  orders 
defined,  is  written  as 


K.K..3K  + 

*(&1  w]  + 


(2-139) 


where  the  derivatives  must  be  evaluated  at  the  reference 


value  denoted  by  the  subscript  (o)  . 

Now,  making  the  assumption  that  the  deviations  S3, 

Sv  ,  are  so  small  that  product  terms  are  essentially  zero,  we 
obtain 


(2-140) 


92 


Using  this  method  and  provided  that  the  above  assumption 
is  valid,  we  succeed  to  express  the  perturbation  quantity  SK 
as  a  linear  combination  of  the  perturbation  parameters  from 
which  it  depends  upon. 

People  in  the  field  like  to  denote  the  partial 
derivatives  as  Kj)  ,  ,  R*  »  etc.,  instead  of  (SK/aa  )0  , 

(9t/3(i  )„  ,  (  9k /9v  )a  r  etc.  In  addition,  according  to  the 

convention  met  in  Section  2.6  of  dropping  the  S,  we  obtain: 

K=K{iafKpfi*K,v+  (2-141) 

where  ,  Kv ,  are  called  stability  derivatives. 

Fortunately,  not  all  the  stability  derivatives  in  an 
expansion  of  a  force  or  a  moment  have  a  value.  Some  of  them 
are  zero  or  essentially  zero. 

Expansion  then  of  forces  and  moments,  in  terms  of  the 
significant  stability  derivatives  yields: 

For  longitudinal  equations: 

4rX  =  Xuu  +  +  X5S  (2-142) 

~l  =  +  Zww  +Zi,  *  +  Z&S  (2-143) 

1 

2  M  *  Mu  u  +  f  Myi)  t  Mc^  +  Ms  3  ( 2-144) 

where  it  is  clear  that 

Y  1  SX  , 

n  %-  G  tc-  (2-145) 
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and  denotes  the  control  input  due  to  control  surfaces 
deflections. 

For  the  lateral -directional  equations: 


^  Y  *  Yv  v  +Yv*  +YpP  *Xr  vYsS  (2-146) 

J-  L  -  L,j(i  ♦  ♦  LpF  +Lrr  (2-147) 


N  *  ♦  Nj  p  f  Np  p  +  Mr  r  +  Ns£  (2-148) 

where  it  is  clear  that 

Y  j_  51  y  L 

v  =  m  'hv  >  C  =  «i  9f  (2-149) 

Therefore  the  equations  of  motion  become: 

Longitudinal  Equations 

■u  +  XyU-Xy^w  +  cjc0s^,0  -  XgS 

-Zyu  +.  w  -  2W  w  -  Zyo  ii  -Do<^  -Z<^  +  sin^.0  »  ZjS  (2-15  0) 

”  Mt>  D  ~  Myy  w  —  M  nj  w  4 

Lateral  Dire'  ..onal  Equations 

v  -  Yy  v  -  Yj,  v  -  Yp  p  *  l>„  r  -  Yr  r  -  ^  cos  -  Y§  <5 

~L^  “  4  £  +p  "  ^pP  -  Ix^x„  r  -  r  =  Ls5 

-  M(?(i  -  N(j  ^  f  -Npp  +  r-Nrr  =  N5  J  (2-151) 

with  auxilliary  relations  of  Equation  2-136. 

It  is  more  convenient  to  Laplace  Transform  the  equations 
and  arrange  them  in  matrix  form.  Further  augmenting  the 


'  .  >  •  •  i- 
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equations  by  the  auxilliary  relations,  using  Equation  2-130 
for  the  lateral  directional  set  and  dividing  the  Y-equation 
by  we  obtain  the  form 

Longitudinal  Equations 


5  -  Xo 

-x* 

0 

8C0S& 

V 

*6 

-Zu 

S -sZ ii  “Z*i 

asmfr 

W 

Zb 

-M, 

-sM*  -Mw 

s-Mi 

0 

% 

MS 

0 

0 

i 

-s 

9 

0 

(2-152) 


Lateral -Directional  Equations 


4-sYv  -Yv 

X 

*-Y r 

T 

-  £  0 

s-Lp 

~s  Ixx  ~Lr 

0  0 

-s  -  Njs 

-*£~N  P 

s-Mr 

0  0 

0 

i 

0 

-5  sfsin^1 

0 

0 

l 

O  -s(cos£o) 

where 


0 

Y/ 

P 

is 

r 

- 

NS 

£ 

0 

0 

:  (2-153) 


Y?*  =  Yf/U.  (2-154) 

Y?  =  Y-/U,  (2-155) 

Y£  =  Y5 1  Vo  (2-156) 

We  will  further  sirplify  the  equations  of  motion  in  the 
following  chapters  in  a  form  that  is  more  suitable  for 
obtaining  aerodynamic  transfer  functions,  and  computer 
solutions. 

Our  interest  first  will  be  focused  on  examining  the 


aerodynamic  stability  derivatives.  They  will  be  left  in 
their  dimensional  form  rather  than  non-dimensional. 
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CHAPTER  3 

NON-DIMENSIONAL  STABILITY  DERIVATIVES 

3.1  INTRODUCTION 

Stability  is  a  very  important  concept  in  successful 
vehicle  flight.  We  have  encountered  the  terms  "stability" 
and  "equilibrium"  but  we  did  not  emphasize  particularly  those 
concepts. 

In  the  last  chapter  we  mathematically  obtained  equations 
of  motion  and  linearized  them  about  a  referenced  straight 
flight  which  was  assumed  to  be  an  equilibrium  condition. 

But  what  is  equilibrium?  Equilibrium  is  a  state  of  rest 
or  uniform  (unaccelerated  motion)  in  a  straight  line  where 
all  forces  and  moments  acting  on  a  body  are  balanced  out. 

Stability  is  a  property  of  an  equilibrium  state.  We  will 
consider  two  types  of  stability  which  are  of  great  interest 
in  flight  vehicle  dynamics,  defined  as  follows:  Static 
stability  is  defined  by  the  direction  and  magnitude  of  the 
initial  tendency  to  return  to  equilibrium  condition  after  a 
small  disturbance  has  occurred.  It  is  distinguished  by 
positive  static  stability  if  the  body  tends  to  return  to  its 
equilibrium  state,  negative  static  stability  if  the  body 
diverges  from  its  equilibrium  state,  and  neutral  static 
stability  if  the  body  remains  in  the  disturbed  position. 


I  / 
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As  an  example/  the  pendulum  considered  in  Section  1.3  has 
a  positive  static  stability  at  its  position  but  has  a 
negative  static  stability  at  its  6*180*  position/  while  its 
motion  remains  invariant  if  the  pivot  is  moved  in  any  other 
position  and  is  said  to  have  a  neutral  static  stability. 

Dynamic  stability  determines  the  resulting  motion  in  time 
of  the  body  if  initially  disturbed.  It  is  distinguished 
similarly  in  positive,  negative,  and  neutral  dynamic 
stability. 

Various  conditions  of  dynamic  behavior  were  considered  in 
Chapter  1,  in  analyzing  spring  mass  damper  systems. 
Figure  3-1  summarizes  all  possible  time  responses  of  any 
dynamical  system. 

Dynamic  stability  implies  the  existence  of  static 
stability,  but  the  opposite  is  not  necessarily  true. 

Now  the  geometry  of  the  vehicle  comes  to  play  a  unique 
and  important  role  in  specifying  any  particular  mode  of 
dynamical  motion.  All  surfaces  exposed  to  the  airstream 
contribute  forces  and  moments  about  the  center  of  gravity, 
which  are  expressed  as  we  have  seen  in  the  last  chapter  as 
linear  combinations  of  stability  derivatives. 

The  physical  meaning  of  each  of  the  stability  derivatives 
is  straight  fo-ward.  Thus,  Xy  is  the  rate  of  chance  of  the 
force  generated  in  the  x-direction  due  to  a  v-velocity 
perturbation,  divided  by  the  mass  of  the  vehicle,  while  Xy  u 


cimat  ditflaanmtni 
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is  the  force  generated  in  the  x-direction  due  to  a  v-velocity 
perturbation  divided  by  the  mass  of  the  vehicle. 


In  this  chapter  we  will  give  more  insight  in  evaluating 
the  stability  derivatives  in  terms  of  the  vehicle's  geometry 
and  airstream  properties. 

3.2  LONGITUDINAL  STABILITY  DERIVATIVES 

Longitudinal  stability  derivatives  are  encountered  in  the 
longitudinal  equations  and,  hance,  they  uniquely  specify  the 
X-force,  the  Z-force  and  the  M-pitching  moment. 

Dimensional  analysis  approach  indicates  that  we  can 
express  the  aerodynamic  forces  and  moments  in  terms  of  their 
corresponding  non-dimensional  coefficients  which  are  usually 
deduced  from  wind  tunnel  tests  of  models. 


Figure  3-2 

Forces  and  Moments  Considered  in  Longitudinal  Motion 


In  particular,  we  can  express  the  longitudinal  forces  and 


moments  as 


lift 

l  =  cl 

ipV'S 

(3-1) 

drag 

ipV‘S 

(3-2) 

pitching  moment 

M  =  Cm 

ipV-'Sc 

(3-3) 

where  CL  ,Cj>  ,CM  are  the  corresponding  lift,  drag,  and  pitching 
moment  coefficients,  p  the  air  density,  V  the  velocity  of  the 
vehicle,  S  is  a  reference  area  usually  taken  the  wing  area, 
and  c,  a  reference  length,  usually  taken  the  length  of  the 
mean  aerodynamic  chord. 

The  Longitudinal  Stability  Derivatives  can  be  evaluated 
aa  follows: 

1.  The  Xu  Derivative 

Is  the  rate  of  change  in  the  X-force  due  to  a 
y -velocity  perturbation  divided  by  the  mass  m. 


Y  i  W 

'  m  3u  {3_4) 

According  to  Figure  3-2,  the  X-force  is  expressed  as: 
X  =  L  snrt#  -  D  cos#  +  T  cos  <xT  (3-5) 

where  T  denotes  the  thrust. 

For  a  small,  «T  essentially  zero,  and  making  use  of 
Equation  3-2  we  obtain: 


(3-6) 
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The  quantity  Cj>u  is  essentially  zero  for  low  subsonic 
speed,  but  sometimes  near  the  critical  roach  number 
(0.8<M<1.0)  where  a  large  increase  in  drag  occurs,  Ct>u 
obtains  a  large  positive  value. 

2.  The  1m  Derivative 

It  is  the  rate  of  change  in  the  Z-force,  due  to  a 
u-velocity  perturbation  divided  by  the  mass  m. 


Zu 


L 

'Wl 


(3-14) 


According  to  Figure  3-2,  the  Z-force  is  expressed  as: 
Z  =  W cos9  -  Lcos«  -Dsina  -  Tsmar  (3-15) 

which  is  simplified  to: 

Z  =  W  cos ®  -  L  (3-16) 

Neither  W  or  0  depend  on  v - vel oc ity .  Since  by 

similar  approach  as  in  the  Xu  derivative,  we  obtain: 


where 


/"  Ve  '"clQ 
1,1  1  2"^ 


(3-17) 

(3-18) 


The  quantity  Ciu  is  essentially  zero  for  low  subsonic 
speed,  but  sometimes  near  the  critical  mach  number  may  reach 
a  considerable  value  and  its  sign  can  be  changed  depending  on 
the  airframe  geometry,  mach  number,  dynamic  pressure  and 
aeroelastic  properties. 

For  subsonic  straight  level  flight  (#,=  0) 
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(3-19) 


£  W 

L  =  pV/S/2 

and  Zv  is  expressed  simply  as 

2  =  - 

Lv  Ve  (3-20) 

3.  The  Mu  Derivative 

It  is  the  rate  of  change  in  the  pitching  moment  M  due 
to  a  t»-velocity  perturbation,  divided  by  the  moment  of 
inertia. 


s-  It  OCm 

LMv  =  7  "37  (3-23) 

The  quantity  Cmw  is  very  sensitive  upon  the  variation 
of  the  u -velocity  and  its  sign  can  change  depending  upon  the 
geometry,  aeroelastic  properties,  raach  number  and  dynamic 
pressure.  Thrust  can  also  affect  for  propeller  driven 

aircraft. 

4.  The  Xw  Derivative 

It  is  the  rate  of  change  in  X-force,  due  to  a 
w-velocity  perturbation  divided  by  the  mass  m. 
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A  w-velocity  perturbation  causes  a  small  change  in 
the  angle  of  attack  « .  Thus,  in  this  case  we  cannot  neglect 
the  angle  of  attach  ot  which  is  approximately  given  by: 


w  =  — 

U, 


Therefore,  the  X-force  is  expressed  as 

X=  L  smoi  -Dcos<*  +T 


(3-25) 

(3-26) 


*  *  l  (CL  ~  Cv  LOS  a)  +T  (3-27) 

Neglecting  the  effect  of  the  w-velocity  on  the  thrust 
T,  by  differentiating  Equation  3-27  with  respect  to  w  and 
dividing  by  the  mass  m,  we  obtain 

v  J  pUo  s  ,  _  \ 

+CS  sma  -Cdetcos«)  (3_28) 

Now  assuming  *  small  we  express 


(3-29) 


where 

r  -  ^ 

'  07  (3-30) 

The  quantity  Cp^  represents  the  change  in  drag 
coefficient  due  to  angle  of  attack  changes.  When  the  angle 
of  attack  increases,  the  total  drag  increases  and  this 


quantity  has  normally  a  positive  value.  It  is  affected 
primarily  by  the  wings,  and  a  little  by  the  tail  and  the 
fuselage.  In  aircraft  dynamics  it  is  not  very  important 
parameter  mainly  because  the  changes  in  angle  of  attack  are 
small  and,  hence,  its  value  is  very  small. 

5.  The  Zw  Derivative 

It  is  the  rate  of  change  in  Z-force  due  to  a 
w-velocity  perturbation  divided  by  the  mass  m. 


7  =  1 

(3-31) 

The  Z-force  simillarly  is 

Z  =  W  cos  0  -  L  cos  ot  ~D  smoi 

(3-32) 

or 

Z  =  W  cos  9  -  i  pV2S  (Q  cos*  -Cj,  sw a) 

(3-33) 

Neither  W  or  0  depend  on  w-velocity, 

since  by 

differentiating  Equation  3-33  with  respect  to  w  and 

dividing 

by  the 

mass  m  we  obtain: 

7  <  2Z  f 

l«  ^  97  "  "  C-Q  since  .C^coscy  ♦  C*  cos*  ^  s.not) 

(3-34) 

Now,  assuming  *  small  we  express 

Zw  — 

(3-35) 

where 

c  -  ^ 

(3-36) 

105 


The  quantity  C  Ltl  represents  the  change  in  CL 
coefficient  due  to  angle  of  attack  changes.  Graphically,  it 
is  well-known  as  "lift  curve  slope",  as  shown  in  Figure  3-3. 


Figure  3-3 

Lift  Characteristic  Curves  for  Different  AR's 

The  changes  in  G  ,  are  linear  functions  of  ot  ,  below 
the  stall.  Aspect  ratio  (A«),  is  defined  simply  as  the  ratio 
(span/chord)  or  in  more  general  case  as: 

A2  _  CsiwO2  _ 

Cwiwjarea)  "  £  (3-37) 

This  ratio,  determines  the  slope  of  the  lift  curve.  As  shown 
in  Figure  3-3,  low  A^'s  produce  smaller  lift  coefficient  G 
for  the  same  angle  of  attack  <* . 

The  0*  derivative  is  always  positive  in  the  linear 
region.  It  is  affected  primarily  by  the  wings,  while  the 
tail  and  fuselage  contribute  a  little  in  aircrafts,  but 
largely  in  missiles.  Aeroelastic  effects  in  addition  may 
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cause  bending  or  twisting  of  the  wings  that  change  the  angle 
of  attack  and  hence  the  value  of  CLgl  . 

6.  The  Mw  Derivative 

Is  the  rate  of  change  in  the  pitching  moment  Mf  due 
to  a  w-velocity  perturbation,  divided  by  the  Iyy  moment  of 
inertia. 


M*  ■  is 


(3-38) 


A  w-velocity  perturbation  is  associated  with  an  angle 


of  attack  which  as  we  have  seen  is  given  by 

^  */l). 


(3-39) 


When  the  angle  of  attack  changes  from  its  equilibrium 
position,  the  lift  of  the  wing  changes  in  the  same  sense,  and 
that  produces  a  positive  or  negative  pitching  moment 
depending  on  the  center  of  gravity  location.  In  addition, 
the  lift  of  the  tail  changes  and  that  contributes  to  the 
total  pitching  moment.  These  contributions  together  with  the 
fuselage  contribution  make  up  this  derivative. 

We  desire  a  change  in  angle  of  attack  to  produce  a 
pitching  moment  that  tends  to  restore  the  original  angle  of 
attack.  This  requires  a  zero  pitching  moment  in  the 
undisturbed  reference  flight,  a  condition  known  as  "trim 
condition" . 

The  pitching  moment,  as  we  have  seen,  is  expressed  in 
terms  of  the  non-dimentional  pitching  moment  coefficient  as 
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(3-40) 


M  *  ipV*ScCH 

Because  a  small  change  in  w-velocity  doesn't  affect 
the  V  velocity,  by  differentiating  Equation  3-40  with  respect 
to  w  we  obtain: 


_  J_  P\J*£C 

Ci.i  Ov  O./  -ii  Ow  01 1  ^  Cl-/ 


(3-42) 


Ow  Qa  *2*/  ~  ’5*"  "  2Uo  '  2ot  (3-41) 

and  finally  by  letting  V-V<,  we  get 

p  S  Vo  c 

=  C“*  (3-42) 

The  center  of  gravity  location  plays  a  very  important 
role  in  derivative.  In  fact  Cu„  is  proportional  to  the 

distance  between  the  center  of  gravity  and  the  aerodynamic 
center,  the  center  where  the  lift  force  acts. 

If  the  center  of  gravity  coincides  with  the 
aerodynamic  center.  Cm#  is  zero;  if  it  is  ahead  of  the 
aerodynamic  center.  Cm*  is  negative  and  the  airframe  is  said 
to  be  statically  stable;  and  if  it  is  aft  of  the  aerodynamic 
center,  Cm„,  is  positive  and  the  airframe  is  said  to  be 
statically  unstable.  The  more  ahead  the  center  of  gravity  is 
located,  the  steeper  the  C"*  curve  becomes.  All  the  above 
cases  are  shown  in  Figure  3-4. 

It  is  clear  that  a  stable  situation  like  the  center 
of  gravity  at  B  or  A,  an  increase  in  angle  of  attack  from  <*» 
to  <*  produces  a  negative  pitching  moment  as  Figure  3-4  shows, 
which  tends  to  reduce  the  angle  of  attack. 


m 


Figure  3-4 

c..  for  Different  Center  of  Gravity  Locations 


Figure  3-4  also  shows  that  for  stable  situations,  for 
nearly  small  angles  of  attack,  where  the  lift  is  nearly  zero 
a  mechanism  must  be  provided  to  generate  a  positive  cm,  . 
Here  comes  the  effect  of  the  tail  which,  being  a  negative 
angle  of  attack,  generates  a  small  lift  force  which  though 
multiplied  by  a  larger  lever  arm  produces  a  balancing 
positive  pitching  moment. 

In  other  designs  we  find  the  tail  being  ahead  of 
wings  with  positive  angle  of  attack.  This  is  a  canard 
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configuration  most  usually  found  in  missiles ,  also  being  the 
mechanism  to  generate  a  positive  . 

In  Figure  3-4  curves  are  shown  as  linear 

functions  of  v.  This  is  not  true  for  large  variations  in  cxr 
but  can  be  nearly  true  for  small  variations  around  . 

7.  The  Derivative 

Is  the  rate  of  change  in  Z-force,  due  to  a  rate  of 
change  in  w-velocity  perturbation  divided  by  the  mass  m 


7  1  1L 

'  m  ‘  9w  (3-43) 

The  relation  o w/Uo,  leads  us  to  the  expression 

1L  - 

=  Do  Sir  (3-44) 

According  to  Equation  3-16  the  Z-force  is  simplified 
to 

Z  -  W  cos  0  -  L  (3-45) 

Neither  W  or  0  depend  on  rate  of  change  in  angle  of 
attack  of  ,  since  differentiating  Equation  3-45  with  respect  to 
o'  we  obtain 

%  _  _2L 

Qi  (3—46) 

where  from  Equation  3-1 


7  ■  ^  _  (  3L  -  i  i  ui ^  '^~L 

**  in  /ml),  9 i  /wit  2  P  ^ of  (3—48) 
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To  form  a  non-dimensional  coefficient,  we  multiply 


distribution  on  the  wing,  tail  and  fuselage,  does  not  adjust 
itself  instantaneously  to  its  equilibrium  value.  This,  a 
change  in  lift  and  drag  occurs  to  the  sudden  changes  in  angle 
of  attack. 


Its  sign  is  positive  for  low  speeds  but  can  be 
positive  or  negative  on  high  speeds  depending  on  aeroelastic 
effects. 


Also,  a  sudden  change  in  angle  of  attack  alters  the 
downwash  field  as  seen  by  the  tail  and  a  time  delay  is 
elapsed  before  the  tail  senses  this  alteration. 

All  those  effects  are  combined  to  give  rise  in  the 
derivative. 


8.  The  Mw  Derivative 

Is  the  rate  of  change  in  the  pitching  moment  M,  due 
to  a  rate  of  change  in  w-velocity  perturbation  divided  by  the 
moment  of  inertia. 

‘  L,y  fa  (3-52) 

This  derivative  'Jfysw  can  be  expressed  as: 


"  V .  Qit  (3-53) 

According  to  Equation  3-3 


M  -  j  p  V  Sc  Gvi 


(3-54) 


Since  V  is  independent  of  the  rate  of  change  in  angle 
of  attack  *,  by  differentiating  Equation  3-54  with  respect  to 
v  we  obtain: 


„  ipV'Sc 

Z  V  Oi 


(3-55) 


Therefore, 


M*  .-L =  1 


ipV7Sc 


%  '3*  V.lw  UaIw  2  r  %  (3-56) 

To  form  a  non-dimensional  coefficient,  we  multiply 
and  divide  Equation  3-56  by  c/2Uo  so  that  letting  V  *  Uo,  we 
get 

Ws_L..ipU*Sc  -i-  - 

Vofyy  2  -?U>  ^(yc/2V.) 

pSc* 


M. 


(3-57) 


or 


M, 


4% 


Cm, 


(3-58) 
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where 


r 

^  ‘-9CK/2U.)  (3-59) 

For  quantity  Cm*  represents  a  change  in  pitching 
moment  coefficient  due  to  a  variation  in  the  rate  of  change 
of  the  angle  of  attack  «. 

The  change  in  lift  that  is  produced  by  sudden  changes 
in  angle  of  attack  discussed  already,  multiplied  by 
corresponding  lever  arms,  generates  changes  in  pitching 
moment  and  that  is  why  Cm*  arises. 

Same  lag  effects  in  the  downwash  field  discussed, 
contribute  to  this  derivative,  as  well  as  aeroelastic 
effects,  that  also  disturb  the  air  stream  of  the  wing  and 
tail,  so  that  the  pressure  distribution  over  the  aerodynamic 
surfaces  does  not  adjust  instantaneously. 

The  sign  of  this  derivative  is  negative  and  is 
relatively  important  in  longitudinal  dynamics.  Large 
absolute  values  of  this  derivative  are  desired  because  it 
causes  damping. 

9.  The  Zq  Derivative 

Is  the  rate  of  change  in  Z-force,  due  to  a  q-pitching 
velocity  perturbation  divided  by  the  mass  ra 


to: 


7  i  2L 

\  (3-60) 

According  to  Equation  3-16,  the  Z-force  is  simplified 


£  =  W  cos  9  -  L  (3-61) 
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Neither  W  or  9  depend  on  the  pitching  velocity  (q) , 
since  differentiating  Equation  3-61  with  respect  to  q,  we 
obtain: 


(3-62) 


where  from  Equation  3-1 


%  2  p  “H 


(3-63) 


Therefore: 


7  -  —  ^  _  ■<  2k  -  _  i.  1  flV/2c 
H  *  'k  "*•  3q  '*«>  2^  ' 


^  ^  2r*  ~  3%  (3-64) 

To  form  a  non-dimensional  coefficient  we  multiply  and 
divide  Equation  3-64  by  c/2Uo  so  that  letting  V=Uo  we  get: 


where 


■  i  ipUl$  - —L— 

*  2  "  2Vo  /2Va 

Za  «.  c 


c  -  r^~L 

L%  =  gC^c/2U.) 


(3-65) 


(3-66) 


(3-67) 


The  quantity  0.^  represents  a  change  in  lift 
coefficient  due  to  a  variation  in  the  pitching  velocity. 

If  the  airframe  pitches  about  its  center  of  gravity , 
the  angle  of  attack  changes  and  a  change  in  left  of  wing, 
tail  and  fuselage  is  developed  that  gives  rise  to  this 
derivative. 


If  has  a  small  positive  value  in  low  speed  flight  but 
may  have  large  values  positive  or  negative  in  high  speeds  due 
to  aeroelastic  effects. 

Curved  longitudinal  flight  paths  give  rise  also  to 
centrifugal  forces  on  all  components  of  the  airframe  that  due 
to  aeroelastic  effects  cause  a  bending  of  aerodynamic 
surfaces  and,  hence,  a  change  in  lift. 

In  the  past  without  computer  solutions  available, 
this  derivative  was  considered  zero  to  simplify  the  analysis. 

10.  The  Mq  Derivative 

Is  the  rate  of  change  in  the  pitching  moment  H,  due 
to  a  pitching  velocity  perturbation,  divided  by  the  Iyy 
moment  of  inertia. 


M.  .i-  .21 

* 

(3-68) 

According  to  Equation  3-3, 

M  =  ipV2ScCM 

(3-69) 

Since  V  is  independent  of  the  pitching  velocity  q,  by 
differentiating  Equation  3-69  with  respect  to  q,  we  obtain: 


.  ipv2s  c  ^ 
H  2  p  H 


Therefore, 


M, -f-  lr  =  -4- 


(3-70) 


(3-71) 
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To  form  a  non-dimensional  coefficient,  we  multiply 
and  divide  Equation  3-71  by  c/2Uo,  so  that  letting  V=Uo  we 
get: 


M.-  4"  p  U>2  S  c  -£-■ 
%  a  i-  r  011- 


(3-72) 


where 


pV»$  c2 

=  "Til-  Cmi 


=  ox:^c/2ti) 


(3-73) 


(3-74) 


The  quantity  represents  a  change  in  pitching 
moment  coefficient  due  to  a  variation  in  the  pitching 
velocity  q. 

The  change  in  lift  that  was  produced,  due  to  pitching 
velocity  perturbations  discussed  already,  multiplied  by 
corresponding  lever  arms,  generates  changes  in  pitching 
moment  and  that  is  why  C  rises. 

The  tail  is  the  main  contributor  to  this  derivative 
which  is  also  referred  to  as  "pitching  damping  derivative". 
Aeroelastic  effects  and  centrigugal  forces  generated  due  to 
curved  longitudinal  flight  paths  also  contribute  to  this 
derivative. 

In  longitudinal  dynamics  is  very  important.  It 
has  a  negative  value,  and  the  larger  it  is,  the  more  the 
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damping  is  increased;  a  condition  that  is  highly  desired.  In 
high  speed  flights  may  have  positive  or  negative  value 
depending  on  aeroelastic  effects.  It  is  also  noted  that 
tailless  airplanes  have  a  poor  damping. 

3.3  LATERAL  -  DIRECTIONAL  STABILITY  DERIVATIVES 

Lateral -Directional  stability  derivatives  are  encountered 
in  the  lateral  directional  equations  andr  hence,  they 
uniquely  specify  the  Y-force,  the  L-rolling  moment  and  the 
N-yawing  moment. 

Dimensional  analysis  approach  indicates  that  we  can 
express  the  lateral-directional  aerodynamic  forces  and 
moments  in  terms  of  their  corresponding  non-dimensional 
coefficients  as: 


Y-force 

Y 

,C^pV25 

(3-75) 

L-rolling  moment 

l 

=  ce  J  pv*st 

(3-76) 

N-yawing  moment 

KJ 

=  C*  i  p\Z*Sb 

(3-77) 

Cy  t  ce  ,  c, 1  are 

the 

corresponding  Y-force, 

rolling 

moment  and  yawing  moment  coefficients,  and  b  is  a  reference 
length,  usually  taken  the  wing  span. 

The  Lateral-Directional  Stability  Derivatives  can  be 
evaluated  as  follows: 

1.  The  Yv  Derivative 

Is  the  rate  of  change  in  Y-force  due  to  a  v-velocity 
perturbation  divided  by  the  mass  m: 
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(3-78) 


Y  _  4  3Y 
v  ff*  9v 

A  v-velocity  perturbation  causes  a  small  change  in 
the  angle  of  sideslip  j2  which  is  expressed  according  to 
Equation  2-133,  approximately  as: 

p  s  */U>  (3-79) 

so  that 


Y  .  j_  2L 

v  ‘  *0.  '  (3-80) 

Assuming  that  the  V-velocity  is  independent  of  the 
angle  of  sideslip  fi  ,  differentiating  Equation  3-75  with 
respect  to  fi ,  we  obtain: 


21 

V 


.  i  «»!s  2E» 

2  *  ty 


Therefore,  letting  V=Uo  we  get: 


(3-81) 


Y  PUoS  r 

Tv  ‘  2m  (3-82) 

where 

C9/i  =  ^V9/S  (3-83) 

The  quantity  represents  the  change  in  the  Cy  side 
force  coefficient  due  to  a  change  in  the  angle  of  sideslip  fi. 

The  main  contributor  of  this  derivative  is  the 
vertical  tail.  Secondarily  the  fuselage,  and  also  the  wings 
contribute  a  little. 
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Figure  3-5  shows  the  forces  that  the  airframe 
experiences  due  to  a  v-velocity  perturbation.  First  an  angle 
of  sideslip  builds  up  and  all  aerodynamic  surfaces  are 
sensing  a  change  in  the  approaching  free  stream. 
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Figure  3-5 

Forces  Experienced  Due  to  a  v-Velocity  Perturbation 

As  a  consequence,  a  force  Ff  due  to  the  fuselage  and 
a  force  Fvi  due  to  the  vertical  tail  appears.  The  wings  also 
experience  a  different  lift  force.  For  the  configuration 
shown  in  Figure  3-5,  the  right  wing  experiences  more  lift 
than  the  left  wing  because  it  senses  a  greater  V-normal 
component  of  the  approaching  free  stream.  Those  lift  forces 
may  contribute  to  a  small  change  in  the  Y-force. 

All  these  contributions  make  a  rise  to  derivative 
which  is  very  important  in  lateral -directional  dynamics.  It 
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introduces  a  damping  that  tends  to  eliminate  the  angle  A, 
because  the  resultant  Y-force  appears  acting  on  a  point  aft 
of  CG  location  and  that  generates  the  appropriate  moment  to 
do  this. 

Its  value  is  negative  and  it  might  be  reasonable  to 
desire  a  large  value,  for  the  £  angle  to  die  out  quickly. 
Although  high  maneuverability  airplanes  must  have  a  small 
negative  value  for  this  derivative,  if  they  were  to  perform 
easy  roll  maneuvers  and  banked  turns. 

2.  The  Lv  Derivative 

Is  the  rate  of  change  in  L-rolling  moment  due  to  a 
v-velocity  perturbation  divided  by  the  Ixx  moment  of  inertia. 


,  1  3L 

Lv  *  I**  9v 

It  is  also  clearly  expressed  as: 


(3-84) 


Lv  - 


1 


*31 


UT**  W  (3-85) 

Differentiating  Equation  3-76,  with  respect  to  we 


obtain 


V2S  b  — 

V  If 


Therefore,  letting  V=Uo  we  get: 


(3-86) 


where 


C  <^t 

^  (3-88) 

The  quantity  represents  the  change  in  the  Cg 

rolling  moment  coefficient  due  to  a  change  in  the  angle  of 
sideslip 

It  is  usually  referred  to  as  the  effective  Dihedral 
Derivative. 

Dihedral  is  often  used  as  a  mechanism  to  improve 
lateral-directional  stability.  Figure  3-6 (a)  shows  a  head-on 
view  of  an  airplane  whose  wings  are  tilted  up  at  some 
dihedral  angle  with  respect  to  the  horizontal. 

We  have  already  discussed  that  due  to  a  v-velocity 
perturbation,  according  to  the  configuration  shown  in 
Figure  3-5,  where  a  positive  sideslip  angle  (1  is  shown,  right 
wing  generates  more  lift  than  the  left.  This  is  true  for  the 
swept  wings.  If  the  wings  are  not  swept,  the  mechanism  to 
produce  the  same  effect  is  the  dihedral.  Usually  even  swept 
wings  are  given  a  small  dihedral  to  improve 
lateral-directional  stability. 

To  understand  the  dihedral  effect,  consider  the 
aircraft  shown  in  Figure  3-7  with  non-swept  wings  but  with 
some  dihedral  angle.  The  aircraft  is  distributed  by  an  angle 
of  sidesl.  p/2. 
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Figure  3-6 

Dihedral  Effect  on  Lateral -Directional  Stability 


Figure  3-7 

Dihedral  Effect  on  Non-Swept  Wings 
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The  perpendicular  components  of  the  free  stream 
velocity  are  equal  but  the  angle  of  attack  that  the  right 
wing  projects  to  this  perpendicular  component  is  greater  than 
the  angle  of  attack  that  the  left  wing  projects,  as  an  effect 
of  the  dihedral  angle. 

As  a  consequence,  the  right  wing  generates  more  lift 
than  the  left  wing  as  shown  in  Figure  3-6(c).  Those  lift 
forces  contribute  to  a  change  in  the  rolling  moment 
coefficient. 

Figure  3—6 ( d )  and  (c)  show  the  effect  of  the 
high-wing  placement  which  has  the  same  effect  as  the  dihedral 
angle.  It  is  clearly  shown  that  in  the  low  wing  placement 
configuration,  the  fuselage  does  not  permit  the  vertical 
component  of  the  free  stream  to  act  upon  the  right  wing  to 
produce  a  stabilizing  moment,  whereas  in  the  high  wing 
placement  configuration  the  fuselage  effect  is  nearly  zero 
and  that  makes  a  stabilizing  rolling  moment  to  show  up. 

In  addition,  the  force  acting  on  the  vertical  tail 
shown  in  Figure  3-5,  multiplied  by  the  corresponding  vertical 
arm  relative  to  the  x-longitudinal  axis,  contributes  also  to 
the  same  ef  feet . 

is  nearly  always  negative  in  sign  and  the 
corresponding  rolling  moment  L,  is  negative  for  a  positive 
sideslip  angle  .  It  is  a  very  important  derivative  in 
lateral  directional  dynamics  because  it  causes  a  damping  and 
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improves  the  stability.  Although  high  negative  values  of 
this  derivative  will  tend  to  decrease  the  maneuverability  of 
the  airplane  in  roll  movements.  In  missiles  a  low  negative 
value  of  this  derivative  is  desired. 

3.  The  Nv  Derivative 

Is  the  rate  of  change  in  N-yawing  moment  due  to  a 
v-velocity  perturbation  divided  by  the  Izz  moment  of  inertia. 


w  1  W 

N  y  »  —  ■  T7 — * 

Jn  9v 


(3-89) 


It  is  also  clearly  expressed  as: 


N  -  J- 
Ny  U.I„  W 


(3-90) 


obtain: 


Differentiating  Equation  3-77 ,  with  respect  to  ft. 


we 


,  l  pv'sb  i? 

2  r  9£ 


Therefore,  letting  V=(Jo  we  get: 


(3-91) 


p  u»  S  [> 

v  '  Tj IT  Ci,(*  (3-92) 

where 

Cvp  =.  (3-93) 

The  quantity  On*  represents  the  change  in  the  Cn 
yawing  moment  coefficient  due  to  a  change  in  the  angle  of 
sideslip  ft* 


124 


It  is  usually  referred  to  as  the  Weathercock 
Derivative.  The  main  contributor  of  this  derivative  is  the 
vertical  tail  and  secondarily,  the  fuselage.  If  the  airframe 
experiences  a  v-velocity  perturbation  and  an  angle  of 
sideslip  £  appears,  the  free  relative  stream  strikes  the 
fuselage  and  the  vertical  tail  in  such  a  manner  that  the 
force  Ff  ,  shown  in  Figure  3-5,  produces  a  destabilizing 
yawing  moment  and  the  force  a  stabilizing  moment.  Wings 
contribute  a  little  in  the  positive  (stabilizing)  sense.  The 
total  value  of  this  derivative  is  positive,  signifying  static 
directional  stability. 

The  derivative,  is  very  important  in 

lateral-directional  dynamics.  The  vertical  tail  must  be 
properly  designed  to  provide  a  sufficient  weathercock  effect. 
At  large,  sideslip  angles  a  stall  may  occur  due  to  flow 
separation  and  a  catastrophic  sideslip  divergence  may  result. 
A  low  aspect  ratio  prevents  stalling  but  decreases  stability. 

Usually  a  large  positive  value  of  the  Gy  derivative 
is  desired,  by  extending  the  vertical  tail.  A  large  value 
helps  the  pilot  in  good  coordinated  turns  and  prevents 
sideslip. 

4.  The  Yy  Derivative 

Is  the  rate  of  change  in  Y-force,  due  to  a  rate  of 
change  in  v-velocity  perturbation,  divided  by  the  mass  m. 

Y  31 

*  m  %  (3-94) 
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c  - 

H  ”  ^C^b/ZOo)  (3-99) 

The  quantity  represents  the  change  in  the 

coefficient  with  variation  in  rate  of  change  of  the  sideslip 
angle  /*. 

If  the  sideslip  angle  /£  changes  rapidly,  the  pressure 
distribution  on  the  vertical  tail  and  fuselage  does  not 
adjust  itself  instantaneously  to  its  equilibrium  value. 

This  effect,  as  well  as  aeroelastic  effects,  give 
rise  in  derivative  but  its  value  is  very  small  and  is 

usually  neglected. 
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5.  The  Lv  Derivative 


Is  the  rate  of  change  in  L-rolling  moment,  due  to  a 
rate  of  change  in  v-velocity  perturbation  divided  by  the  Ixx 
moment  of  inertia. 


,  .  (  9L 

*  'Ixx  <3v  (3-100) 

It  is  also  clearly  expressed  as: 

I  •  4  'SL 

v  =  aixx  ’  W  (3-101> 

Differentiating  Equation  3-76,  with  respect  to  /3 ,  we 

obtain: 


% 


To  form  a  non-dimensional  coefficient 


—  -  i  nV3S  h  ~  ^ 

0/S  2  ?  2 U.  Ufa /2V.) 

Therefore,  letting  V=Do 


(3-102) 


(3-103) 


v  '  4Jxx  e<i  (3-104) 

where 

r 

'  Q((3b/2U»)  (3-105) 

The  quantity  represents  the  change  in  the  Cp 

rolling  moment  coefficient  with  variation  in  rate  of  change 
of  the  sideslip  angle  {!• 


The  change  in  the  side  forces  on  the  vertical  tail 
and  the  fuselage  produced  due  to  sudden  changes  in  the 
sideslip  angle  ft,  discussed  already/  multiplied  by 
corresponding  lever  arms  provide  a  total  change  in  the  Q 
rolling  moment  coefficient. 

This  derivative  has  a  very  small  value  because  the 
lever  arms  are  very  small  and  it  is  usually  neglected. 

6.  The  Ny  Derivative 

Is  the  rate  of  change  in  the  N-yawing  moment/  due  to 
a  rate  of  change  in  v-velocity  perturbation  divided  by  the 
Izz  moment  of  inertia. 


4  3N 


la  Bv 

It  is  also  clearly  expressed  as: 


(3-106) 


1 


9N 


"  V^z  3(*  (3-107) 

Differentiating  Equation  3-77,  with  respect  to  /3,  we 


obtain: 


— r  =  -  pv  S  b  7 rr- 

2f>  fy 

To  form  a  non-dimensional  coefficient: 


(3-108) 


9N 


2(1  2  r  oil 


9Cn 


2Ve  /3C<ib/2U.) 

Therefore,  letting  V»Uo  we  get: 

Hi  .  if*  a, 

A  Im  <* 


(3-109) 


(3-110) 
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where 


»  *  9C/ib/2a)  (3-111) 

The  quantity  represents  the  change  in  the  Gi 

yawing  moment  coefficient  with  variation  in  rate  of  change  of 
the  sideslip  angle  0. 

The  change  in  the  side  forces  on  the  vertical  tail 
and  the  fuselage  produced  due  to  sudden  changes  in  the 
sideslip  angle  /2  ,  discussed  already#  multiplied  by 
corresponding  lever  arms  provide  a  total  change  in  the  G 
yawing  moment  coefficient. 

The  lever  arm  of  the  vertical  tail  this  time  is  long 
and  hence#  this  derivative  may  be  considerable.  Positive 
values  of  this  derivative  increase  the  damping  in  the  yawing 
motion  and  are  desirable  in  general. 

7.  The  Yp  Derivative 

Is  the  rate  of  change  in  the  Y-force#  due  to  a 
p-rolling  velocity  perturbation  divided  by  the  mass  m. 

Y  _  i.  f>! 

]P  "  m  9p  (3-112) 

Differentiating  Equation  3-75#  with  respect  to  p#  we 

obtain: 


3Y 


(3-113) 
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To  form  a  non-dimensional  coefficient 

H,  -  L Jl 

9p  "  2  P  2Uo  3Cpb/2Lt )  (3-114) 

Therefore,  letting  V=Uo 

Y  p  S  t> 

Yp  =  M/P  (3-115) 

where 

C  ^ 

'  VCpb/2U„)  (3-116) 

The  quantity  represents  the  change  in  the  Cy 

Y-force  coefficient  with  variation  in  the  rolling  velocity  p. 

When  an  aircraft  has  a  roll  velocity  p,  a  force  is 
acting  on  the  vertical  tail  that  usually  opposes  the  roll 
movement.  This  consists  the  main  contribution  to  this 
derivative  and  is  small  enough  to  be  neglected. 

The  force  acting  on  the  vertical  fin  may  be  affected 
by  the  sidewash  from  the  wing  that  gives  a  slightly  different 
direction  of  the  airstream. 

8.  The  Lp  Derivative 

Is  the  rate  of  change  in  the  L- rolling  moment,  due  to 
a  rolling  velocity  p  perturbation  divided  by  the  Ixx  moment 
of  inertia. 


■  _  4  9L 

P  Jxx  9p 


(3-117) 
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Differentiating  Equation  3-76,  with  respect  to  p,  we 

obtain: 


OL 

Qp 


To  form  a  non-dimensional  coefficient. 


(3-118) 


'cJL  i  ,  b  QQ 

Op"  =  2  P  b  2Uo  9fpb/2U0  (3-119) 

Therefore,  letting  V»CJo 

pUoSb- 

p  '  4Ixx  ?P  (3-120) 

where 

C,  .  % 

P  3  O(pb/2UO  (3-121) 

The  quantity  Qp  represents  the  change  in  Q  rolling 
moment  coefficient  with  variation  in  the  rolling  velocity  p. 

This  derivative  is  also  referred  to  as  "roll  damping 
derivative",  because  it  has  usually  a  negative  or  damping 
effect  coming  mainly  from  the  wings  and  secondarily,  from  the 
vertical  and  horizontal  tails. 

When  the  aircraft  rolls,  the  force  acting  on  the 
vertical  tail,  shown  in  Figure  3-8,  multiplied  by  the 
corresponding  lever  arm,  produces  an  opposing  rolling  moment. 
In  addition,  on  the  down-going  wing,  the  effective  angle  of 
attack  is  increased  and  on  the  up-going  wing,  decreased. 
This  makes  the  lift  on  the  down-going  wing  greater  and  thus 
an  opposing  rolling  moment  is  produced.  Similar  effect,  on  a 
smaller  scale,  is  produced  by  the  horizontal  tails. 
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Figure  3-8 

Forces  Due  to  Rolling  Velocity  (p) 

This  is  true  as  long  as  the  down-going  wing  is  not 
flying  near  the  stall  due  to  a  large  angle  of  attack 
experienced  in  some  cases.  In  such  cases  a  positive  effect 
will  result  in  autorotation. 

The  derivative  is  quite  important  in 

lateral-directional  dynamics.  The  value  of  this  derivative 
directly  affects  the  ailerons  design.  A  small  negative  value 
of  this  derivative  will  result  in  a  better  aileron  input 
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response  while  a  large  negative  value  will  make  the  airframe 
heavily  damped  in  roll  maneuvers. 

9.  The  Np  Derivative 

Is  the  rate  of  change  in  the  N-yawing  moment,  due  to 
a  rolling  velocity  9  perturbation,  divided  by  the  Izz  moment 
of  inertia. 

M  _  J_  2! 

P  r  1*Z  '  9p  (3-122) 

Differentiating  Equation  3-77,  with  respect  to  p,  we 

obtain: 

=  i  pV^Sb 

9p  2  r  9p  (3-123) 

To  form  a  non-dimensional  coefficient, 

9?  *  2  P  P  sq,  2(?bl2U,l  (3-124) 

Therefore,  letting  V=Uo 

N  -  f”tSb’  c 

Np  '  C”f  (3-125) 

where 


r  _ 

*P  *  9Cpb/2u.)  (3-126) 

The  quantity  Cnp  represents  the  change  in  C*  yawing 
moment  coefficient  with  variation  in  the  foiling  velocity  p. 

This  derivative  usually  has  a  negative  or  positive 
effect  coming  from  the  wings,  the  vertical  tail  and 
secondarily,  from  the  horizontal  tail. 


> 
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When  the  aircraft  rolls,  the  force  acting  on  the 
vertical  tail,  shown  in  Figure  3-8,  multiplied  by  the 
corresponding  lever  arm  produces  a  yawing  moment  towards  the 
same  side  of  the  roll  motion,  i.e.,  roll  to  starboard  causes 
yaw  to  starboard.  This  contributes  a  positive  effect  and 
depands  upon  the  vertical  tail  geometry  and  the  sidewash  from 
the  wings. 

In  addition,  the  contribution  of  the  wings  can  be 
considered  to  be  made  of  two  parts,  one  opposite  the  other, 
the  total  effect  being  negative.  When  the  airframe  rolls,  on 
the  down-going  wing,  the  effective  angle  of  attack  increases 
and  on  the  up-going  wing  decreases,  as  shown  in  Figure  3-8. 
As  a  consequence,  the  lift  and  drag  are  greater  on  the 
down-going  wing  than  on  the  up-going  wing.  The  projections 
of  the  drag  forces  on  the  xcy  plane  usually  produce  a  yawing 
moment  towards  the  same  side  of  the  roll  motion.  This 
produces  a  positive  effect  while  the  projections  of  the  lift 
forces  produce  a  yawing  moment  towards  the  opposite  side  of 
the  roll  motion,  and  this  produces  a  negative  or  damping 
effect. 

The  overall  value  of  this  derivative  is  usually 
negative,  i.e.  roll  to  starboard  produces  a  yaw  to  port.  At 
higher  speeds  the  numerical  value  decreases  and  may  become 
positive . 
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This  derivative  is  fairly  important  in 
lateral-directional  dynamics,  as  far  as  roll  maneuvers  and 


their  effect  to  yaw. 

10.  The  Yr  Derivative 

Is  the  rate  of  change  in  the  Y-side  force,  due  to  a 
yawing  velocity  perturbation  r,  divided  by  the  mass  m. 


Y  -  L  JL 

lr  =  m  Qr  (3-127) 

Differentiating  Equation  3-75,  with  respect  to  r,  we 

obtain: 


2S..i,V‘s  3l 

>  2  r  ° 

To  form  a  non-dimensional  coefficient. 


_  i  0  y2s  JL  ^ 

Or  “  2  r  2Vo  OCrb /2V.) 

Therefore,  letting  V=0o 


where 


Yr 


pUoSb 

At* 


C\*r 


(3-128) 


(3-129) 


(3-120) 


'  Xrb/20,)  (3-131) 

The  quantity  Cyr  represents  the  change  in  Cy  Y-force 
coefficient  with  variation  in  the  yawing  velocity  r. 

When  the  aircraft  has  a  yaw  velocity  r,  as  Figure  3-9 
shows,  a  force  is  acting  on  the  fuselage  that  usually  has  a 
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negative  or  damping  effect  on  the  yaw  motion  and  a  force  is 
acting  on  the  vertical  tail  that  usually  has  the  same  effect. 


Figure  3-9 

Forces  Due  to  Yaw  Velocity  (r) 

Both  contributions  are  small  enough  to  be  neglected, 
therefore,  this  derivative  is  of  little  importance  in 
lateral-directional  dynamics. 

11.  The  Lr  Derivative 

Is  the  rate  of  change  in  the  L-rolling  moment,  due  to 
a  yawing  velocity  r  perturbation,  divided  by  the  Ixx  moment 
of  inertia. 

I  _L  JL 

L,r  s  I*,  (3-132) 
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Differentiating  Equation  3-76 ,  with  respect  to  r,  we 


obtain: 


—  =  - 

9 r  2  '  9r 

To  form  a  non-dimensional  coefficient. 


2L  i  v^u,  _L 

9r  *  2  P  bb  2Uo  9C^b/ 2Ub) 

Therefore,  letting  V=Uo 


I  pu=st!  , 

ir  =  ~ IT  C(r 


where 


(3-133) 


(3-134) 


(3-135) 


,  .  *1 

'  ^(rb/2 U.)  (3-136) 

The  quantity  Qr  represents  the  change  in  Q  rolling 
moment  coefficient  with  a  variation  in  the  yawing  velocity  r. 

When  the  aircraft  is  yawing  about  its  vertical  axis, 
the  outer  wing  moves  faster  than  the  inner.  As  a 
consequence,  the  lift  on  the  outer  wing  is  greater  than  that 
of  the  inner  and,  hence,  a  rolling  moment  is  produced  towards 
the  same  side  of  the  yaw  motion.  This  contributes  a  positive 
effect,  i.e.  yaw  to  starboard  produces  roll  to  starboard. 

In  addition  to  this  major  wing  contribution,  the 
force  acting  on  the  vertical  tail  multiplied  by  the 
corresponding  lever  arm,  produces  the  same  rolling  effect. 
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The  total  value  of  this  derivative  is  usually 
positive,  but  it  is  not  so  important  in  lateral-directional 
dynamics. 

12.  The  Nr  Derivative 

Is  the  rate  of  change  in  the  N-yawing  moment,  due  to 
a  yawing  velocity  r  perturbation,  divided  by  the  Izz  moment 


of  inertia. 


Nr  .  -L  .  *L 
Izz  Or 


(3-137) 


Differentiating  Equation  3-77,  with  respect  to  r,  we 


obtain: 


DM  i  ,  or* 

FT  =  7  Pv  7— 

Qr  2  r  9r 

To  form  a  non-dimensional  coefficient, 

^  .  i  pv'sb  i  _ 

Or  2  r  20,  V.rb/20.) 

Therefore,  letting  V*Do 


(3-138) 


(3-139) 


(3-140) 


where 


r  (K1\ 

'BCrb/20.)  (3-141) 

The  quantity  Qnr  represents  the  change  in  G,  yawing 
moment  coefficient,  with  a  variation  in  the  yawing 
velocity  r. 


138 


This  derivative  is  known  as  the  "yaw  damping 
derivative",  because  it  usually  has  a  negative  effect. 

The  main  contributor  is  the  vertical  tail.  When  the 
aircraft  is  yawing,  the  force  acting  on  the  vertical  tail 
multiplied  by  the  corresponding  long  lever  arm,  produces  a 
considerable  yawing  moment  that  opposes  the  initial  yawing 
disturbance . 

In  addition,  the  drag  on  the  outer  wing  is  greater 
than  that  of  the  inner  wing.  That  produces  a  yawing  moment 
that  opposes  the  initial  yawing  disturbance.  The  force  on 
the  fuselage  also  contributes  an  opposing  effect. 

This  derivative  is  very  important  in 
lateral-directional  dynamics.  Its  total  effect  is  negative 
and  large  values  are  desired  for  more  effective  damping. 

3.4  CONTROL  SURFACE  DEFLECTIONS 

We  have  discussed  so  far  all  the  important  stability 
derivatives  that  rise  due  to  the  airframe  geometry.  They 
uniquely  specify  the  forces  and  moments  about  the  center  of 
gravity  provided  that  they  are  multiplied  by  the  correct 
factor  that  depends  upon  the  airstream  and  airframe 
characteristics. 

A  stability  analysis  is  possible  now  for  determining  the 
longitudinal  and  lateral  directional  modes  of  motion, 
corresponding  to  an  initial  disturbance,  i.e.  homogeneous 
solution  to  the  differential  equations. 
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In  addition^  some  means  of  control  in  longitudinal  and 
lateral  directional  motion  are  available  to  the  human  or 
automatic  pilot.  They  are  shown  in  Figure  3-10. 
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Figure  3-10 

Basic  Longitudinal  and  Lateral  Control  Systems 


Elevator,  controls  the  longitudinal  motion  of  the 
airframe  providing  nose  up  or  nose  down  moments  as  a  response 
to  corresponding  control  surface  deflections. 

Aileron  and  rudder,  control  the  lateral  directional 
motion  providing  rolling  and  yawing  moments  correspondingly 
as  a  response  to  the  corresponding  control  surfaces 
deflections. 

The  symbols  used  for  the  corresponding  angles  of  rotation 
and  control  surface  deflections  are  shown  below: 
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In  the  longitudinal  and  lateral  directional  equations  of 
motion.  Equations  2-152  and  2-153,  respectively,  we  have 
denoted  by  the  symbol  $  the  control  input  due  to  control 
surfaces  deflections. 

The  right-hand  side  of  the  longitudinal  Equation  2-152  is 
expanded  as  follows: 

=  xn 

Zs  =  Ztj  (3-142) 

=  Mri  +  Mirt 
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and  the  right-hand  side  of  the  lateral  directional 
Equation  2-153  as: 

YS  =  YS 

4  =  Lg  * Lg  * L5  (3-143) 

Ns  =  Ng  *  N5  +  Nj 

In  addition,  to  the  right-hand  side  terms  that  take  care 
of  the  main  control  surfaces,  there  may  exist  terms  rising 
from  the  auxiliary  surfaces  such  as  flaps,  stabilizers,  dive 
brakes,  etc. 

A  general  approach  suggests  to  examine  the  subscripted 
derivatives  instead  of  the  particular  derivatives  dedicated 
for  a  specific  control  surface  deflection,  i.e.  the  total 
effect  by  all  the  possible  combinations  of  the  control 
surface  deflections.  The  approach  is  similar  to  that 
followed  in  Sections  3.2  and  3.3. 

1.  The  *S  Derivative 

Is  the  rate  of  change  in  the  x-force  due  to  a 
control  surface  perturbation  divided  by  the  mass  m. 


y«  =  _  _ 

5  m  'dS  (3-144) 

By  &  we  denote  a  general  control  surface  deflection. 
For  longitudinal  control  it  might  be  elevator,  stabilizer, 
flap,  3lat,  dive  brakes,  etc. 

According  to  Equation  3-27,  the  X-force  is  simplified 


to: 

X  -  -  j  p  V*S  Cp  +  T 


(3-145) 
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We  can  assume  that  neither  V  or  T  depends  on  the 
control  deflection.  Thus  differentiating  Equation  3-145, 


the  elevator  mainly,  in  producing  a  change  in  the  lift  force 
acting  on  the  tailplane.  A  positive  elavator  deflection 
produces  an  increase  of  lift,  i.e.  negative  Z-force. 

It  has  a  very  small  value  on  aircrafts  with  tail  but 
a  large  value  on  tailless  aircrafts. 

3.  The  Derivative 

Is  the  rate  of  change  in  the  M-pitching  moment  due  to 
a  control  surface  perturbation  divided  by  the  Iyy  moment  of 


inertia 


_L  w 

1  yy  35 


(3-155) 
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According  to  Equation  3-3,  the  pitching  moment  is 


expressed  as: 


M  =  ~  p  c  Cm 


(3-156) 


Differentiating  Equation  3-156,  with  respect  to  &  we 


obtain: 


The  change  in  Z-force  produced  by  the  elevator 
control  surface  deflection,  multiplied  by  a  rather  long  arm, 
provides  a  pitching  moment  in  the  negative  sense,  i.e.  a 
positive  elevator  deflection  results  in  a  nose-down 
(negative)  pitching  moment. 

That  is  the  primary  function  of  the  elevator  and  is 


considered  to  be  the  most  important  of  all  the  control 


surface  functions.  It  adjusts  the  angle  of  attack  and  it 
trims  the  vehicle  in  its  longitudinal  motion,  by  providing 
opposite  pitching  moments  to  accommodate  disturbing  moments 
and  center  of  gravity  relocations. 

For  a  design  point  of  view,  the  larger  the  center  of 
gravity  range  is,  the  greater  the  elevator  pitching  moments 
must  be  provided,  i.e.  larger  values. 

4.  The  Yj  Derivative 

Is  the  change  in  the  Y-side  force  due  to  a  5  control 
surface  perturbation  divided  by  the  mass,  m. 


Y  -  1  21 

5  '  (3-160) 

Differentiating  Equation  3-75,  with  respect  to  S  we 

obtain: 


21 


% 

7* 


If  we  let  V**0o  and 

C<)S  r 

we  arrive  at  the  expression: 


(3-161) 

(3-162) 


Y5  =  JLHi£ 

2t*i 

The  quantity  Cyj  represents  the  rate  of  change  in  the 
Y-side  force  coefficient  due  to  a  variation  in  a  lateral 
control  surface  deflection. 


Qs 


(3-163) 


146 


This  derivative  rises  mainly  from  the  rudder.  A 
positive  rudder  deflection  is  defined  to  produce  a  positive 
side  force.  Therefore,  the  value  of  this  derivative  is 
positive  but  small  enough  to  be  neglected. 

Rudder  deflections  also  cause  an  increase  of  the 
total  drag  of  the  vehicle  but  that  has  a  negligible  effect  on 
the  longitudinal  motion,  and  is  neglected. 

5.  The  U  Derivative 

Is  the  rate  of  change  in  the  L-rolling  moment  due  to 
a  c)  control  surface  perturbation,  divided  by  the  Ixx  moment 
of  inertia. 


t 


The  side  force  generated  by  a  rudder  deflection, 
produces  a  rolling  moment  in  the  positive  sense,  i.e.  a 
positive  rudder  deflection  generates  a  positive  rolling 
moment.  In  airplanes,  with  large  vertical  tails,  this  effect 
is  negligible  and  is  often  neglected. 

The  main  contribution  to  this  derivative  is  from  the 
ailerons,  and  it  is  a  measure  of  the  aileron  effectiveness. 
A  positive  aileron  deflection  is  defined  as  the  one  to 
produce  a  positive  rolling  moment,  i.e.  left  aileron  down  and 
right  aileron  up,  therefore,  the  value  of  this  derivative  is 
positive. 

In  lateral-directional  dynamics,  this  derivative  is 
probably  the  most  important  because  it  establishes  the 
maximum  available  rate  of  roll  of  an  airplane  which  is  a  very 
important  consideration  in  combat  environment. 

It  also  plays  an  important  role  in  take-offs,  and 
landings,  where  adequate  and  rapid  lateral  control  is  highly 
desirable. 

Rapid  rate  of  change  of  aileron  deflections  gives 
al3o  a  rise  to  this  derivative  due  to  lag  effects,  i.e.  not 
immediate  adjustment  of  the  airstream. 

6.  The  Ng  Derivative 

Is  the  rate  of  change  in  the  N-yawing  moment  due  to  a 
S  control  surface  perturbation,  divided  by  the  Izz  moment  of 
inertia. 
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Differentiating  Equation  3-77,  wit i  respect  to  S  ,  we 


obtain: 


ON 


*5-Pv‘S‘  ir 


If  we  let  V=Uo 


Cn§  -  dC-n  /  95 

we  arrive  at  the  expression: 


(3-169) 

(3-170) 


v.  .  _£££i_r 


(3-171) 


The  quantity  Cnj  represents  the  change  in  the  C» 
yawing  moment  coefficient  due  to  a  variation  in  a  lateral 
control  surface  deflection. 

The  side  force  generated  by  a  rudder  deflection, 
produces  a  yawing  moment  in  the  negative  sense,  i.e.  a 
positive  or  to  the  left  rudder  deflection,  produces  a 
negative  yawing  moment.  This  derivative  is  a  measure  of  the 
rudder  effectiveness. 


Rapid  rate  of  change  of  rudder  deflections  gives  also 
a  rise  to  this  derivative  due  to  lag  effects. 

This  derivative  plays  an  important  role  in  crosswind 
takeoffs,  and  landings,  and  counteracts  adverse  yaw  in 
rolling  maneuvers. 


Finally,  aileron  deflections  give  also  a  rise  to  this 
derivative  since  each  aileron  introduces  a  drag  force  of 
different  magnitude.  This  effect  may  be  positive  or  negative 
depanding  on  the  rigging  of  the  aileron  and  the  angle  of 
attack  of  the  airframe. 

3.5  DETERMINATION  OF  THE  STABILITY  DERIVATIVES  AND  OTHER 
CONSIDERATIONS 

All  the  stability  derivatives  discussed  so  far  rise  from 
the  flight  vehicle's  geometry,  being  highly  dependent  upon 
the  design  layout.  We  have  referred  mostly  in  this  chapter 
to  the  conventional  airplane  with  horizontal  tail  and  we  have 
seen  that  the  derivatives  are  made  up  of  a  number  of  separate 
contributions. 

We  had  strictly  remained  in  the  dimensional  form  of  the 
stability  derivatives  and  we  simply  expressed  them  in  terms 
of  airstream  characteristics,  airframe  geometry  and 
non-dimensional  derivatives. 

Dimensional  form  has  the  advantage  of  obtaining  real  time 
solutions  rather  than  working  with  normalized  time  mass  and 
inertia  components  that  must  also  be  non-dimensional  in 
considering  that  form. 

The  advantage  of  non-dimensional  derivatives,  i.e. 
independent  of  size  and  airstream  characteristics  is 
maintained  because  dimensional  derivatives  are  expressed  with 
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the  corresponding  non-dimensionals,  some  of  which  though  do 
depend  upon  the  airstream  and  airframe  characteristics. 

Thus,  some  of  the  dimensional  expressions  and  some  of  the 
non-dimensionals  are  constants,  but  they  can  be  expressed  in 
polynomial  form. 

As  an  example,  for  the  METEOR  7: 

Qv  =  -0.052  -  0.009  CL 
c(r  =  0.0295  +  0.159  CL 
Ctflp  =  “0 .042  CL 
Cvj,  =  -0.182  -  0.052  Ci 

For  others,  plots  may  exist  indicating  their  variation 
with  respect  to  Q  or  mach  number. 

During  the  development  of  an  aircraft,  it  is  important  to 
know  the  numerical  values  of  the  derivatives  under  various 
conditions  of  flight.  In  the  first  stage,  theoretical 
estimates  can  be  made  supported  by  wind  tunnel  tests  and 
model  experiments. 

By  then,  simulation  in  digital  computers  can  be  made  to 
obtain  response  characteristics  to  initial  disturbances  or 
control  deflections  and  determine  whether  the  values  of  the 
derivatives  are  satisfactory. 

When  the  aircraft  has  been  built,  flight  tests  are 
carried  out  to  confirm  the  theoretical  with  the  actual 
responses. 
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There  also  exist  flight  test  techniques  for  determining 
the  values  of  the  derivatives. 

Tables  3.1  and  3.2  summarize  the  longitudinal  and  lateral 
directional  derivatives  including  definition,  unit, 
expression  and  auxilliary  relations. 
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TABLE  3-2 


LATERAL  DIRECTIONAL  STABILITY  DERIVATIVES 
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CHAPTER  4 


LONGITUDINAL  DYNAMICS 


4.1  INTRODUCTION 

Having  specified  all  the  vehicle's  dynamic  properties,  we 
will  now  examine  the  dynamic  motion,  i.e.  the  vehicle's 
response  to  disturbances  and  control  inputs. 

Since  any  flight  vehicle  is  a  dynamic  system,  described 
by  differential  equations  of  time,  we  must  expect  in  a 
particular  situation,  oscillatory  or  exponential  converging 
or  diverging  response. 

Any  particular  response  will  be  associated  by  its 
corresponding  characteristic  root  location  in  the  complex 
plane,  as  we  have  clearly  seen  in  the  dynamic  response  of  the 
spring  mass  damper  system. 

The  characteristic  polynomials  in  the  Laplace  variable 
(s)  involved,  are  usually  of  the  third  degree  or  higher  and 
thus  hard  to  solve  by  hand. 

In  the  past,  a  lot  of  approximations  were  introduced  and 
even  reduction  of  airframe's  degrees  of  freedom  technique  was 
required  to  analyze  coarsely  a  particular  mode  of  motion,  by 
reducing  in  effect  the  degree  of  the  characteristic 
polynomial. 

In  this  chapter,  we  will  analyze  the  longitudinal  motion 
of  a  conventional  airplane.  We  will  first  discuss  the 
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Laplace  Transform  solution  by  obtaining  the  roots  of  the 
associated  characteristic  polynomial.  Then  we  will  apply  the 
reduction  of  degrees  of  freedom  technique  to  isolate  and 
analyze  a  particular  mode  of  motion.  Finally,  the  equations 
of  motion  will  be  rearranged  in  a  form  that  is  more  suitable 
for  computer  solution-.. 

As  an  example,  the  longitudinal  dynamic  motion  of  the 
aircraft  shown  in  Figure  4-1  will  be  analyzed  with  the 
parameters  indicated. 

4.2  CHARACTERISTIC  POLYNOMIAL  OF  LONGITUDINAL  MOTION 

Longitudinal  equations  of  motion  are  farther  simplified 
by  neglecting  the  effect  of  the  Zq  and  Zw  derivatives. 
Experimental  results  have  shown  that  the  effect  of  neglecting 
them  is  only  of  secondary  importance.  Further  a  straight 
level  flight  is  assumed  as  the  reference  flight  condition 
(i.e.  >0)  ,  hence.  Equation  2-152  can  be  reduced  to: 
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The  characteristic  polynomial  is  found  by  evaluation  of 
the  above  determinant.'  The  Equation  4-1  can  be  reduced  in  3 
equations  of  the  original  w,w,q,  longitudinal  variables  by 
avoiding  to  write  down  the  definition  equation  for  q,  but 
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that  doesn't  simplify  the  hand  calculations.  The  determinant 
comes  out  to  be : 


As4  +  Bs3  +-Cs2  +  T)s-t-£  =o  (4-2) 

where 

A-i  (4-3) 

B *  ~  - U»Mi  -  Zw  -Xu  ( 4—4 ) 

C  =  Zm  -  14  Mw  -  Xw  Zu  +■  Xo  C  +  14  MW  *■  Zw)  ( 4—5 ) 

*  "  Xu  (  Z*  -14M*»)  +  Zu(XwM^  f<j  MW)  -  My  (I4Xi*  -^)  (4—6) 

E  *  |  (ZuMw  -  MuZ#)  (4-7) 


Equation  4-2  is  referred  to  as  the  longitudinal  stability 
quartic.  The  roots  in  this  case  for  any  conventional  vehicle 
come  out  to  be  complex  with  negative  real  parts,  which  are 
associated  with  two  oscillatory  convergent  motions. 

The  overall  longitudinal  motion  of  the  vehicle  is  a 
superposition  or  a  combination  of  these  two  oscillatory 
modes. 

One  of  the  modes  usually  is  associated  with  a  pair  of 
complex  roots  located  near  the  origin  of  the  negative  complex 
plane.  This  mode  is  called  "Phugoid  motion"  and  is  a  fairly 
lightly  damped  oscillation  with  long  periodic  time. 

The  other  mode  usually  is  associated  with  a  pair  of 
complex  roots  located  far  from  the  origin  of  the  negative 
complex  plane.  This  mode  is  called  "Short  period  motion"  and 
is  a  heavily  damped  oscillation  with  very  short  periodic 
time. 
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The  longitudinal  oscillatory  modes  are  shown  in 
Figure  4-2. 


For  the  given  example  we  find 


A  -  1.0 

(4-8) 

B  *  4.2177 

(4-9) 

C  »  18.2962 

(4-10) 

D  -  0.1814 

(4-11) 

E  -  0.0722 

(4-12) 

and  the  characteristic 

polynomial  is 

S4  +  4.2177  ss  *■  18-  2962  Sf+  0.(814  s  f  0.0722  =0 

which  has  roots 

(4-13) 

-0.0045  ±  j 

0.0627 

(4-14) 

-2.1043  ±  j 

3.7184 

(4-15) 

15  9 


The  pair  of  roots  -0.0045  ±  j  0.0627  correspond  to  the 
phugoid  motion  and  the  pair  of  roots  -2.1043  ±  j  3.7184 
correspond  to  the  short  motion. 

The  fourth  degree  polynomial  can  be  factored  in  two 
second  degree  polynomials.  Bach  one  of  the  second  degree 
polynomials  is  associated  with  the  corresponding  mode. 

Further  it  is  more  convenient  to  put  each  one  of  the 


polynomials  in  the  familiar  form 


s2+2Su)r,*>  +(Uv?  (4-16) 

so  that  the  natural  frequency  ,  the  damping  5  are  figured 
out  by  inspection  while  the  damped  frequency  can  be 


determined  by: 


uy  =  co*  C  4  -  5* ) 


(4-17) 


If  we  denote  with  the  subscript  (SP)  the  short  period 
quantities  and  with  the  subscript  (PH)  the  phugoid 
quantities,  the  characteristic  polynomial  may  be  written  in 


the  form: 


(s*+  25s,  U)nsp  S  +0^  )(s2t  S  tCOnw)  =0 


(4-18) 


Those  quantities  and  their  relation  is  shown  in  the 
complex  plane  in  Figure  4-3. 

The  roots  of  the  characteristic  polynomial  are  clearly 


shown  in  Figure  4-3  to  be 


-  5S,  (Ons,  ±  j 


(4-19) 

(4-20) 


or 

'  5pu  U>n m  i  j  U>nw  (  ( "  po )  ( 4-21 ) 

~5s»u>r)SP  *  ju>«sp  C  4-5*sp)'te  (4-22} 


Figure  4-3 

Longitudinal  Roots  and  Characteristic  Quantities 
in  the  Complex  Plane 

Other  quantities  of  importance  in  the  longitudinal 

l 

|  dynamics  are:  The  phugoid  period  denoted  by  Tvu  ,  the  short 

|  period  period  TSP  given  by: 

I  TP)I  =  2n/uJdf#4  (sec}  (4-23) 

Tsp  =  2n/(4<sp  Csec)  (4-24) 

and  the  time  for  half  amplitude  for  phugoid  and  short  period 
[  given  by: 

f 


> 
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ivfp(4  = 


0.693/ 

5w»  tl)BM 


(4-25) 


,  0.69  5i 

’  SsrUnsr  (4-26) 

The  Equations  4-25  and  4-26  come  from  the  fact  that  the 


oscillatory  response  decays  exponentially  described  by  the 
-8c  t 

envelope  e  where  Re  is  the  real  part  of  the 

characteristic  root.  Theretore,  to  obtain  half  of  the 
amplitude  the  relation 

=  0.5  (4-27) 

must  exist.  By  taking  the  natural  logarithm  of  both  sides  we 
obtain  the  above  relations. 

For  our  example,  the  characteristic  polynomial  of 
Equation  4-13  can  be  factored  in  the  form 

'sz  +  0.0090  s  +  0.00395)  (s2  +  4.2086s  +  18.2546)  (4-28) 

From  this,  by  inspection  we  can  determine: 


Wnw  *  (0.00395)</2 

=  0.0628  rad/sec 

(4-29) 

WnSP  *  (18.2546)"* 

=  4.2725  rad/sec 

(4-30) 

0.0090 

—  H  D7\7 

3 »  0.0625 

—  \J.\J  l  1  / 

(4-31) 

r  4.2086 

bsp  9*4.2725 

=  0.4925 

(4-32) 

and  also 


0.6951 

0.0717  *  0.06 28 


(55.95  sec 


(4-33) 


,  o.  G95i 

\<l,  —  - —  -  0.55 

*SP  0.4925  *  4  2725 

SCC 

(4-34) 

*  0. 0628(1-0. 07172)^ 

»  0.0627  rad/sec 

(4-35) 

»  4. 2725(1-0. 4925a)'/2 

*  3.7184  rad/sec 

(4-36) 

T”“  *  0.0627  =  W'21  « 

(4-37) 

Ta>  *  3.7184  *  scc 

(4-38) 

It  is  worthwhile  to  note  that: 


(a)  The  damped  frequency  is  indicated  by  the 
imaginary  part  of  the  corresponding  characteristic  root. 

(b)  The  real  part  of  the  characteristic  root 
indicates  the  exponentially  decaying  envelope  of  the 
corresponding  oscillatory  motion,  i.e.  the  damping. 

(c)  The  following  relations  hold 


3 

(4-39) 

C 

=  (Urip^  +  A  5^  CUnSP  +  ^nsp 

(4-40) 

'D 

=  2  Uitf  Jf  +  2S«.  COrtjp  u)n  ^ 

(4-41) 

E 

PM  ^SP 

(4-42) 

The  two  types  of  dynamical  longitudinal  oscillations  are 
shown  in  Figure  4-4. 

The  forward  velocity  v  and  the  pitch  angle  ©  are  usually 
excited  in  the  phugoid  mode  while  the  w  velocity  and,  hence, 
the  angle  of  attack  *  remains  nearly  unchanged. 
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Phugoid  IsngltwdiiMi  oscillation 


ftrt«d  cttgk  a(  dUir  wr*rk>» 

Short- period  bngitnd»M<  octfllaHo* 


Figure  4-4 

Dynamic  Longitudinal  Oscillations 

In  contrast,  the  angle  of  attack  is  heavily  excited  in 
the  short  period  mode  as  well  as  the  pitch  angle  © ,  while  the 
•forward  velocity  remains  nearly  unchanged  in  that  mode. 

The  pilot  can  generally  control  the  phugoid  oscillation 
although  because  of  his  slow  reaction  it  is  possible  to 
worsen  to  unstable  due  to  out  of  phase  action  of  his  control. 
Phugoid  oscillation  must  be  highly  considered  in  take-orfs 
and  landings.  Short  period  is  usually  out  of  pilot  control, 
but  is  of  importance  in  automatic  control. 

The  requirements  that  an  aircraft  should  satisfy  are 
written  in  appropriate  manuals.  It  is  usually  imposed  that 
the  oscillation  should  not  take  more  them  one  cycle  to  damp 
to  half  amplitude. 


164 


Flight  test  can  be  carried  out  to  determine  the 
oscillations  characteristics  under  different  conditions  of 
flight.  Gyroscopic  instruments  can  be  placed  and  recordings 
of  the  longitudinal  variables  can  be  made  possible  to  obtain. 
The  aircraft  is  trimmed  to  the  selected  speed  and  height  and 
longitudinal  oscillation  is  initiated  by  am  abrupt  elevator 
deflection. 

4.3  PHUGOID  AND  SHORT  PERIOD  APPROXIMATIONS 

The  fact  that  the  angle  of  attack  and  hence,  the 
w-velocity  remains  nearly  unchanged  in  the  phugoid  mode,  gave 
the  idea  of  determining  the  phugoid  quadratic  by: 

(a)  eliminating  the  w  terms. 

(b)  eliminating  the  pitching  moment  equation  since  no 
considerable  pitching  moment  variation  occurs,  and 

(c)  letting  q=s9,  obtaining  thus  the  determinant: 
s-X0  9 

-Zu  -Oos  (4-43) 


The  characteristic  phugoid  quadratic  comes  out  to  be: 


s2  ’ Xu  5  ‘  =  ° 


(4-44) 


For  example,  given,  the  characteristic  quadratic  is 
2 


s  +0. 0097s  +  0.0047  =  0 
which  has  roots 

-0.0049  ±  j  0.0681 


(4-45) 


(4-46) 
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These  roots  are  approximately  equal  to  the  exact  phugoid 
roots. 

From  Equation  4-44  it  is  clear  that  the  natural  frequency 
of  the  phugoid  is  heavily  dependent  upon  the  Zu  derivative 
while  the  damping  of  the  phugoid  is  dependent  upon  the  Xu 
derivative.  It  is  also  obvious  that  high  speeds  tend  to 
reduce  the  phugoid  frequency  and  to  increase  the 
corresponding  period. 

Since  for  subsonic  straight  level  flight  the  Zu 
derivative  is  simply  expressed  as: 

Zu  -  -  2 C|  /  Ui>  ( 4—47 ) 

the  following  approximate  relations  hold  for  the  phugoid 


>» 

(4-48) 

Do  Xu 

3  Vi  £ 

(4-49) 

Aircraft  designers,  in  order  to  achieve  larger  values  of 
phugoid  damping,  must  increase  the  value  of  the  Xu  derivative 
but  this  derivative  is  proportional  to  the  drag. 
Fortunately,  in  take-offs  and  landings  where  phugoid  motion 
is  undesired,  high  drag  devices  such  as  flaps  are  used 
anyway . 

A  more  precise  approximation  for  the  phugoid  quadratic  is 
given  by  leaving  the  w  terms  and  the  Hu  term  which  in  modern 
crafts  is  never  zero.  This  is  a  three  degree  of  freedom 
phugoid  approximation. 
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S-^u 

x« 

-Zu 

S-  Zw 

-  Mo 

-M« 

(4-50) 


The  charactecistic  phugoid  quadratic  becomes  then 

2  [  v  Mu(UoXw-j)  1  3  (  n  Mu  ,  \ 


(4-51) 


which  gives 


^  ”=  f  tt  ^  lv  ' 

r  Mu  (UoXw  -g)  1 

^  „  i-  u  DaMw  J 

"  '  afifV.it  7  V2 

2l  u, 


(4-52) 


(4-53) 


Short  period  approximations  can  be  made  possible  by: 

(a)  eliminating  the  v-equations  since  very  small 
v-velocity  change  occurs. 

(b)  setting  all  X-force  derivatives  to  zero  since 
very  small  X-force  variation  occurs. 

(c)  letting  q*s6,  obtaining  thus  the  determinant: 


s-Zw  -Do 

-sMw-Mvm  s-M^  (4-54) 

The  characteristic  short  period  quadratic  comes  out  to 


+  (-UoMw  -Mf-Zwjs  +  (-U0Mw  +  MjZw)  =  0  (4-55) 
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For  the  example  given,  the  characteristic  quadratic  is: 

s2  +  4.208s  +  18.2556  =  0  (4-56) 

which  has  roots 

-2.104  ±  j  3.7187  (4-57) 

These  roots  are  approximately  equal  to  the  exact  short 
period  roots. 

In  summary,  for  typical  flight  conditions,  the  phugoid 
mode  can  be  considered  to  consist  of  changes  in  the 
u-velocity  and  the  9-pitch  angle,  while  the  short  period  mode 
can  be  considered  to  consist  of  changes  in  the  w-velocity 
(angle  of  attack)  and  the  pitch  rate  q. 

For  the  short  period  mode,  the  following  approximate 
relations  hold: 

(JUtJjp  =  CM^Zw  -U>M*)4'2  (4-58) 

y  -(U0Mw+M<t+Zw) 

5sP  =  2  CM^Z*  (4”59) 

A  critical  condition  occurs  when  the  center  of  gravity  is 

in  that  position  for  which 

-  U>Mw  (4-60) 

For  this  case  the  characteristic  equation  becomes: 

s  £  s  *  (-  UoMw  -  -  ZM  )J  =o  ( 4-61 ) 

which  has  roots 

0  j  U0M  w  +  +  Zvu  (4—62) 

The  non-zero  root  is  usually  negative  which  yields  in  a 

non-oscillatory  and  heavily  damped  exponential  response. 
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4.4  SOLUTION  TO  THE  LONGITUDINAL  EQUATIONS 


The  longitudinal  equations  form  a  set  of  linear  first 
order  differential  equations  with  constant  coefficients. 
They  are  accompanied  by  a  set  of  equal  numbered  initial 
conditions,  i.e.  initial  deviations  from  the  equilibrium 
flight  conditions,  and/or  a  number  of  control  inputs. 

Many  ways  have  been  developed  for  solving  sets  of 
differential  equations.  Laplace  Transform  method  is  usually 
the  simplest  although  it  is  not  suitable  for  computer 
solutions. 

In  the  development  of  the  solution,  we  would  like  to 
rearrange  the  equations  in  the  matrix  form: 

x(-t)  =  A  y  (t)  +  B  u[t)  (4-63) 

where  *(t)  is  a  n-dimensional  time  varying  column  vector 
named  state  vector,  ji(t)  is  a  r-dimensional  time  varying 
column  vector  named  control  input  vector,  A  is  a  nxn 
constant  matrix  and  a  is  a  nxr  constant  matrix. 

Taking  the  Laplace  Transform  of  Equation  4-63  according 
to  theorem  3  of  Appendix  A,  we  obtain: 

s  X(s)  ^  *(*)  +  A  X  6)  Ufs) 
or 

(sI-A)X(s)  *  x(o)  +  U(s) 

or 

X(s)  *  (sI-A)  x(o)  4-  (si  -  A)  B  U  (s) 
where  1  is  the  nxn  identity  matrix  and  *(0)  is  the 
initial  condition  column  vector. 
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(4-64) 

(4-65) 

(4-66) 


The  first  term  of  the  transformed  solution  corresponds  to 
the  initial  conditions  and  is  the  homogeneous  solution  and 
the  second  term  corresponds  to  the  control  input  and  is  the 


particular  solution. 

By  inverse  Laplace  Transforming  Equation  4-66  we  obtain 
the  time  solution: 

X(t)  =[Z1CsI-Ay1  ]  x(o)  +  i’1  [(sI-AY^Ucs)]  (4-67) 

To  implement  this  solution,  in  a  digital  computer,  time 
domain  formulation  is  followed. 

The  key  is  to  multiply  the  matrix  differential 
Equation  4-63  by  exp(-dt)  obtaining  thus: 

£-K  [  ift)  -  A  i  M]  =  e  il  ?  u«) 


or 


(4-69) 


Integrating  between  0  and  t.  Equation  4-69  gives 


(4-70) 


or 


Art-r) 

Z  d  u(t)dr 


(4-71) 


The  exponential  term  is  expanded  in  Taylor  series  as: 

e-1  ♦  jr  Kakl  *•  -«■£  Akik  (4-72) 

In  reality,  a  series  truncated  to  ten  or  twenty  terms  is 
often  quite  sufficient  for  normal  accuracy  requirements. 
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The  computer  program  determines  the  values  of  the  state 
variables  at  time  intervals  selected  by  the  user. 

We  now  go  back  to  the  untransformed  longitudinal 
equations  and  try  to  arrange  them  in  the  form  of  the  matrix 
Equation  4-63 . 

For  convenience,  the  unt ransf orraed  equations  are 


replicated  below: 

V  -  (Xu)u -CXw)w  f|©  =(X$)S  (4-73) 

-i.Zo)v  +  */  -(Zw)tAI  -( U>)^  =  (Zs)S  (4-74) 

-(Mu)u  -(Mw)w  +  (£  S  CM\s)S  (4-75) 

®  *  %  (4-76) 


Equation  4-73  is  written: 

u  -  (Xu)u  ^-(k'w)w  -(|)0 -+Cys)§ 
Equation  4-74  is  written: 

\!j  r  (2u)v  +(2w)w  fCUj^  +Czs)& 
Equation  4-75  is  written: 

r  (Mu)  u  +  (Mw)u>  +  (Mw)w  +0^5)5 

or  substituting  Equation  4-7  8  for  we  obtain: 

^  =  (Mu  u  +■  C Mw  f  M»i» Zw)  +  ( MwU>)^  +‘CMj+  MwZs)5 

or  if  we  define: 

Mu  =  Mu  +  Mi  Zu 
M  m  -  M«  Mw  Zw 
*  MwU, 

M5*  a  Ms  Mw  Zb 


(4-77) 

(4-78) 

(4-79) 

(4-80) 

(4-81) 

(4-82) 

(4-83) 

(4-84) 


we  will  have: 

i 


CMu)u  +CMZ)m  ♦  CMS)  5 


(4-85) 
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It  is  further  desirable  to  include  the  height  of  the 
vehicle  as  another  variable  of  interest.  The  equation  for 
the  height  is  obtained  from  the  original  Z-equation 
Equation  2-134  which  for  g.  =0  becomes 

-  Z  -  «  -U>4  (4-86) 

4W  v 

The  left-hand  side  of  this  equation  is  just  the  z-axis 
acceleration,  i.e. 

-ctz  =  -  H  «•  vu  -  U,  ^  (4-87) 

Integrating  the  last  equation  with  respect  to  time  we 
obtain 

k  =  -  W  +Uo9  (4-88) 

the  constant  of  integration  being  the  initial  w-velocity 
perturbation  which  is  zero.  Equation  4-88  describes  the 
height  perturbation  around  the  nominal  value. 

Collecting  Equations  4-77,  4-78,  4-85  and  4-88  in  matrix 
form  we  obtain  the  so  called  state  variable  representation  of 
the  longitudinal  dynamics  model,  the  state  variables  being, 
u,w,q,0,  and  h. 


V 

*1>  X*  0  0 

u 

w 

Zu  Zw  U#  o  0 

nv 

2 

% 

= 

o 

o 

h’ 

% 

0 

o 

o 

o 

o 

0 

0 

It 

o 

I 

o 

O 

k 

0 
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State  variable  representation  form  is  perfectly  suitable 
for  computer  solution.  Many  ordinary  differential  equation 
solver  software  routines  are  available  for  solution.  Most  of 
these  are  provided  with  plotting  capabilities  of  the  state 
variables. 

The  solution  that  follows  corresponds  to  the  given 
example  with  initial  conditions  and  control  inputs  indicated 
in  the  summary  of  specifications  tables. 

It  can  be  seen  from  the  plots  that  the  v-velocity,  the 
0-pitch  angle  and  the  h-height  are  heavily  excited  in  the 
phugoid  motion,  while  the  w-velocity  and  the  q-pitch  rate  are 
heavily  excited  in  the  short  period  motion. 

For  the  phugoid  approximation  the  following  model  may  be 
used: 

it  Xu  O  u  X5 

9  =  - Zulu i  o  o  ■  9  +  Z5/Uo  S  (4-90) 

h  0  Vo  o  K  o 

while  for  the  short  period  approximation: 


w  Zy  U>  *>  is 

£ 

^  ”  Mj  ^  +  Nfa 

4.5  LONGITUDINAL  AERODYNAMIC  TRANSFER  FUNCTIONS 


(4-91) 


Up  to  now  we  went  through  the  analysis  of  the  dynamical 
behavior  of  the  flight  vehicle,  by  making  use  of  the 
mathematical  model  that  adequately  characterizes  the 
longitudinal  characteristic  motions.  In  any  dynamical  system 


m 


two  types  of  variables  exist,  the  input  and  the  output 
variables.  The  input  variables  influence  the  outp.ut 
variables  in  a  manner  determined  by  the  dynamics  of  the 
system. 

The  objective  of  this  section  is  that  of  determining  the 
nature  of  those  input  -  output  relations  called  transfer 
functions. 

For  the  longitudinal  model,  the  input  variable  was  the 
elevator  control  deflection  &  and  perhaps  the  initial 
conditions  or  disturbances  that  can  be  treated  as  inputs, 
while  the  output  variables  were  v,w,q, 0  and  h. 

The  corresponding  transfer  functions  are  obtained  by 
solving  simultaneously  at  a  time,  the  transformed  equations 
of  motion  for  the  output  variable  of  interest,  keeping  all 
initial  conditions  zero. 

The  transfer  function  0(s)/5(s)  is 


s-X« 

0 

*5 

-Zo 

S~2w  -U« 

2i 

-  Mw 

-sMw-M*  s-M^ 

MS 

60s)  _ 

0 

o  I 

0 

Kj) 

As4 

+■  Cj*  E 

where 


As  S1  *•  fig  5  4-  Cfl 

(4-92> 


A*  *  M*  *  MwZs  (4-93) 

'Bg  =  *S(Z0M -M>(Xv*Zw)  (4-94) 

Ce  =  X8(Z.,Mw-ZwMu)+ZS(MwXw-M*Xv)f  MS(z»/X«/ -X*Zk)  (4-95) 
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Ae  (s  +  Z0,)Cs+  z««) 
As,**  Cs2  +Vs  ♦•£ 


(4-96) 


where  29k  ,  z»,  are  the  roots  of  the  numerator  polynomial 
called  zeroes  of  the  transfer  function. 

Similar ily  we  can  derive 


_  S0EJ 

Vs)  "  fa 


S  Cs+Zg,)(s+Zg,) 

As^Ss3*Cs2+:k  +  E 


(4-97) 


where 


wCs)  Aw  s3+-  Bw  s1  +Cw s  +  !D») 

scs)  '  ’7^7iI*7cI23Tri_ 


Aw  -  Z§ 


'Bw  *  /jZu  -ZsC^u  +  M^  +  USU, 

c*>  =r  XS  Cu„  Mu  -Zu  )  *  ZS  Xw  4^  -Ck«3  Xu 

3w  =  ^  (  ZSMU  -M5ZU) 


(4-98) 


(4-99) 

(4-100) 

(4-101) 


(4-102) 


This  transfer  function  may  take  the  following  forms 


Aw  Cs+zw,)(s1^25wc0ws+tui) 
As4  +  BsJ  •*■  Cs2+3s  +  E 


(4-103) 


*i(i)  Aw  Cs+2w<)  C  S+Zio})  f  S+Zu»t) 

S(s)  As4  f  Bsi+-Es24-Ps +  e 


(4-104) 


We  can  therefore  have  three  zeroes,  two  complex  and  one 
real  or  three  real. 

Other  longitudinal  transfer  functions  are: 


Au  +  Bu  S1  *■  Cu  S  +■  4)„ 
A^+Bs^Cs^Ds^E 


(4-105) 
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where 


Ay  =  XS 

%u  -  -  C  Zv  +  +  UoMvg )  +  Z 5  Xw 

Cy  =  XS  C  Zw  -  Do  M*i)  -  Zi  (yw  f^Mi)  4.  MS  (0oXw-£) 

Vv  (M&Zw  -2s  Mw) 

This  transfer  function  may  also  take  the  forms 

IXS)  Ay  CS4-Zy,)  Cs24-2?uCUuS-I-CUu1) 

SfeT  =  A^+BsJ+CsHDSfE 


(4-106) 

(4-107) 

(4-108) 

(4-109) 


(4-110) 


Au  Cs+ZujCs+ZyJCs+Zyj) 

As^Bs^CsVk  +  E 


Similarily 


Kb)  AkS3  f  3k  &2 »  4$  4-Pk 
to)  “  s(  As'  +  '&Us'aPs+e) 

where 

Ah  =  -7; 

3)k  =  -  X5 Zy  +  Zs  (  M(^  (k  M  i»  *-  Xy ) 

Cw  *  VSZo  CM^aMw)-26(ky(M^U,Mi»)-aMw]  -  MSU.M* 

3k  s  *  ^sU»(2w  My  -MwZy)  ■*■  Zj  [  My  (U«Xw-£  )  -&|U"X<;] 

+  Mj  [  U,  Zw  Xu  -  Zy  CU,Xw-^)J 

and  may  take  for  forms 

h(s)  _  4Cs»OCs**25htifrs  fcuh2 ) 

^  ~  s  CAs4  f'Bs’i-Cs^Ps^E) 


(4-111) 


(4-112) 


(4-113) 

(4-114) 

(4-115) 


(4-116) 


(4-117) 


kb)  .  Ak  Cs+Zk,)(s*4,X5+2Kt) 

SA*1  5  (As4*  Bs*+  Cs2  ) 


(4-118) 
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Another  transfer  function  of  interest  is  that  relating 
the  z-axis  acceleration  to  the  control  input,  i.e. 


az(s)  _  -s2kv  Ak  s  (s*zh()  C  s1*  25h^Hs  ) 

'  5(s)  =  (4-119) 

or 

q2(s)  _  Jih  s  Cs*-zhl)Cuzhl)Cs*zi,l) 

S(s>  "  y\s4+'Bs^Csa+Ds^-& 

(4-120) 

Those  transfer  functions  are  used  in  automatic  or  manual 
flight  control  systems,  since  they  specify  the  input  -  output 
relationships  of  the  variables.  Por  example,  to  find  the 
open  loop  response  ©(t)  due  to  a  unit  step  elevator  input  we 
use  the  transfer  function  Equation  4-96 


SCO .  [  4£1  Ks)  =  ie) 

^  S(s)  J  -^s4+3«^+Csl*l>sf'E 


(4-121) 


Taking  the  inverse  Laplace  Transform  of  Equation  4-121  we 
obtain  the  time  response. 


4.6  LONGITUDINAL  EQUATIONS  IN  NON-DIMENSIONAL  SYSTEMS 

The  equations  of  motion  of  a  flight  vehicle,  can  be 
written  in  a  number  of  different  forms,  depending  on  the 
particular  axes  system  selected  for  the  definition  of 
variables  and  the  stability  analysis. 

We  have  just  described  the  body  axes  system  in  the 
dimensional  form,  in  which  most  stability  work  is  carried  out 
[Ref.  1] . 
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In  most  aerodynamic  problems  it  is  more  convenient  to  use 
non-dimensional  coefficients  to  represent  the  parameters 
involved.  By  this  non-dimensional  system,  the  equations  of 
motion  can  be  described  in  the  body  axes  or  wind  axes  system 
as  well. 

The  non-dimensional  coefficients  have  the  advantage  that 
the  effect  of  speed,  size  and  air  density  are  automatically 
accounted.  However,  it  is  required  to  express  the  mass  and 
the  time  in  non-dimensional  forms  as  well.  Therefore,  one 
cannot  obtain  real  time  solutions. 

In  the  non-dimensional  body  axes  system  the  longitudinal 
equations  of  motion  are  given  as  reference  [Ref.  2], 

2 p s  -  Cx„  -Cx*  Cu  u  o 

3Q„-CZu  2(js-s  Cz„ -r-i«  -  2|js  •  #  -  CZn  -n  (4-122) 

~Cmu  -S  Cm^  -  Grig,  (.,  S2  -  s  9  s  Gn^  Cw„ 

where  the  non-dimensional  mass,  inertia  and  time  are 
expressed  as: 

| /  =  ™/pSe  (4-123) 

2  being  the  chavacteristic  length  taken  as  the  half  of  the 
aerodynamic  chord 


In  the  non-dimensional  wind  axes  system  the  longitudinal 
equations  are  [Ref.  3] 


jCCv-Q)  jCt 

V 

0 

C<-  s  +  jCl«  -s 

* 

*-2rt 

®  s  Cm*  *  Gw*  s  Ciri£  -  s*Jgi 

9 

where  the  non-dimensional  mass,  inertia  and  time  are 
expressed  as: 


p  =  w  /p  S  i 


(4-127) 


r  =  w/p  SUo 


(4-128) 
(4-129) 
(4-13  0) 
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TABLE  4-1 


SPECIFICATIONS  FOR  COMPUTER  SOLUTION  OF  THE 
LONGITUDINAL  EQUATIONS  DEDICATED  FOR  PHUGOID 
MOTION  DUE  TO  INITIAL  CONDITIONS 


A.  Variables  and  Initial  Conditions 
u  =  5.0  ft/sec 
w  ■  2.5  ft/sec 
q  =  0.05  rad/sec 
0  *  0.075  rad 
h  =  10.0  ft 
t  *  0.0 


B.  Stability  Derivatives  and  Constants 


Xo 

=  -0.0097 

Zo 

=  -0.0955 

Mu  =  0.0 

g 

=  32.174 

Xw 

=  0.0016 

Z  wv 

*  -1.43 

Mw  =  -0.0235 

Uu 

-  660.0 

x4 

«  0.0 

Zs 

=  -69.8 

Mi,  =  -0.0013 

Mt  =  -1.92 
M4  =  -26.10 


C.  Special  Functions 

My  *  Mu  +■  Z u 

M  J  »  Mw  ♦  Mi  Zjh 
+  Mw 0. 

c  Mj  *  Ui»  Zj 


D.  Control  Input 
§  =  0.0 


E.  Derivatives 

u  =  XyD  +  -|9 

to  »  2u  w  ♦  Z«y  m  +■ +  25  S 

r  Mo  U+  M*U  W  f-  h  5 

®  *  % 

(i  s  -w  +•  Uo0 

F.  Outputs 

17,  w,  q,  9  ,  n  vs  time  at  interval  0.0625 
end  calculation  when  time  >  250.0 
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o 

o 


U  VS  TIME 


X-SCflLE=5. 00E+01  UNITS  INCH. 

Y-SCHLE=2 . QGE+Q i  UNITS  INCH. 

Figure  4-5 

^-Velocity  vs  Time  for  Specifications  of  Table  4-1 
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X-SCRLE=5. 00E+01  UNITS  INCH. 
T-SCRLE=2. GOE  +  OO  UNITS  INCH. 


Figure  4-6 

w-Velocity  vs  Time  for  Specifications  of  Table  4-1 


THETA  VS  TIME 


X-SCRLE=5 . 00E+0 1  UNITS  INCH. 
Y-SCRLE=5 . 00E-02  UNITS  INCH. 


Figure  4-8 

® -Pitch  Angle  vs  Time  for  Specifications  of  Table  4-1 
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H  VS  TIME 


X-SCALE=5 . 0QE+0 1  UNITS  INCH. 
T-SCALE=5. 00E+02  UNITS  INCH. 

Figure  4-9 

h-Height  vs  Time  for  Specifications  of  Table  4-1 
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TABLE  4-2 


SPECIFICATIONS  FOR  COMPUTER  SOLUTION  OF  THE 
LONGITUDINAL  EQUATIONS  DEDICATED  FOR  SHORT  PERIOD 
MOTION  DUE  TO  INITIAL  CONDITIONS 


A.  Variables  and  Initial  Conditions 
v  =  5 .0 

w  *  2 .5 
q  =*  0.05 
6  ■  0.075 
h  =  10.0 
t  »  0  .0 

B.  Stability  Derivatives  and  Constants 


Xu  a  -0.0097 

Zu  a  -0.0955 

Mu 

=  0.0 

g  =  32.174 

X-  =  0.0016 

Zw  *  -1.43 

Mw 

a  -0.02355 

va  a  66  0.0 

Xj  a  0  . 0 

00 

• 

a\ 

vo 

1 

ll 

•o 

Mo, 

a  -0.0013 

n\ 

a  -1.92 

Ms 

a  -26.10 

Special  Functions 

D. 

Control  Input 

%  »  0.0 

Mu  *  Mu  +  MwZu 
Mm  n  Mw  *  Mil  Zmt 
4*  M  ir  IX 

Ml  •  Ma  <■  Mm  z§ 


E.  Derivatives 

i>  =  Xy  u  >  x»>  «*•  -^0  X5  5 
«»  -  ZuU  +  Zl*>  +U»(^fZs§ 
cj^  r  Mu  u  *  «/  *  ^  +  ill  5 

h  *  -ul  +  U.0 

0  l 

F.  Outputs 

uf  w,  q,0,  h  vs  time  at  interval  0.00125 
end  calculation  when  time  >5.0 
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X-SCPLE=5 . OQE-O 1  UNITS  INCH. 

Y-SCRLE=2. OOE+OO  UNITS  INCH. 

Figure  4-10 

D-Velocity  vs  Time  for  Specifications  of  Table  4-2 


W  VS  TIME 


X-SCRLE=5 . OOE-O 1  UNITS  INCH. 
Y-SCRLE=2. QQE+OO  UNITS  INCH. 

Figure  4-11 

w-Velocity  vs  Time  for  Specifications  of  Table  4-2 
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a  VS  TIME 


X-3CfllE=5 . OOE-O 1  UNITS  INCH. 

Y-5CflLE=2.  00E-Q2  UNITS  INCH. 

Figure  4-12 

q-Pitch  Rate  vs  Time  for  Specifications  of  Table  4-2 


i  on 


THETH  VS  TIM 


X -SCRLE=5 . 00E-01  UNITS  INCH. 

Y-SCRLE=2 . Q0E-G2  UNITS  INCH. 

Figure  4-13 

0 -Pitch  Angle  vs  Time  for  Specifications  of  Table  4-2 


190 


'-V  1 »  ••  «s  . 


H  VS  T I  Ml 


X-SCRLE 
T  -SCALE 


h-Height 


5.00E-01  UNITS  INCH. 

5 . 00E+0 1  UNITS  INCH. 

Figure  4-14 

vs  Time  for  Specifications  of  Table  4-2 


TABLE  4-3 


SPECIFICATIONS  FOR  COMPUTER  SOLUTION  OF  THE 
LONGITUDINAL  EQUATIONS  DUE  TO  INITIAL 
CONDITIONS  AND  CONTROL  INPUTS 


A.  Variables  and  Initial  Conditions 
v  *  5.0 
w  *  2 .5 
q  =*  0.05 
9  »  0.075 
h  »  10.0 
t  »  0.0 


g  *  32.174 

Va  *  66  0.0 


B.  Stability  Derivatives  and  Constants 

Xu  ■  -0.0097  Zu  3  -0.0955  Mu  3  0.0 

X*  «  0.0016  Zm  -  -1.43 

Xs  *  0.0  Z5  “  “69.8 


M* 

Mi 

M. 


M* 


-0.0235 
-0.0013 
-1.92 
-26  .10 


C.  Special  Functions 


Mu  -  Ho  *  Mi  Zu 

Hw  .  M w  +■  Mi  Zmi 

-  M3  *  Z5 


D.  Control  Input 

S  =  0.01  for  0<t<10 .0 

-0.01  for  10.0<t <20.0 
5  »  0.0  for  t>20 .0 


E.  Derivatives 


u  -  Xu  u  *  Xw  w  -^0  + 

i*>  s  "Zu  u  ♦  Z«  **  +  Do  y  ±  23  5 

^  *■  Wj  v  ^  +•  MS  S 

it  »  w  + 

9  '  4 


F.  Outputs 


q ,  g  ,  h  vs  time  at  interval  0.0625 
end  calculation  when  time  >  25  0.0 
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i 


i 


'  -  :  ^  -L.  l=  2  .  juE^O  i  oN!~:  I  :i-h. 

Figure  4-15 

v-Velocity  vs  Time  for  Specifications  of  Table  4-3 
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w-Velocity  vs  Time  for  Specifications  of  Table  4-3 


-Pitch  Rate  vs  Time  for  Specifications  of  Table  4-3 
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-  til. 


Figure  4-19 

h-Height  vs  Time  for  Specifications  of  Table  4-3 


CHAPTER  5 


LATERAL  DIRECTIONAL  DYNAMICS 

5.1  INTRODUCTION 

In  this  chapter  we  will  analyze  the  lateral  directional 
motion  of  a  conventional  airplance,  by  folliwng  the  same 
treatment  as  in  the  preceding  chapter  on  longitudinal 
dynamics. 

In  the  lateral  directional  case  also,  any  particular 
response  will  be  associated  as  in  the  longitudinal  case,  by 
its  corresponding  characteristic  root  location  in  the  complex 
plane. 

Further  simplified  set  of  equations  will  be  developed 
that  apply  to  the  individual  modes  of  motion.  Finally, 
computer  solutions  will  be  emphasized  as  well. 

As  an  example,  the  lateral -directional  motion  of  the  same 
aircraft  as  in  the  preceding  chapter  will  be  analyzed  with 
the  lateral-directional  parameters  indicated  below: 
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aUilude  Ml)  20  000 

weight  (lb)  30.500 

wath  'number  0.638 

air  speed  Cf*A«)  060 


Yv  *  -0.0829  Lp  = 

Y5  >  7.656  Ip  = 

Lr  = 

4  « 

^5  * 


4.77 

-  3  SS 

■4.695 

Np  -  -0.0025 

0.4776 

Nr  =  -00967 

27.25 

N$  -  -QSI5 

0.G66 

N*  *  -4.385 

Ixz/r*»  =  00665 


Ixz/I„  »  0.0370 


5.2  CHARACTERISTIC  POLYNOMIAL  OF  LATERAL-DIRECTIONAL 
MOTION 

Lateral-directional  equations  of  motion  are  further 
simplified  by  neglecting  the  effect  of  Y*  ,  Yp  ,  Yr  ,  L*,  and  N* 
derivatives  which  are  of  no  importance  in  lateral  directional 
dynamics.  Furthermoref  a  straight  level  flight  is  assumed  as 
the  reference  flight  condition  (i.e.  =0)  .  Hence, 

Equation  2-153  is  reduced  to; 


*-Yy 

0 

1 

-9/lfc 

0 

/* 

it 

S-Lp 

0 

o 

p 

4 

-Np 

S-N, 

0 

0 

r 

MS 

0 

1 

0 

-s 

0 

* 

0 

0 

0 

l 

0 

-s 

0 
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The  characteristic  polynomial,  as  known,  is  found  by 
evaluation  of  the  above  determinant  which  comes  out  to  be: 


,Csi  fDs  +  E  *0 

where 

(5-2) 

(5-3) 

(5-4) 

c  =  Np  +  Lp(Yv  +  Nr)  *  Mp  Yy  -Lrj  +  Yv(y^-Lt  *Nrj  r 

^  (5-5) 

-  -NpLp  ♦  Yv  ( KlpLr  -If  Mr)  +  Nplf  -  £  (t*  +  ^  Np)  (5_6) 

E  *  £  Cu  Nr  -  M^Lr)  (5-7) 

Equation  5-2  is  referred  to  as  the  lateral -directional 

stability  quartic.  The  roots  in  this  case  for  any 

conventional  vehicle  are  two  real  negatives  or  positives 
which  are  associated  with  two  convergent  or  divergent  motions 
and  one  pair  of  complex  roots  associated  with  an  oscillatory 
motion. 

The  overall  lateral  directional  motion  of  the  vehicle  is 
a  superposition  or  a  combination  of  all  three  modes. 

The  complex  pair  of  roots  correspond  to  the  as  named 
Dutch  Roll  oscillation,  while  the  large  real  root  corresponds 
to  the  Roll  Subsidence  mode  and  the  small  root  to  the  Spiral 
mode. 
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For  the  numerical  example  given,  we  find: 


A  = 

0.9975 

(5-8) 

B  - 

1.867 

(5-9) 

C  « 

3.6841 

(5-10) 

D  * 

6.2637 

(5-11) 

E  * 

-0.0086 

(5-12) 

and  the  characteristic  polynomial  is: 

0*5975  s4  +  -1867  s1  +  3  634 J  s2  *  6.2657  s  -  O.OO&S  -o 

(5-13) 

which  has  roots 

• 

-0.0465  ±  jl.8784 

(5-14) 

+0.014 

(5-15) 

-1.7801 

(5-16) 

The  pair  of  roots,  -0.0465  ±  jl.8784,  correspond  to  the 
Dutch  Roll  mode,  while  the  root  +0.0014  corresponds  to  the 
Spiral  mode  and  the  root  -1.7801  to  the  Roll  Subsidence  mode. 

Further,  it  is  more  convenient  to  put  the  characteristic 
polynomial  in  the  form: 

$  +cunBt)(s-  ir*)(s-rw')  (5-17) 

so  that  the  natural  frequency,  the  damping  and  the  real  roots 
are  obtained  by  inspection  while  the  damped  frequency  of  the 
Dutch  Roll  can  be  determined  by: 

(<-$«)"■  l5-181 

The  root  plotting  in  the  complex  plane  will  appear  as  in 

Figure  5-1. 
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Figure  5-1 

Lateral -Directional  Roots  and  Characteristic 
Quantities  in  the  Complex  Plane 

The  imaginary  roots  of  the  characteristic  polynomial  can 
be  written: 

t  j  (5-19) 

or 

(5-20) 

Other  quantities  of  importance  in  the  lateral -directional 
dynamics  are:  the  Dutch  Roll  period  denoted  by  given  by: 

~  2n/w<^  (5—21) 

the  half  amplitude  for  the  Dutch  Roll  given  by: 


* 


(5-22) 
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the  half  amplitude  for  the  Roll  Subsidence  mode  (if  stable) 


*s 


0  6931 


(5-23) 


and  the  half  amplitude  for  the  Spiral  mode  (if  stable) 


0.693) 


For  our  example,  the  characteristic  polynomial  of 


Equation  5-13  can  be  factored  in  the  form: 

U2  +  0095s  +  5.5505)  (  s-  0.00l4)(s  +  i780l)  (5-25) 

From  this  by  inspection  we  can  determine: 

0 *  (  3.5505)  1  =  1879  rnd/sec  (5-26) 

=  0.095/2-  1879  =  0  0247  (5-27) 

lib  #  *  0.G93I/0 0247- 1.879  •=  14.9051  sec  (5-28) 

s  0.6931  /; J50J  -  0  3894  sec  (5-29) 


Many  of  the  basic  ideas  involving  longitudinal  stability, 
also  apply  to  lateral  and  directional  stability.  Lateral  and 
directional  stability  are  interrelated  due  to  cross  coupling 
effects  of  roll  and  yaw  motion,  that  give  rise  to  the  three 
modes  of  motion  discussed  so  far. 

In  the  Spiral  mode,  the  airframe  seems  directionally 
nearly  stable,  with  no  sideslip  but  it  is  excited  in  a  banked 
turn.  The  unstable  mode  is  found  usually  in  large  finned 
airplanes  with  no  dihedral,  so  that  the  side  forces  developed 
tend  to  turn  the  plane  into  the  relative  wind  but  since  the 
outer  wing  travels  faster,  generates  more  lift  and  the  plane 
will  roll  to  still  a  higher  bank  angle.  The  flight  of  the 


203 


airplane  is  a  heavy  banked  turn  of  ever-decreasing  radius. 


i.e.  a  tightening  spiral.  The  stable  mode  is  found  in 
airplanes  with  dihedral  where  the  flight  is  a  slight  banked 
turn  of  gradually  increasing  radius.  The  spiral  mode  is 
illustrated  in  Figure  5-2. 


Figure  5-2 

Spiral  Stable  and  Unstable  Modes 

In  the  Roll  Subsistence  mode  the  airframe  seems  laterally 
nearly  stable  with  no  bank  but  it  is  excited  in  a  sideslip 
turn.  The  unstable  mode  is  present  when  side  forces  generate 
yawing  moments  that  continue  to  increase  the  initial  sideslip 
disturbance.  This  condition  may  continue  until  the  airplane 
is  broadside  to  the  relative  wind.  Airplanes  with  successive 
weathercock  stability,  usually  have  this  mode  stable  where  an 
initial  sideslip  tends  to  die  out.  The  Roll  Subsistence  mode 
is  illustrated  in  Figure  5-3. 


Figure  5-3 

Roll  Subsistence  Stable  and  Unstable  Modes 

Dutch  Roll  is  a  motion  which  has  the  characteristics  of 
the  Spiral  convergent  and  Roll  Subsistence  divergent  mode. 
It  is  found  in  airplanes  with  strong  lateral  stability  but 
weak  directional.  If  a  sideslip  disturbance  occurs  and  the 
plane  is  yawing  in  one  direction,  rolling  moments  are 
developed  in  a  counter  motion,  i.e.  yaw  to  the  left  is 
causing  roll  to  the  right  and  the  plane  wags  its  tail  from 
side  to  side.  Dutch  Roll  mode  is  illustrated  in  Figure  5-4. 

The  damped  frequency  of  the  Dutch  Roll  is  indicated  by 
the  imaginary  part  of  the  complex  root,  while  the  real  part 
indicates  the  damping.  To  reduce  the  effects  of  Dutch  Roll, 
an  increase  in  the  directional  and  a  decrease  in  the  lateral 
stability  is  required. 


Figure  5-4 
Dutch  Roll  Mode 


5.3  DUTCH  ROLL,  SPIRAL,  AND  ROLL  SUBSISTENCE  MODE 
APPROXIMATIONS 


By  analogy  with  the  longitudinal  short  period  and  phugoid 
approximations,  we  can  reduce  the  degrees  of  freedom  to 
obtain  Dutch  Roll  Spiral  and  Roll  Subsistence  mode 
approximations. 

Dutch  Roll  is  a  motion  which  is  characterized  by 
negligible  variations  of  roll  and  sum  of  the  rolling  moments 
nearly  zero.  The  Dutch  Roll  quadratic  is  obtained  from  the 
matrix  Equation  5-1,  by:  (a)  eliminating  the  4>,  and  p 
equations,  and  (b)  setting  all  rolling  moment  derivatives  to 
zero  obtaining  thus,  the  determinant: 


s-Yv 


1 

5-Nr 


(5-30) 


The  characteristic  Dutch  Roll  quadratic  comes  out  to  be: 

sl  -  (Yv+  Mr)  s  +  (  Yv  NJr  f  )  «  O  (5-31) 


2 


For  the  example  given,  the  characteristic  Dutch  Roll 


quadratic 

is 

Sa  +  0.1786  s  ♦  5.5579  =  o 

(5-32) 

which  has 

the 

roots 

-0.0893  ±  jl .8 841 

(5-33) 

These  roots  approximate  the  exact  Dutch  Roll  roots. 

From  Equation  5-31  it  is  clear  that  the  frequency  of  the 
Dutch  Roll  is  heavily  dependent  on  the  weathercock  derivative 
^  and  the  damping  on  the  yaw  damping  derivative  Hr  and  the 
derivative. 

All  the  derivatives  involved  in  the  Dutch  Roll  mode  are 
dependent  mainly  from  the  vertical  tail,  and  secondarily,  on 
the  wings  and  fuselage.  Thus,  proper  care  has  to  be  taken 
for  the  vertical  tail  as  far  as  the  Dutch  Roll 
specifications.  It  is  usually  required  that  Dutch  Roll 
oscillations  should  be  damped  to  half  amplitude  within  one 
cycle. 

A  first  approximation  to  the  Spiral  and  Roll  Subsistence 
mode  is  obtained  by  expanding  the  determinant.  Equation  5-30, 
by  the  roll  damping  term  (s-Lp)  or  s(s-Lp)^.  The  determinant 
takes  the  form  of : 

s-Yv  o  « 

0  s(s-Lp)  o  (5-34) 

-ty  0  S-Nr 
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(5-35) 


The  characteristic  equation  comes  out  to  be: 

s(s-lp)  [  S2-  (Yv  +  Mr  )s  +  (  YyNr  +  N/O]  -O 
Here,  the  free  5  corresponds  to  the  Spiral  mode,  the 
(s-Lp)  corresponds  to  the  Roll  Subsistence  mode,  and  the 
quadratic  to  the  Dutch  Roll  mode. 

By  this  approximation,  for  the  given  example,  the  Spiral 
root  is  zero,  the  Roll  Subsistence  root  is  -1.695  and  the 
Dutch  Roll  mode  as  previously  -0.0893  ±  jl.8841. 

A  more  precise  approximation  for  the  Spiral  and  Roll 
Subsistence  modes,  is  obtained  by  expanding  the  determinant 
0  -9/U.  i 

-Lp  s(s-lp)  -Lr  (5-36) 

-N|t  s-Np  s*Mr 

This  lead  to 


Subsistence  roots 


5.4  SOLUTION  TO  THE  LATERAL-DIRECTIONAL  EQUATIONS 

Following  the  same  approach  as  in  the  preceeding  chapter, 
we  first  try  to  arrange  the  lateral -directional  equations  in 


the  form 

*(t)  -  A  xtt)  *§  u(l)  (5-41) 

For  convenience,  the  untransformed  equations  are 
replicated  below: 

£-(Yv)f  +  r  -  C9/0.)  4>  =  W)S  (5-42) 

-Of)  ft  +p-(Lp)p  -(-fxz/ixy)r  -  (ir)r  =(Ls)5  (5-43) 

-(W^)/3  -  (lxr/lzz)  p  -  (Np)  p  +  r  -  (Nr)r  ?iNs)S  (5-44) 

4»  =  f  (5-45) 

V  =  r  (5-46) 

Equation  5-42  is  written 

fi  =  (Y V)(S  -r  t-(Vu,)4>  +  CY*)$  (5-47) 

Equation  5-43  is  written 

p  =  (Lf j)(S  +  CLp) f>  +  (Ixr/jxx)f  ♦  (Lr)r  *■ (Lf)1!  «-CLs)5  (5-48) 

Equation  5-44  is  written 


r  =  f  CWl„)f  +•  CNp) p  +  CNr)r  +(MS)5  4/5")$  (5-49) 

In  Equations  5-47  ,  5-48  and  5-49,  the  control  inputs  are 
expanded  in  terms  of  aileron  (?)  and  rudder  (?)  deflections, 
as  indicated  by  Equation  3-143  neglecting  the  dotted 
derivatives  L§  ,  JLjj  . 
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Substituting  Equation  5-49  into  Equation  5-48  we  obtain 


P  =  U  *  Tl%»  N#  g  +  Lp  »  IxMxx  Np  +  Lr  *■  Wl,,  Mr 

4-  ^lxviwIa  <-  ^^“/ixxTzZ  4  -  I^/lxx  Izz 

+  ^  »  Iy>/lxx  Mg  C  f  4  +  ^j4r»  Ms  ^ 

(5-50) 

Similarly  substituting  Equation  5-48  into  Equation  5-49 
we  obtain 


f  Huhkl£.  tk^kkf  ,. 

~  XI/lxxIzz  4 ~  '^/ixylzz  4  -  ^  **/:Iy*Jzz 

+  Ms  is  -  +  Ns+^/izzU  ^ 

^-TWIvkI«  ^  +  4-I2«/IxxIzz 

(5-51) 

We  now  define  for  convenience 

,*  4  ‘  r%,  tv 

f  '  <5-52> 


/  =  VHe.Np 
P  <-IVx«i« 


etc. ,  so  that  Equation  5-50  and  5-51  can  be  written 

y  s(tf)fi*CCT)r  +(£)r  Mf)S  +UV*  (5-54) 

r  =  f(^)^CNf)p  *Wj)F  +CNS)5  (5-55) 
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Collecting  Equations  5-47,  5-54,  5-55,  5-45  and  5-46  in 
matrix  form  we  obtain  the  state  variable  representation  of 
the  lateral  directional  dynamics  model,  the  state  variables 
being  p,  r ,  <P  and  V . 


The  solution  that  follows  corresponds  to  the  given 
example  with  initial  conditions  and  control  inputs  indicated 
in  the  summary  of  specifications  tables. 

For  the  Dutch  Roll  approximation  the  following  model  may 
be  used 

Yv  -L  o  ji  0  Yj 

Np  Nr  O  r  +  N|  Ns  ^  (5-57) 

o  i  o  ry  o  o 

which  includes  all  the  parameters  of  the  determinant. 
Equation  5-30,  while  for  the  Roll  Subsistence  and  Spiral 
approximation 


1* 

0 

0 

f 

v 

h 

r 

- 

Y 

Np 

<P 

0 

1 

0 

0 

-l  S/tt  o 

lr  0  0 

Nr  o  o 

0  O  o 

l  o  o 


o  y| 

? 

Jig  is 

r 

+ 

o  o 

♦ 

1  C?  0 

I 

5 


(5-58) 
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which  includes  all  the  paramet  rs  of  the  determinant. 
Equation  5-36 . 


5.5  LATERAL-DIRECTIONAL  AERODYNAMIC  TRANSFER  FUNCTIONS 

Following  a  same  approach  as  in  the  preceeding  chapter, 
we  can  develope  lateral-directional  aerodynamix  transfer 
functions. 

For  the  lateral-directional  model,  the  input  variables 
were  the  aileron  control  deflection  (?)  and  the  rudder 
control  deflection  (5)  while  the  output  variables  were  (5,  p,  r, 
<t>  and  . 

Thus 

Yj  0  f  -j/u  0 

4  s-Lp  -s^-U  o  o 

Hi  S-Nr  o  o 

0  i  0  -s  o 

_  0  0 _ i _ o  -S 

S(5)  s[  As4  +Bss+Csa4.Ds+e3 


sS%|a  s2+qts  +  T>fl 
As *  +  1**  +  Cs*+7>s  *  E 


where 

Ap  *■  y$  O*  kIzz) 

.  -Vf  [  If  <•  N,  ♦  H,  *  lr]  -  %  4  -M  J 

■=  Y'i  C  LpNr-NpLr)  «-Nalp  -L$Np 
'"^S  *  (  MjLr 


(5-59) 

(5-60) 

(5-61) 

(5-62) 

(5-63) 
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This  transfer  function  usually  takes  the  form 
(. l(s )  _ 

SC*)  =  As^Bs’fCsHDj  +  E  (5“64) 

Other  lateral -directional  transfer  functions  of  interest 

are 


8ft  >  As4f  Bs*rCs*  *'Ds*e 

(5-65) 

where 

- 

4  +T"/jxy  Ms 

(5-66) 

^  • 

Ys*  C  Lp  *  *%,  Mu )  -  4  (Mr+Yr)  *  NS  Ur  -  Yr ) 

(5-67) 

^  - 

YT  ClrN(*  -Lp  Mr)  +  U(YvNr  +  Np)-  Ns  CLp  *■  Y*l_r) 

(5-68) 

and  takes  the  form 

fro  -4  (s2  + 25*0^5  +  u>*) 

5(0  As^+Bs^CsMis  +  E 

\ 

(5-69) 

Also 

r(sl  _  Ar  53+  8r  $*+G  S  +Dr 

Sft)  J(sSBss+a^Tb+E 

(5-70) 

where 

Ar  - 

MS  *  X«/Tlz  LS 

(5-71) 

%  - 

Yfi  CN(J  +  I)a/jn  ijt)  +  4(Np-I*%2  Yyj  -  Ms  (Yr+Lp) 

(5-72) 

Cr  = 

Ys  C  4  Mp  -  Mp  Lp)  -LjYrMp  +■  NgYrlp 

(5-73) 

3 

CuN,t  -  M$i(0 

(5-74) 

and  takes  the  following  forms 


ffc)  _  Ar  C S»2r,)  (sa»25rCOrS4-tUr2) 
is4f3sJ+Cs2».3>svE 


(5-75) 


or 

r(£)  Ar(S*-Zr.)(s*Zr1)(sfZ  r,") 

Sfe)  ^s4<-^sJ4-Cs;*-^s*-E  (5-76) 

It  is  further  desirable  to  include  the  transfer  function 
relating  the  acceleration  on  the  y-axis  to  the  control  input. 
The  acceleration  on  the  y-axis  is  given  from  the  original 
V-equation,  Equation  2-135,  which  for  =0  becomes 

dy  =  v  4-U>r-<g<t>  (5-77) 

and  the  corresponding  transfer  function  is: 

dj  Cs)  Aa  S4  +  Ua  33  +  Q,  S2  +  D„  S  4-  Ea 

As44.'Bs*4-Cs1  +  'DsfE  (5-78) 

where 


A«  x  (5-79) 

'fla  -  UF  *  Mr  4  Ixz/jXJtNp  +  V )  (5-80) 

Ca  =  YS  (  N(,  +  Lp  Mr  -Np  Lr  *  If )  -U.Yv  (*%, U  ►  Mj)  <  5-81 ) 

T)«=Ys[Nplf-ii|i  -W/tCLp^4I%),)]^.Yv[Ls(I-Mf)+Ni(^x^)]  (5-82) 

Ea  *  g  [Y^  Of  N,  -tyU)  *  Yv(KlSU-L5Nr)]  (5-83) 

and  takes  the  form 

A®  (  StZa,)  (s+Za,  Ks+Za,)  Cs+Z«^) 

^  ( 5-84 ) 


5.6  LATERAL-DIRECTIONAL  EQUATIONS  IN  NON-DIMENSIONAL 
SYSTEMS 

In  the  non-dimensional  body  axes  system  the 
lateral -directional  equations  of  motion  are  [Ref.  2]: 
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ty*  -C9,» 

S*tA  -S^f  -51,  -(.?r  •  <P  =  |»| 

-C*||  -S1!*  -SG* p  Stc  -C«r  r  C«|  Or»5 

where  the  non-dimensional  mass,  inertia  terms  and  time  are 


(3 

0  c„s 

• 

4) 

= 

Qf  ^5 

r 

C»j  0r>5 

(5-85) 


expressed  as 

|i  r  w/pse  (5-86) 

2  being  now  the  half  of  the  wing  span  (E=  V2 ) 

lA  =  Iww /f  S€*  (5-87) 

tc  =  1*2 /pS«*  (5-88) 

i,  --  Wpse3  (5-89) 

-l’*'  *  */&  (5-90) 

In  the  non-dimensional  wind  axes  system  the  lateral 
directional  equations  are  [Ref.  31 : 


-2s  +  Cfy  2S 

ct  £  0 

i«<fr 

-S‘I«  .  V  ^  .  (S-M) 

p0? 

jCi)r-sXz 

1  sGip  0  Gn? 

where  the  non-dimensional  mass,  inertia  terms  and  time  are 
expressed  as 

fi  -  0T7  /pSb  (5-92) 


lx  *  a 

Iz  -  2 


TV 

Izz/n* 


X  =  rm  /pSU« 

t*  •  Vr 


(5-93) 

(5-94) 

(5-95) 

(5-96) 
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TABLE  5-1 


SPECIFICATIONS  FOR  COMPUTER  SOLUTION  OF  THE  LATERAL- 
DIRECTIONAL  EQUATIONS  DUE  TO  INITIAL  CONDITIONS 


A.  Variables  and  Initial  Conditions 
/3  =  0.05 
P  *  0.05 
r  *  0.075 
i>  =  0.1 
ti>  a  0 .2 
t  a  0 .0 


B. 


Stability  Derviatives  and  Constants 


Yv 

Y5 


-0.0829 
7.656 


I* 

Lp 

Lr 

Lf 

L, 


-4.77 

N/. 

3 

3  .55 

Ixz/Ixx  =  0.0663 

-1.695 

NP 

3 

-0.0025 

Ixz/Izz  =  0.0370 

0.1776 

Nr 

3 

-0.0957 

g  =  32.714 

27  .25 
0.666 

NS 

Ns 

* 

-0.615 

-1.383 

v.  =  66  0.0 

C.  Special  Functions 


Y*  -  V s/U, 

*  <-  (  JWiwIm) 

Lp,  »  [  Lp  +(WlJ^]  /D 

^  *  [  ty*  +  ^T”/l«  )  Lp]  /D 

Lp  =  [  Lp  Np]  ID 

=  [Np  +  0*/jzz)  Lp]  /D 

L*  *  [  Lr  +  CI«/lx>)  Nr]  /D 

M*  *  [  Mr  +n«/l„)  Lr  J  /  J> 

L*5  *  [  L?  *  (  Wl**)  N5]  /D 

*  [Mg  *(Wj«  )  Lgl  /j> 

4  =  [  /5  +  N5]/3 

Nj  *  [  M$  *(Wl«)  L5]/D 

D.  Control  Inputs 
l  »  0.0 
5  a  0  .0 


TABLE  5-1  (Continued) 


E.  Derivatives 

1$  *  yy  (J  -  r  +C  ?/U>)  <$>  ♦  5 

'P  a  (l  +  f  *"  L?  r  +  If  ? 

r  *  ^Npp  ♦  N?r  +  N*  f  +  Ny  5 

4»  =  ? 

T}>  *  r 

F.  Outputs 

/if  Pf  r,  4>,ip  vs  time  at  interval 
end  calculation  when  time  >  50.0 


.0125 
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Figure  5-5 

Sideslip  Angle  vs  Time  for  Specifications  of  Table  5-1 
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Vi 


.  w \  i  ~  3  IN C  n . 

Figure  5-7 

Rate  vs  Time  for  Specifications  of  Table  5-1 
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Figure  5-8 

$-Bank  Angle  vs  Time  for  Specifications  of  Table  5 
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■jJcJ  uu  1  JUcJ  j  o  w  j  L  4 


TABLE  5-2 


SPECIFICATIONS  FOR  COMPUTER  SOLUTION  OF  THE  LATERAL- 
DIRECTIONAL  EQUATIONS  DUE  TO  INITIAL  CONDITIONS  AND 

CONTROL  INPUTS 


A.  Variables  and  Initial  Conditions 

=  0.05 
p  *  0.05 
r  ■  0.075 
9  »  0.1 
tp  »  0.2 
t  -  0.0 

B.  Stability  Derivatives  and  Constants 


Y*  *  -0.0829 

Lm 

X 

-4.77 

N/» 

3 

3.55 

Ixz/Ixx  =  0 

Y5  »  7.656 

Lp 

* 

-1.695 

NP 

* 

-0.0025 

Ixz/Izz  *  0 

U 

s 

0.1776 

Nr 

3 

-0.0957 

y  =  32.174 

Li 

3 

27.25 

N, 

= 

-0.615 

=  66  0.0 

l5 

= 

0.666 

n5 

3 

-1.383 

C.  Special  Functions 


yJ  =  y5  /Uo 


T> 

«  4-C 

I«  /lx* 

$ 

• 

+  C  Wlw)  fy  ]  /3 

¥ 

-'[N(J  «■ 

(Win 

)lf]  /t> 

<? 

=  [Lp 

♦  CWlx»)  NpJ  /p 

n;  . 

:  [  Np  + 

( Ixr  /jzz 

)Lf3/» 

c 

=  [Lr 

t  Cl«/jyx)  Nr]  /!> 

Nr  ■ 

'  [  ^  r 

(Txz/Ja 

3  Lr]  /D 

4 

•  CLi 

+  N?]  /P 

Ng  • 

r  [  N? 

Lf  3  /J> 

4 

3  t  L5 

+  (it*/!**)  JD 

*  L  ^5 

♦  d«/j„)Ls  j/d 

D.  Control  Inputs 

5  *  0.01  for  0<t<2.5 
5  =  -0.01  for  2.5<t<5.0 
5  *  0.0  for  t>5.0 


5  =  0.0  for  0<t<15 .0 
5  =  0.1  for  15 ,0<t<17  .5 
5  =  -0.1  for  17  .5<t<20 .0 
's  =*  0.0  for  t>20 .0 


0663 

0370 
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TABLE  5-2  (Continued) 


E.  Derivatives 

(4  r  Y„  /i  -  r  *  (9/u,)  4»  +  Y5  5 
?  *  If  ft  +  Lp  P  *  Lr  r  +  if  *  <-  i-5  5 
r  =  Nl^f «  v  Np  P  +  Nr  r  <-  N  j  F  <■  5  ? 

F.  Outputs 

fz  p#  r,  <j 1,1)  vs  time  at  interval 
end  calculation  when  time  >  50.0 


4>  *  f 
2f>  =>  r 

5.0125 
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X  —  SCRLE  =  1 . 00E  +  Q1  UNITS  INCH. 

V-SCflLE=5. 00E-02  UNITS  INCH. 

Figure  5-10 

0 -Sideslip  Angle  vs  Time  for  Specifications  of  Table  5-2 


i'' 
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p 


V 5  TIME 


X-SCRLE= 1 . OOE+0 1  UNITS  INCH. 
r-SCRLE=2. 00E-01  UNITS  INCH. 


Figure  5-11 

p-Roll  Rate  vs  Time  for  Specifications  of  Table  5-2 


-002  -001  000  001  002 


X-SCRLE= 1 . CGE+0 1  UNITS  INCH. 

T-SCRLE=1 . 00E-01  UNITS  INCH. 

Figure  5-12 

r-Yaw  Rate  vs  Time  for  Specifications  of  Table  5-2 
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'•  -tu at* 


cn  co 


FI  VS  TIME 

CRLE=1 . 00E+01  UNITS  INCH. 

CRLE=2. 00E-01  UNITS  INCH. 

Figure  5-13 

<i>-Bank  Angle  vs  Time  for  Specifications  of  Table  5-2 
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X-SCALE=1.0OE-O1  UNITS  INCH. 
T-SCflLE=2.00E-01  UNITS  INCH. 


Figure  5-14 

TV -Yaw  Angle  vs  Time  for  Specifications  of  Table  5-2 
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CHAPTER  6 


STATE  VARIABLE  FEEDBACK  LAW 


6.1  INTRODUCTION 

The  automatic-flight  control  functions  of  an  aircraft  are 
provided  by  an  automatic  flight  control  system  which  is 
referred  to  as  autopilot.  Autopilot  statilizes  the  dynamic 
response  of  the  aircraft  providing  the  best  response  to 
disturbances  and  control  inputs.  Today,  almost  all  types  of 
aircrafts,  helicopters  and  missiles  are  equipped  with 
autopilots  with  complication  depending  on  the  particular 
application. 

The  equipment  comprising  the  autopilot  includes  sensors, 
controllers,  estimators,  computer  and  actuators.  To  damp  the 
vehicle  responses  to  disturbances  and  control  inputs,  the 
rates  of  the  angular  motion  in  roll,  pitch  and  yaw  are  sensed 
by  gyroscopic  instruments  and  fed  back  to  the  flight  control 
computer.  The  computer  combines  the  feedback  information 
with  the  desired  pilot  response  and  computes  commands  which 
are  applied  to  the  actuators  that  deflect  the  control 
surfaces . 

In  this  chapter  we  will  discuss  the  state  variable 
feedback  method  as  a  means  to  control  the  vehicle. 


^,1:, _ _ 
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6.2  THE  FEEDBACK  CONCEPT 


The  negative  feedback  attempts  to  modify  the 
characteristic  equation  of  the  system  so  that  one  can  get  the 
specified  by  the  requirements  response. 

Consider  as  an  example  that  the  roll  motion  of  an 
aircraft  is  described  by  the  equation 

j,  =  Lp?  +  L§$  (6-1) 

The  open  loop  transfer  function  is 

_  Ls 

S(s)  *  s-Lp  (6-2) 

Suppose  that  the  p-variable  is  continuously  measured  and 
fed  back  as  shown  in  Figure  6-1. 


Figure  6-1 

Simple  Feedback  System  for  Roll  Stabilization 
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The  closed  loop  transfer  function  is 

—  =  WCs-L  p)  s  IS 

4  +  KFLi  /a-lp)  s+CtpLS-Lp)  (6-3) 

It  is  clear  that  the  single  characteristic  root  changed 
from  Lp  to  Kp-Lg-Lp.  Usually  Lf>  is  negative  and  Ls  positive. 
So  the  system  will  have  a  faster  response  indicated  by  the 
greater  negative  root  of  the  characteristic  polynomial. 

In  automatic  control  systems  all  the  variables  of  the 
system  are  fed  back  with  corresponding  gains  K's  chosen  to 
give  a  particular  desired  response. 

In  general,  the  system  is  described  in  matrix  form  as 


=  A  x  ^  u 


with  initial  condition  vector 


X  (O')  =  X, 

Applying  the  feedback  law 


U  =  -  <  X  4-  U, 


where 


(6-4) 


(6-5) 


(6-6) 


K  *  [K4  K2  K,  ...  K„ ]  (6-7) 

is  the  feedback  gain  matrix  and  the  desired  input,  one 


obtains  the  system 


=  (A-  8k  )  x  +  B 


(6-8) 


The  matrix  A-BK  now  determines  the  ’■.ew  characteristic 
roots  of  the  system.  Note  that  the  characteristic  roots  of 
the  system.  Equation  6-4,  are  obtained  by  the  determinant 


I  -  A 


(6-9) 
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while  for  the  system,  Equation  6-8,  by  the  determinant 


,1  -  (A-BlO  1  =  O 


(6-10) 


Since  the  response  of  the  system  is  determined  by  the 
location  of  the  characteristic  roots,  any  desired  response 
will  specify  the  characteristic  polynomial  of  the  system 
which  if  equated  with  Equation  6-10  the  feedback  gains  are 
obtained. 

Feeding  back  all  the  variables  means  that  somehow  all  of 
them  must  be  measured  or  estimated.  Usually  not  all  of  the 
variables  are  measured  but  the  ones  that  are  not  measured  are 
estimated.  For  example,  the  pitch  angle  can  be  obtained  by 


integration  of  the  q-pitch  rate. 

e  -  r  q 


(6-11) 


Even  unstable  systems  can  be  stabilized  with  state 
variable  feedback  provided  that  the  system  is  controlable. 

Most  times  the  feedback  gains  are  selected  such  that  for 
the  system  to  meet  a  performance  criterion.  In  this  case 
optical  control  methods  are  used  to  select  the  gains. 

In  the  next  pages  a  computer  solution  is  shown  of  the 
given  example  for  the  longitudinal  and  lateral  directional 
case  using  the  state  variable  feedback  concept. 
Specifications  for  computer  solution  are  shown  in  Tables  6-1 
and  6-2. 
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TABLE  6-1 


SPECIFICATIONS  FOR  COMPUTER  SOLUTION  OF  THE  LONGITUDINAL 
EQUATIONS  DUE  TO  INITIAL  CONDITIONS,  CONTROL  INPUTS  AND 
STATE  VARIABLE  FEEDBACK  STABILIZATION 


A.  Variables  and  Initial  Conditions 
v  *  5.0  ft/sec 

w  *  2 .5  ft/sec 
q  *  0.05  rad/sec 
&  »  0.075  rad 
h  *  10.0  ft 
t  »  0.0 

B.  Stability  Deriviatives  and  Constants 


X„  *  -0.0097  Zu  *  -0.0955 

X*  »  0.0016  Z*  *  -1.43 

X&  *  0.0  Zs  »  -69.8 

C.  Special  Functions  D. 

My  ar  My  ■*  MwZu 

Mui  *  *  tJjiLi 

m  *■  14(2,  It 
Mi*  -  Ms  ♦-  Um  Za 

E.  Derivatives 

o  *  yuu  +yu*'«  -^9  +  5 

ii  r  2u  o  +■  Z\u  **  4-  ll>  ZS  S 
^  *  Mtr  0  +  Mw  W  4.  ^  v  M3  S 

e  a  % 
k  *  - 1»  +•  U*9 


Mu  *  0.0  g  *  32.174 

M*  =  -0.0235  v.  =  660.0 
M*  =  -0.0013 
M^  =  -1.92 
Ms  »  -26.10 

Control  Inputs 

=  0.01+  f  for  0<t<10.0 
a  -0.P1+  f  for  10 ,0<t<20 .0 
a  0.0  +  f  for  t>20 .0 

f  *  -0.005  u  -  0.0005  m  ♦  0. 022. 

+  0.005  9  +  0. OOOOOO 1  Vl 


F.  Outputs 

v,  w,  q,  9  ,  h  vs  time  of  interval  0.0625 
end  calculation  when  time  >  250.0 
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X-SCRLE=5 . OOE  +  G 1  UNITS  INCH. 

T-SCRLE=2. OOE+OO  UNITS  INCH. 

Figure  6-2 

v>-Velocity  vs  Time  for  Specifications  of  Table  6-1 
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W  VS  TIME 

X-SCfiLE=5. OQE+Oi  UNITS  INCH. 
Y-SCflLE=2 . OOE+O 1  UNITS  INCH. 

Figure  6-3 

w-Velocity  vs  Time  for  Specifications  of  Table  6-1 
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Q  VS  TIME 


X-SCPLE=5 . 00E+0 I  UNITS  INCH. 

Y-5CPLE=5. 00E-02  UNITS  INCH. 

Figure  6-4 

q-Pitch  Rate  vs  Time  for  Specifications  of  Table  6-1 
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-005  000  005  010 


CO  CO 


CRLE=5 . 00E+0 1  UNITS  INCH. 

CflLE=5. OGE^Ol  UNITS  INCH. 

Figure  6-6 

h-Height  va  Time  for  Specifications  of  Table  6-1 


TABLE  6-2 


SPECIFICATIONS  FOR  COMPUTER  SOLUTION  OF  THE  LATERAL- 
DIRECTIONAL  EQUATIONS  DUE  TO  INITIAL  CONDITIONS,  CONTROL 
INPUTS  AND  STATE  VARIABLE  FEEDBACK  STABILIZATION 


A.  Variables  and  Initial  Conditions 
/S  a  0.05  rad 

p  «  0.05  rad/sec 
r  «  0.075  rad/sec 
fi  *  0.1  rad 
V  a  0 .2  rad 
t  »  0,0  sec 

B.  Stability  Derivatives  and  Constants 

Yy  a  -0.0829  Lp  a  -4.77  Np  »  3.55  Ixz/Ixx  =  0.0663 

Ys  »  7.656  Lp  =  -1.695  NP  a  -0.0025  Ixz/Izz  a  0.0370 

Lr  =  0.1776  Nr  a  -0.0957  q  *  32.174 

L*  a  27.25  N|  a  -0.615  u,  »  660.0 

L5  =  0.666  Ns  »  -1.383 

C.  Special  Functions 


Yj  =  y5/Uo 

D  a  4  -  (  l\z  /  Ixx  Izz  ) 

CJ  *  [  LjS  Wl«0  N|i  ]/v  N?  •  [  Hf  +  CY«/i„  )  in  3  /J> 

fp  =  [  Lp  +CYw/Ixx)  Np]  / d  NjT  a  [  Up  f  Lp3  ID 

L*r  *  [  Lr  *(te/j„)Nr]/V  Nr  -  [  Mr  +Ll«z/lu)  l r] 

-  C  L?  t  C  I*z./jKl()  NgJ/j)  Nf  =•  [  Nf  f  CIkz/jS7  )  Lg]/D 

*  L  L$  «■  CIxt/jxx)  M53/3>  5  L^l/T) 


D.  Control  Inputs 

?  »  0.01  +  |f  for  0<t<2 .5  S»  0.0  +  cf  for  0<t<15.0 

5  «  -0.01  +  5f  for  2 .5 <t<5 .0  T  =  0.1  +  «f  for  15.0<t<17.5 

f  »  0.0  +  If  for  t>5 .0  a  -0.1  +  5 f  for  17.5<t<20.0 

*»  0.0  +  5f  for  t>20 .0 

|f  a  ~O.0<<8  /i  -o.zll  p  -0,209  r  -0  324  <p  -0,  2G67  d> 

^f  *  0.0515  £  +  0.050S  p  «.  0,495  r  ►  O.O^SZ  <p 
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TABLE  6-2  (Continued) 


E.  Derivatives 

(5  =  Y y  jS  t  r  +  (9/y„)  <$>  +  Y5  5 

-P  -  +  +  l^.r 

r  =  £  ♦  Mp  p  +  Mr  r  ♦  |  v  N5  5 

4*  -  ? 

lj>  =■  r 

F.  Outputs 

(*»  Pr  r,  4»  r  ^  vs  time  at  interval  0.0125 
end  calculation  when  time  >  50.0 
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X-SCflLE= 1 . 00E+0 1  UNITS  INCH. 
Y-5CPLE=5. 0CE-02  UNITS  INCH. 


Figure  6-7 

(3 -Sideslip  Angle  vs  Time  for  Specifications  of  Table  6-2 


P  VS  TIME 


X-SCRLE=1.00E+01  UNITS  INCH. 
T-SCRLE=  l.OOE  —  01  UNITS  INCH. 


Figure  6-8 

p-Roll  Rate  vs  Time  for  Specifications  of  Table  6-2 


R 


VS  TIME 


X-SCRLt=l . OOE-Oi  UNITS  inch. 
T-SCRLE=5. 00E-02  UNITS  INCH. 

Figure  6-9 

r-Ya w  Rate  vs  Time  for  Specifications  of  Table  6-2 
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X-SCr)LE  =  1  .  OOE^O  1  UNITS  INCH. 
T-SCPLE=1 . OOE-Oi  UNITS  INCH. 


Figure  6-10 

4>-Bank  Angle  vs  Time  for  Specifications  of  Table  6-2 


X -SC RLE = 1 . GOE^O i  UNITS  INCH. 

T-SCRLE=1 . GCE-Oi  UNITS  INCH. 

Figure  6-11 

V-Yaw  Angle  vs  Time  for  Specifications  of  Table  6-2 
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APTER  7 


SPECIALIZATION  TO  MISSILES 


7.1  INTRODUCTION 

The  preceding  -hapters  have  been  oriented  primarily 
towards  the  dynamics  of  the  conventional  piloted  airplanes. 
The  equations  of  motion  for  the  missile  are  still  valid  since 
it  can  be  thought  as  an  automatically  controlled  aircraft. 

Missiles  are  found  in  different  geometry  considerations. 
Cruise  missiles  usually  have  no  rotational  symmetry,  i.e. 
rotation  through  360"  is  required  along  the  missile  axis  to 
leave  the  external  shape  unchanged.  Configurations  with  180“ 
rotational  symmetry  are  seldom  encountered  while  120°  or 
higher-order  rotational  symmetry  are  more  typical. 

In  this  chapter  we  will  consider  the  cruciform  missile, 
i.e.,  with  90°  rotational  symmetry  to  take  the  advantage  of 
symmetries.  This  configuration  is  shown  in  Figure  7-1.  This 
missile  may  be  considered  as  typical  for  today's  large 
existing  missiles. 

In  the  next  sections  we  will  consider  the  simplification 
to  the  equations  of  motion  due  to  symmetry  and  the 
complication  of  concerning  the  coupling  effect  of 
longitudinal  and  lateral-directional  motions.  The  last 
happens  because  strong  rolling  moments  are  produced  when  the 
missile  is  simultaneously  pitched  and  yawed. 
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Cruciform  Missile  Configuration 


For  this  particular  model  we  will  assume  that  control 
inputs  are  coming  from  the  vertical  and  horizontal  fin 
deflectors  of  the  tail.  A  deflection  of  the  vertical  fins 
stands  for  a  rudder  deflection  while  a  deflection  of  the 
horizontal  fins  stands  for  an  elevator  deflection  .  These 
are  the  only  existing  control  inputs.  Roll  control  is 
obtained  by  differential  deflection  of  one  vertical  or 
horizontal  fin. 


7.2  SYMMETRY  CONSIDERATIONS 

For  the  missile's  equations  of  motion,  the  reference  body 
axes  system  shown  in  Figure  7-1  will  be  used  as  in  the 
airplane  case,  and  straight  level  horizontal  flight  will  be 
assumed  as  the  nominal  flight  condition. 

Considering  the  existing  symmetries,  we  find  first  that 
in  addition  to  the  aircraft  symmetries  the  xy  plane  is  a 
plane  of  symmetry  and  that  makes  the  inertia  tensor  of  the 
missile  having  only  diagonal  terms  and  as  a  consequence  the 
body  axes  coinsides  with  the  principal  axes.  Further 
Iyy  *  Izz. 

The  fact  that  the  missile  is  not  designed  as  an  airplane, 
gives  no  tendency  to  remain  in  the  same  roll  orientation. 
Vertical  and  horizontal  find  deflections  can  be  optimized  to 
give  zero  bank  angle  in  turning  maneuvers.  This  condition  is 
highly  desirable  in  missiles  as  well  as  in  airplanes  because 
the  lifting  forces  are  most  efficiently  generated  and 
lift-to-drag  ratio  is  maximum. 

However,  accidential  actuator  error  signals,  asymmetrical 
loading  of  the  control  and  lifting  surfaces  in  supersonic 
flight  and  atmospheric  disturbances  may  introduce  large  roll 
rates  and  large  roll  rates  cause  cross-coupling  between  the 
longitudinal  and  the  lateral -directional  motion.  This  makes 
the  equations  of  motion  not  separable  in  a  longitudinal  and  a 
lateral-directional  set  which  was  found  so  convenient  in  the 
aircraft  case. 
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We  therefore  find  the  equations  of  motion.  Equation  2-134 
and  2-135,  combined  in  one  set  of  equations  and  modified  as 
follows: 


v  *ge 


(7-1) 


(7-2) 


(7-3) 


ST  Y  *  »  -‘O-r 


(7-4) 


? 


(7-5) 


IiZ 


(7-6) 


with  auxilliary  relations 


?  • 


r  =  V? 


(7-7) 

(7-8) 

(7-9) 


Before  expanding  the  forces  and  moments  in  terms  of 
stability  derivatives  we  note  tht  due  to  symmetry  the  Y-force 
is  the  same  as  the  Z-force  which  makes  v«w  and 

Yv  ■  Zw  (7-10) 

Y;  -  Zw  (7-11) 
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Also  the  N-yawing  moment 
which  makes  r»-q  and 
Yr  ■  -Zq 
Nv  »  -Mw 
Nv  -  -MW 
Nr  «  Mq 

because  Iyy  *  Izz 

In  addition,  for  the 


is  the  same  as  the  M-pitching  moment 

(7-12) 

(7-13) 

(7-14) 

(7-15) 


control  derivatives  one  has 


Y5  »  -zn 


(7-16) 


N5  =  Mn 
N5  *  M* 


(7-17) 

(7-18) 


7.3  EQUATIONS  OF  MOTION  FOR  THE  NON-ROLLING  MISSILE 

If  we  expand  Equations  7-1  to  7-6  in  terms  of  stability 
derivatives  we  obtain 


u  = 

Xu  U  +  Xy,  W 

-^0  +•  Xn  n 

(7-19) 

<ii  = 

Zu  0  +■  2w  w 

+  Ub^  +  2*  n 

(7-20) 

i s 

Mu  u  *  Mw  ui 

Mw  w  +  ^  .Utirt 

(7-21) 

or  substituting  Equation  7-20 

^  MwZu)  u  ♦  (Mw  +  M*Zw)  w  CM^4.  Uw  (Jo)  ^  +(Mn  + jUwZ*)rj  (1-  22) 


*  «  Yu  v 

-  U0r 

♦a*  * 

Y55 

(7-23) 

9  *  Lu  v 

+  Lp  p 

+  Lr  r 

♦  If* 

(7-24) 

r  *  Mv  1/ 

i-  Mv  v 

*  Npp 

+  Mr  r 

4  N5S 

(7-25) 

At  Equation  7- 

-25  w 

e*  would  like 

to  keep 

the  term  N$v 

because 

it  is 

1  the 

same 

in  magnitude 

as  M^S*. 

Substituting 

Equation 

7-23 

for  v 
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(7-26) 

equation  is 


f  =  (Ny^NyYw^  f  NpP  +fMr-NvU>)r  fgMy^  +•  (M2  +  NJvY^)5 
If  the  missile  is  not  rolling  the 
eliminated  as  well  as  the  p-stability  derivatives.  Also,  if 
the  v-velocity  perturbation  is  assumed  zero,  the  equation 
is  eliminated  as  well  as  the  v-stability  derivatives. 
Therefore,  one  obtains  the  following  model  for  the 
non-rolling  missile 


ws  +  f2nw  (7-2?) 

a  (Mw  +  ^wZw)vu  +  +  +  n  +W|j)2n)ri  (7—20) 

V  =  Yy  V  -  Vor  +  Y$$  (7-29) 

r  =■  (  Mv  +•  Yv)v  +  (  Mr 'NvU,  )r  4.  (  N5 My  Y^S  (7-50) 

or  more  conveniently  in  matrix  form 


vu 

Zw 

V. 

0 

0 

w 

Z« 

O 

i 

nC 

M£ 

0 

O 

1, 

4- 

< 

O 

Y 

0 

0 

Yy 

-u. 

V 

0 

Ys 

• 

r 

0 

0 

Nt 

Mr 

y> 

0 

Ms 

(7-31 ) 


This  model  may  be  augmented  with  the  auxiliary  equations 
and  transfer  functions  and  state  variable  feedback  can  be 
applied  for  proper  stabilization. 
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Theorem  1: 
Theorem  2: 
Theorem  3 
Theorem  4: 
Theorem  5 
Theorem  6 


SOME  LAPLACE  TRANSFORMS 


*U)  -t>0 

i{m)f  --  m 

Stt) 

i 

u(4)  Cutlii  slcp) 

V* 

rCt)  or  t  (unil  romp  ) 

Vs2 

42 

Vs3 

V  Csfa) 

V  Csra)2 

cos  m{ 

s/  (  S24-COJ) 

sin  w  4 

"/ (  s^co1) 

3  e at  cos  Cwt  *<p) 

eJ*4+avja>  +  CJVs*a-jco 

SOME  LAPLACE  TRANSFORM  THEOREMS 
*  a  F,C»)  <•  b  F2(s; 

l  [  K  f«)  J  =  < 
l  {  dfCV/dt  |  =  sPfr)-f(<>) 

i  s  Sa  R?)  -  sffo)  -  d^lo 

jt  { civet)/**  j  -  sj  m  -  -  s  dr/di L  -  Midi1 

l  j  fo*  -rctjdt  5  -  f(s)A 
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